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INTRODUCTION 


A system  of  computer-assisted  stress-strain  data  reduction  analysis,  storage, 
and  retrieval  was  described  in  the  two  previous  parts  of  this  three-part  series 
of  reports.1*2  In  this  report,  the  FORTRAN  listings  of  the  computer  programs 
are  presented. 

FORTRAN  symbolic  routines  and  subroutines  of  four  files  catalogued  in  the 
UNIVAC  1106  machine  are  listed.  The  file  names  and  functions  are  as  follows: 

R*PICPTS  for  raw  data  reduction 

R*EV4LR0  for  data  analysis 

R*RETRIEVE  for  accessing  the  data  stored  in  the  DATABANK  file 

R*FIXBANK  for  updating  the  LOCATER  and  DATABANK  files. 

Two  data  files,  R*L0CATER  and  R*DATABANK,  are  not  listed. 

Two  of  the  files,  R*PICPTS  and  R* RETRIEVE,  use  graphics  subroutines  which 
are  part  of  the  TEKTRONIX  PLOT-10  package.  These  subroutines  are  not  listed  and 
neither  are  a group  of  subroutines  in  a BUFFPK  package  which  interfaces  the 
PLOT-10  subroutines  with  the  particular  UNIVAC  1106  executive  system  we  use. 

Both  the  PICPTS  and  the  RETRIEVE  programs  are  interactive,  requiring  entries 
by  the  user  in  response  to  instructions  appearing  at  the  terminal.  Provisions 
have  been  included  in  both  programs  to  detect  a number  of  different  types  of  in- 
put data  errors  and  to  allow  the  user  to  make  immediate  corrections  at  the  terminal. 

One  objective  of  the  developed  computer  procedures  was  to  make  in-house  data 
readily  available  and  hence  more  useful  to  the  Engineering  Standardization  Division 
in  their  continuous  efforts  to  improve  standards  and  tolerances  for  material. 


1.  PAPIRNO,  R.  P.  AMMRC  Automated  Stress-Strain  Data  Analysis,  Storage,  and  Retrieval  ■ Part  I:  Data  Analysis  Program. 
Army  Materials  and  Mechanics  Research  Center,  AMMRC  TR  79-16,  April  1979. 

2.  PAPIRNO,  R.  P.  AMMRC  Automated  Stress-Strain  Data  Analysis,  Storage,  and  Retrieval  - Part  2:  Data  Entry  and  Retrieval 
Procedures.  Army  Materials  and  Mechanics  Research  Center,  AMMRC  TR  79-17,  April  1979. 
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LISTING  OF  R*PICPTS 


The  PICPTS  file  consists  of  the  following  symbolic  elements: 

Element  Function 

MAIN  Main  program  for  reducing  raw  stress-strain  data 

PICK  Subroutine  for  digitizing  a curve  using  the 

TEKTRONIX  4662  interactive  plotter 

BUFFPK  Subroutine  package  for  interfacing  the  TEKTRONIX 

PLOT-10  package  to  the  UNIVAC  1106 

START  Control  routine  to  prepare  for  punching  output  and 

to  start  execution  of  program 

STOP  Control  routine  to  punch  output  at  a remote  UNIVAC 

9300  terminal 

PRINT  Control  routine  to  punch  PICPTS  output,  execute 

main  program  of  EVALRO  and  print  EVALRO  output. 

Punch  and  print  operations  at  a remote  UNIVAC 
9300  terminal. 

MAPI  Control  routine  used  to  PREP  the  file,  form  the 

absolute  (name:  PICKEM),  PACK  the  file,  and 
print  a list  of  file  elements. 

TEKTRONIX  PLOT-10  subroutines  which  are  called  but  for  which  no  listings  are 
given  are  as  follows: 

1.  In  the  element  MAIN:  INITT,  TERM,  CHRSIZ,  ANM0DE,  TOUTST,  M0VABS,  TSEND, 
BELL,  NEWPAG,  DCURSR,  and  FINITT. 

2.  In  the  element  PICK:  TINSTR. 


The  PLOT- 10  subroutines  are  referenced  from  the  file  DA0*TEX  and  from  the 
BUFFPK  package  when  the  absolute  program  is  created. 

PICPTS  is  executed  from  a TEKTRONIX  4014  graphics  terminal  coupled  to  a 
TEKTRONIX  4662  Interactive  Plotter.  User  instructions  and  a comprehensive  de- 
scription of  the  program  output  are  given  in  Reference  3. 

Listings  of  MAIN  are  given  on  pages  3-8  and  PICK  on  page  9.  The  four  con- 
trol routines  are  listed,  in  the  order  given  above,  on  page  9. 


3.  PAPIRNO,  R.  P.  PICPTS  Computer  Program  for  Digitizing  Testing  Machine  Autographic  Records:  User  Instructions. 
Army  Materials  and  Mechanics  Research  Center,  AMMRC  TN  78-6,  June  1978. 
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1 c 

a C •*  THIS  IS  THE  MAIN  PROGRAM  OF  THE  FILE  PICPTS  (ABSOLUTE:  PICKEM) 

3 c 

4 C 

SC*  FILE  ALSO  INCLUDES  SUBROUTINE  PICK  • GRAPHIC  SUBROUTINES  NOT  LISTED 

6 C * GRAPHICS  SUBROUTINES  FROM  TEKTRONIX  PLOT  10  PACKAGE 

7 C 

8 C • INTERACTIVE  PROGRAM  FOR  DATA  REDUCTION  OF  STRESS-STRAIN  TEST  RESULTS 

9 C * EXECUTED  ON  A TEKTRONIX  4014  COUPLED  TO  A TEKTRONIX  4662  DIGITAL  PLOTTER 

10  C 

1 1 C 

12  C 

13  DIMENSION  I P LTON( 3),IPLT0F(3),S( 99 ) > E ( 99 ) , ICHGE ( 99 ) , 

14  1 SS ( 20 ) , EE( 20 ) , 1 1 ( 20 ) , I ICHGE (20 ) 

15  DATA  IPLTON/27,65,69/ 

16  DATA  IPLTOF/27,65,70/ 

17  CALL  I NI TT ( 30 ) 

18  CALL  TERM( 2 . 1 024) 

19  CALL  CHRSI Z ( 3 ) 

20  C 

21  C 

22  C 

23  C * * * READ-IN  OF  IDENTIFYING  DATA  • * * 

24  C 

25  CALL  ANMODE 

26  NRUN=0 

27  NPRINT  *1 

28  C 

30  C 

31  1 PRINT  200 

32  MX- 1 

33  REA0(5,100,ERR«5) , SNO 

34  2 PRINT  201 

35  MX  = 2 

36  READ) 5,101 , ERR-5 ) , PPL , PULT , ARE A , ELONG , RA , RCH , ISET 

37  PRINT  202 

38  READ  102.ANS 

39  IFlANS.EO. 'YES' ) GO  TO  10 

40  3 PRINT  203 

41  MX *3 

42  READ(5,103,ERR«5), WMTL ■ XMT  L , YMT  L ■ ZMTL , KTEMP , KEDOT 

43  4 PRINT  204 

44  MX-4 

45  READ(5, 104,ERR*5) . BOCK , NPACE , NDATE , TYPE 

46  GO  TO  10 

47  C 

49  C 

50  C * * * * * INPUT  ERROR  CORRECTION  ROUTINE  * * * * * 

51  C 

52  5 PRINT  400 

53  READ  102, ANS 

54  I F ( ANS . EQ. 1 YES  1 ) GO  TO  ( 1 , 2 , 3 , 4 , 1 1 , 32 , 33 , 37 ) ,MX 

55  GO  TO  60 

56  C 

57  C 

58  C 

59  C * * * * CALIBRATION  OF  X-Y  PLOTTER  * * • • 

60  C * * IF  PREVIOUS  CALIBRATION  IS  VALID  THIS  ROUTINE  IS  OMITTED  * 

61  C 

62  10  PRINT  205 

63  READ  102, ANS 

64  I F ( ANS . EG . 'YES' ) GO  TO  15 

65  11  PRINT  206 

66  MX-5 

67  READ(5. 101 , ERR-5) , CALS, CALE 

68  XSTOP-CALE/10. 
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69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
1 02 

103 

104 

105 

106 

107 

108 

109 

110 
1 1 1 
112 
113 


YST0P=CALS 
PRINT  207 
KOUNT  *0 

13  KOUNT  =KOUNT+1 
CALL  BELL 

I F ( KOUNT . LE • 3 ) GO  TO  13 
CALL  TOUTST ( 3 , IPLTON) 

CALL  MOVABS (50,50) 

CALL  TSENO 
CALL  PICK(IE.IS) 

I CALOE- I E 
ICALOS=IS 
LOE  = 1 E/4 
LOS= I S/4 

CALL  MOVABS ( LOE , LOS ) 

CALL  TSENG 
CALL  PICK(IE.IS) 

EUNIT=CALE/( IE-ICALOE) 

SUNI T = CALS/ ( IS-ICALOS) 

CALL  TOUTST (3, IPLTOF) 

KOUNT  = 0 

14  KOUNT =KOUNT+1 
CALL  6ELL 

I F ( KOUNT . LE . 3 ) GO  TO  14 
C 



c 

C » * » * » * DIGITIZING  ROUTINE  **••** 

C * * * * * TWO  STEP  PROCESS  FOR  EACH  POINT  **•**•• 

C * » STEP  1 DIGITIZE  CURVE  ZERO  » * STEP  2 DIGITIZE  CURVE  POINT  * 
C 

15  CALL  NEWPAG 

CALL  DCURSR( IC, IXD, IYD) 

CALL  ANMCDE 
DO  1C  1=1,99 
ICHGEI I ) =0 

IFU.GT.20)  GO  TO  16 
1 1 CHGE ( I ) =0 

16  CONTINUE 
NEW  = 0 
NCHGE*0 
SOLD=0. 

E0LD=0 . 

IOUT=0 

K0UT=0 


114 

115 

116 

117 

118 
1 19 
120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 


K = 0 
1 = 0 
N = 0 

PRINT  208 
K0E=  LOE 
K0S=  LOS 
KE  = L0E+50 
KS  = LOS 
KOUNT  = 0 

17  KOUNT  = K0UN1 +1 
CALL  BELL 

I F ( KOUNT . LE . 3 ) GO  TO  17 
CALL  T0UTST(3, IPLTON) 

18  1=1+1 

19  CALL  :.:OVABS(KOE,KOS) 
CALL  TSENO 
CALL  PICK(IE.IS) 

KOE=  I E 
KOS= 1 S 

CALL  MOVABS ( KE , KS ) 

CALL  TSEND 
CALL  PICK (IE, IS) 

E( I )=EUNIT*( IE-KOE) 
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137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 
140 
1 49 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 
1 78 

179 

180 
181 
182 

183 

184 

185 
166 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 


S ( 1 ) * SUN  I T* ( 1 S-KOS ) 

EOUT  = A BS ( E ( 1 ) ) 

SOUT=ABS(S(I)) 

K0E-K0E/4 

K0S=K0S/4 

I F ( EOUT . LT . XSTOP . AND . SOU  T . GT . YSTOP ) GO  TO  20 
I F ( I . GE  . 50  • AND  • MOD  (1,5)*  EQ . 0 ) GO  TO  20 
I F ( I . GE . 90 ) GO  TO  20 
K E = I E /'  4 
KS=IS/4 

I F ( I .GE.50.AND.K.EQ.0)  GO  TO  20 
IF(K.GT.O)  GO  TO  21 
GO  TO  18 

20  N=I 

IF( EOUT. LT. XSTOP. AND . SOUT. GT .YSTOP ) N-I-1 

21  CALL  10UTST(3, IPLTOF) 

CALL  TSENG 

KOUNT  *0 

22  KOUNT  * KOUNT  + 1 
CALL  BELL 

I F ( KOUNT . LE . 3 ) GO  TO  22 
C 

C 

C********  WRITE-OUT  OF  DATA  *•**•••** 

C 

30  CALL  NEWPAG 

CALL  DCURSR ( IC, IXD, I YO) 

CALL  ANMODE 
PRINT  208 , N 
READ  102, ANS 

I F ( ANS . EO. 1 NON  1 ) GO  TO  32 
PRINT  210 

I F (ANS . EQ, 1 ALL  1 ) PRINT  21 1 , ( I , S( I ) , E ( I ) , ICHGE( I ) . 1-1 ,N) 

I F ( ANS . EO . 1 LAS  1 ) PRINT  2 1 1 , N . S ( N ) , E ( N ) , ICHGE ( N ) 

IF(ANS.EQ. 'COR‘ . AND.NEW.GT.O)  GO  TO  31 
IF(KOUT.GT.O)  PRINT  219 

I F ( KOUT . EQ • 1 ) PRINT  2 1 1 , 1 1 ( KOUT ) , SS ( KOUT ) , EE ( KOUT ) , I ICHGE ( KOUT ) 

I F ( KOUT . GT . 1 ) 

1 PRINT  211,1 II(KI),SS(KI),EE(KI),I ICHGE ( K I ) , KI -1 .KOUT) 

GO  TO  32 

31  PRINT  220, NEW 

PRINT  21 1 , NEW, SOLD, EOLD.NCHGE 
PRINT  221 

PRINT  21 1 .NEW.S(NEW) ,E(NEW) .ICHGE(NEW) 

C 



c 

C * 7 OPTIONS:  PUNCH*  ADD*  CORRECT*  DELETE*  RECALIBRATE-  REDIGITIZE*  DISCARD* 
C 

32  PRINT  212 
MX*6 

READ! 5,101 , ERR-5) , IANS 
IF( IANS.EQ. 1 ) GO  TO  50 
1 F ( I ANS . EQ . 2 ) GO  TO  33 
I F ( I ANS . EQ • 4 ) GO  TO  37 
1 F ( IANS . EQ . 5 ) GO  TO  10 
I F ( IANS.EQ. 6)  GO  TO  15 
I F 1 I ANS . EQ . 7 ) CO  TO  60 
C 



C 

C *****  * OPTION  TO  CORRECT  DATA  ••*••• 

C 

K-K+1 

33  PRINT  213 
MX*7 

READ ( 5 , 1 0 1 , ERR-5 ) , I 
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205 

NEW=  I 

206 

NCHGE- I CHGE ( I ) 

207 

SOLO=S( I ) 

208 

EOLD=E( I ) 

209 

I CHGE ( I ) * I CHGE ( I ) + 1 

210 

KE=KOE+£( I )/<4*EUNIT) 

211 

KS-KGS+SI I )/(4»SUNlT) 

212 

KOUNT  = 0 

213 

34 

K0UNT=K0UNT+1 

214 

CALL  BELL 

215 

I F ( KOUNT . LE ■ 3 ) GO  TO  34 

216 

CALL  TOUTST (3, IPLTON) 

217 

GO  TO  19 

218 

C 

220 

C 

221 

C • 

* * * • * option  to  add  data  * 

***** 

222 

C 

223 

35 

CALL  NEWPAG 

224 

CALL  DCURSR(IC.IXD.IYD) 

225 

CALL  ANMODE 

226 

PRINT  214 

227 

KE=K0E+E(N)/(4*EUNIT) 

228 

KS=K0S+S(N)/(4»SUNIT) 

229 

K = 0 

230 

I=N 

231 

KOLiNT  = 0 

232 

36 

K0UNT=K0UNT+1 

233 

CALL  BELL 

234 

I F ( KOLINT  . LC . 3 ) GO  TO  36 

235 

CALL  TOUTST ( 3 , IPLTON) 

236 

GO  TO  18 

237 

C 

< JU 

239 

c 

240 

c * 

* * » * OPTION  TO  DELETE  A POINT  * 

***** 

241 

c 

242 

37 

PRINT  215 

243 

MX  = 8 

244 

RE AD ( 5 , 1 0 1 ,ERR=5) , IOUT 

245 

KOUT  = KOUT  + 1 

246 

SS ( KQUT ) = S( IOUT) 

247 

EE ( KOUT ) =E ( IOUT) 

248 

I I ( KOUT ) = IOUT 

249 

IICHGE(KOUT)*-1 

250 

NX=N- 1 

251 

IF( IOUT.EO.N)  GO  TO  40 

252 

1 F ( IOUT. EQ. NX)  GO  TO  39 

253 

DO  38  I = IOUT , NX 

254 

E( I ) =E( 1+1 ) 

255 

S( I ) =S( 1+1 ) 

256 

I CHGE ( I ) = ICHGE (1+1 ) 

257 

38 

CONTINUE 

258 

GO  TO  40 

259 

39 

E(NX ) =E( N) 

260 

S ( NX ) = S ( N ) 

261 

ICHGE ( NX ) = I CHGE ( N ) 

262 

40 

NsNX 

263 

GO  TO  30 

264 

C 

266 

c 

267 

c • 

* CURVE  IDENTIFIERS  AND  STRESS-STRAIN  DATA 

PUNCHED  ON  CARDS  * 

268 

c 

269 

50 

NRUN=NRUN+1 

270 

IFINRUN.GT. 1 ) GO  TO  51 

271 

PRINT  216 

272 

READ  101 , NCOPY 
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273  PUNCH  302.NPRINT, NCOPY 

274  51  PUNCH  300 , SNO , WMT L , XMT L , YMTL , ZMT L , BOOK , NPAGE , TYPE , NDATE . 

275  1KTEMP.KED0T ,N, PPL, AREA , PULT , E LONG , R A , RCH , I SET 

276  PUNCH  301  , (S( I) ,£( I) , I, I»1 ,N) 

277  C 

279  C 

280  C ♦ * * * • CONTINUE  WITH  ANOTHER  CURVE  OR  STOP  • • * 

281  C 

282  60  CALL  NEWPAG 

283  CALL  DCURSR(IC,IXD,IYD) 

284  CALL  ANMODE 

285  PRINT  217 

286  READ  102, ANS 

287  I F ( ANS . EQ. 1 YES  1 ) GO  TO  1 

288  PRINT  218 

289  CALL  F INI TT ( 0 , 500 ) 

290  C 

292  C 

293  c***»**»  READ  FORMATS  ****•• 

294  C 

295  100  FORMAT ( A6 ) 

296  101  FORMAT ( ) 

297  102  F0RMAT(A3) 

298  103  F0RMAT(2A2.2A4, 15, 12) 

299  104  FORMAT (A4,12,I6,A5) 

300  C 

301  C 

3C2  C 

303  c * * * » » » CRT  WRITE-OUT  FORMATS  «**••** 

304  C 

305  200  FORMAT ( 1 H ,/,T2,  'REFER  TO  PICPTS  DATA  SHEET  FOR  ENTRIES  1 TO  7 ' , 

306  1 // , T2 , 1 IF  SCREEN  ERASES  DURING  PROGRAM  & X-Y  AXES  APPEAR ',/, T2 , 

307  2 ' KEY  SPACE  BAR,  THEN  KEY  RETURN  TO  CONTINUE1, 

308  3//,T2,12( '* 1 ,1X),//,T2, 'ENTRY  1:  SPECIMEN  NUMBER1) 

309  201  F0RMAT(T2, ' ENTRY  2:  PROP  LIMIT  AND  ULT  LOAD  IN  LBS‘,/,T2, 

310  1 1 AREA  IN  S3-IN.  ELONG  AND  REDUCTION  IN  AREA  IN  PCT',/,T2. 

311  2 ' HARDNESS  AND  SET  NO.  * ENTER  -1.  FOR  UNAVAILABLE  DATA'./, 

312  3T2. 1 EXAMPLE:  6000 . , 9780 . . . 05 , 1 0 . 2 , 35 . 6 , 0 . , 22 1 ) 

313  202  FORMAT! T2. ‘ ENTRY  3:  IF  DESCRIPTORS  SAME  AS  LAST  RUN,  ENTER  YES'. 

314  1/.T2,  'OTHERWISE  ENTER  NO') 

315  203  F0RMAT(T2, 'ENTRY  4:  MATERIAL  DESCRIPTORS  ( 2 A2 , 2 A4 ) ' , / , T2 , 

316  1 'TEST  temp  AND  EDOT  EXPONENT  ( 15 , 12  ) ' , / , T2 , 

317  2 ' EXAMPLES  « ST18  52100  20-4') 

318  204  F0RMAT(T2. 'ENTRY  5:  NOTEBOOK,  PAGE,  DATE  « TEST  TYPE',/,T2, 

319  1 ' USE  TENS,  COMP,  OR  SHEAR  FOR  TEST  TYPE  ENTRY', /,T2, 

320  2 ' EXAMPLE:  1 A460207 1 977T ENS ' ) 

321  205  F0RMAT(T2, 'ENTRY  6:  IF  PREVIOUS  CALIBRATION  IS  VALID,  ENTER  YES', 

322  1/.T2, 'OTHERWISE  ENTER  NO') 

323  206  F0RMAT(T2, 'ENTRY  7:  CALIBRATION  VALUES  FOR  STRESS  AND  STRAIN', 

324  1/,T2, 'EXAMPLE  * 6000. ,.02') 

325  207  FORMAT ( T2 , ' * PROCEDURE  FOR  CALIBRATING  PLOTTER  •',/T2, 

326  1 ' 1 - POSITION  THE  CROSS-HAIRS  ON  ZERO  AND  PUSH  CALL',/,T2, 

327  2 ' 2 • WAIT  FOR  CROSS-HAIHS  TO  RETURN  TO  ZERO', /,T2, 

328  3 ' 3 . POSITION  ON  CALIBRATION  POINT  AND  PUSH  CALL’) 

329  208  F0RMAT(T2, '‘PROCEDURE  FOR  DIG  I T I Z I NG* ' , / , T2 , 

330  1 ' EACH  POINT  REQUIRES  TWO  STEPS', /,T2, 

331  2'DIGITIZE  ZERO  THEN  DIGITIZE  CURVE  POINT', /,T2, 

332  3 ' TO  EXIT,  DIGITIZE  ZERO  THEN  DIGITIZE  EXIT  POINT') 

333  209  FORMAT (T2, 'NUMBER  OF  POINTS  DIGITIZED:  N*',I3,/, 

334  1 T2 , ' FOR  A LISTING  OF  ALL  DIGITIZED  DATA  ENTER  ALL',/, 

335  2T2 , ‘ FOR  LAST  POINT  ONLY,  ENTER  LAST',/, 

336  3T2 . ' FOR  PRINT  OF  CORRECTED  POINT  ONLY,  ENTER  CORRECT',/, 

337  4T2 , ' FOR  PRINT  OF  DELETED  POINT(S)  ONLY,  ENTER  DELETE',/. 

338  5T2 , 1 TO  OMIT  ANY  LISTING,  ENTER  NONE.') 

339  210  FORMAT (T2, '**  LIST  OF  DIGITIZED  VALUES  •*',/, T2, 

340  1 ' *NUM  CHGS*  number  OF  TIMES  GIVEN  POINT  HAS  BEEN  CORRECTED',//, 
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341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 

375 

376 

377 

378 

379 

380 

381 

382 

383 
364 

385 

386 

387 

388 

389 

390 

391 
352 
393 


2T2, ' POINT* ,4Xt ‘LOAD* , 3X STRAIN  * 3X, ' NUM 1 , / , T3 , 

3 1 NO . ' ,5X, *LBS'  ,4X, ‘IN/IN* ,4X,  * CHGS * ,/) 

21 1 FORMAT ( T3 , 1 3 , 3X . F7 . 0 , 1 X , F7 . 5 . 3X , 1 2 ) 

212  FORMAT ( 1 H . / , T2  , ' » MAKE  ONE  OF  SEVEN  CHOICES  *’,/,T2, 

1*1.  enter  1 to  PUNCH  DATA  N0W',/,T2, 

2*2.  ENTER  2 TO  ADD  MORE  DATA  POINTS './, T2 , 

3 ' 3 . ENTER  3 TO  CORRECT  ONE  OR  MORE  POI NTS  1 , / , T2 , 

4 1 4 . ENTER  4 TO  DELETE  ONE  OR  MORE  PO I NTS ' , / , T2 , 

5 1 5 . ENTER  5 TO  RECALIBRATE  & REDI GI t I ZE ' , / , T2 , 

6 ' 6 . ENTER  6 TO  REDIGITIZE  THIS  CURVE', /,T2, 

7*7.  ENTER  7 TO  DISCARD  THIS  SET  OF  DATA') 

213  FORMAT ( T2 ,' ENTER  POINT  NO.  TO  BE  CORRECT  ED T2 , 

1 'WAIT  FOR  CROSS-HAIRS  TO  MOVE  TO  ZERO', /,T2, 

2 1 USE  TWO-STEP  PROCEDURE  TO  CORRECT  ONE  POINT*. /T2, 

3 * AUTOMATIC  EXIT  FROM  PLOTTER  AFTER  ONE-POINT  CORRECTION') 

214  F0RMAT(1H  , / , T2 , ' * USE  TWO-STEP  PROCEDURE  TO  ADD  DATA  POINTS 
1/.T2. 'CROSS-HAIRS  WILL  MOVE  TO  ZERO  & THEN  TO  LAST  POINT* , 

2/, T2,  'RESUME  DIGI T I ZING ' , / , T2 , 

3 ' USE  TWO-STEP  PROCEDURE  FOR  EXIT  POINT1) 

215  FORMAT ( T2, * ENTER  POINT  NO.  TO  BE  DE LETED * , / , T2 . 

1 1 DELETED  POINT  WILL  BE  LISTED  WITH  NUM  CHGS--1  ' ) 

216  FORMAT ( T4 , ' I F A PRINTED  PI  REPORT  IS  TO  BE  PREP ARED ' , / , T2 , 

1 'ENTER  THE. NUMBER  OF  COPIES  REQ I RED ' , / , T2  , 

2' IF  NO  REPORT  REQUIRED.  ENTER  1') 

217  FORMAT ( T2 , ' I F YOU  WISH  TO  DIGITIZE  MORE  CURVES, ENTER  YES'. 

1/,T2,  'OTHERWISE  ENTER  NO') 

218  F0RMAT(T2, 'PLEASE  SIGN  OFF  TERMI NAL ' , / , T2 , 

1 ' 1 . IF  ONLY  CARDS  REQD ■ ENTER  ?ADD  R* P ICPTS . STOP ' , / , T2 . 

2 ' 2 . IF  CARDS  & REPORT  REQD:  ENTER  0ADD  R * P I CPTS . PR INT ' , /T2 , 

3*3.  WAIT  FOR  MESSAGE:  READY', /.T2, 

4' IF  PRINT  REQUESTED,  SENT  BY  T32  WILL  ALSO  APPEAR ',/, T2 , 

5 1 4 . THEN  WHEN  PROMPT  (>)  APPEARS  ENTER  PFIN') 

219  FORMAT!/, T2, ' **  FOLLOWING  DELETED  (ORIGINAL  POINT  NOS.):') 

220  FORMAT ( T2 , ' * * VALUES  FOR  POINT', 13,'  BEFORE  CORRECTION:') 

221  F0RMAT(/,T2, 'CORRECTED  VALUES:') 

C 

C 

C*******  PUNCH  FORMATS  ***••*• 

C 

300  FORMAT (A6,2A2,3A4, I2.A5, 16, 15 , 2 12 , F7 . 0 , F6 . 4 , F7. 0 , 3F4 . 1 ,14) 

301  FORMAT (F7.0,F7.5,63X,I3) 

302  FORMAT( 12, 2X, 12) 

C 

C 

C 

C • * * * * * ERROR  DIAGNOSTIC  FORMAT  * * • * • 

C 

400  FORMAT ( T2 , ' * **  ERROR  IN  DATA  JUST  ENTERED  •**',/, T2, 

1 ' TO  CORRECT,  ENTER:  YES',/,T2, 

2 'TO  TERMINATE  EXECUTION  ENTER:  NO') 

C 

END 
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1 c 

2 SUBROUTINE  PICK(IX.IY) 

3 C 

4 C 

5 C • DIGITIZES  LOAD  STRAIN  CURVES  ON  TEKTRONIX  4662  DIGITAL  PLOTTER 

6 C * FULL  SCALE  PLOTTER  RANGE:  X-4095  UNITS  AND  Y*2731  UNITS 

7 C * PLOTTER  COORDINATES  SCALED  TO  LOAD  AND  STRAIN  IN  MAIN  PROGRAM 

8 C 

10  C 

11  DIMENSION  I N ( 7 ) 

12  CALL  TINSTR(7,IN) 

13  IX»0 

14  IVO  5 

15  DO  200  I«1 ,5,2 

16  IDX*MOD( IN( I ) ,32) 

17  I X* I X*32+IDX 

18  IDY=MOD( INI  1+1 ) ,32) 

IS  I Y« I Y*  32+IDY 

20  200  CONTINUE 

21  I X« I X/3 

22  IY-IY/8 

23  RETURN 

24  C 

25  END 


1 9ASG.UP  PUN 

2 9FREE  PUN 

3 9ASG.A  PUN 

4 PBRKPT  PUNCHS/PUN 

5 PXQT  R»PICPTS. PICKEM 


1 PBRKPT  PUNCH* 

2 PFREE  PUN 

3 9SYM.C  PUN , , RM9301 


1 #ASG,A  R«EVALRO. 

2 9BRKPT  PUNCHS 

3 »FREE  PUN 

4 9DA0*UTIL. SUSPEND 

5 9XQT  R*EVALRO. FIT DATA 

6 PADD  PUN. 

7 PDAO*UTIL. RESUME, P RM9301 

8 PSYM.C  PUN..RM9301 


1 PPREP  R»PICPTS. 

2 PMAP.I  CIS. R«PICPTS. PICKEM 

3 IN  R* PICPTS . MAIN 

4 IN  R»PICPTS.8UFFPK 

5 LIB  R»PICPTS. 

6 LIB  OAO*TEX. 

7 END 

8 PPACK  R»P ICPTS . 

9 PPRT.TL  R-PICPTS. 
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LISTING  OF  R*EVALRO 


k 

I 


\ The  EVALRO  file  consists  of  the  following  symbolic  elements: 

L. 

Element  Function 

MAIN  Main  program  for  analysis  of  stress-strain  data 

ORDER  Subroutine  to  arrange  the  input  data  in  increasing 

order  of  strain  magnitude 

LSTSQ  Subroutine  to  perform  a linear  least-squares 

analysis 

RJA75  Subroutine  for  numerical  solution  of  an  equation 

of  the  form:  Z + AXn  + B = 0 

Control  routine  used  to  assign  LOCATER  file  to  run 
and  to  execute  the  program 

Control  routine  used  to  PREP  the  file,  form  the 
absolute  (name:  FITDATA),  PACK  the  file,  and 
print  a list  of  file  elements. 

EVALRO  is  a batch  process  program  executed  from  the  UNIVAC  9300  hard-copy 
terminal  using  punch  card  input  or  from  an  interactive  graphics  terminal  using 
a file  of  punch  card  images.  The  hard-copy  output  on  both  cases  appears  at  the 
UNIVAC  9300  terminal. 

Listings  for  the  main  program  and  the  three  subroutines  are  given  on  pages 
11-25.  The  two  control  routines  are  listed,  in  the  order  given  above,  on  page  26. 


START 

MAPI 
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I 


► 


i_ 


1 c 

2 C * THIS  IS  THE  MAIN  PROGRAM  OF  THE  FILE  R»EVALUO  (ABSOLUTE:  FITDATA) 

3 

4 C 

5 C * THE  FILE  ALSO  CONTAINS  3 SUBROUTINES 

6 C • THE  SUBROUTINES  ARE:  ORDER,  LSTSQ,  AND  RJA75 

7 C * PROGRAM  DOES  BEST-FIT  R-0  ANALYSIS  6 PREPARES  OUTPUT  OATA  FOR  DATA  BANK 

8 C 

9 C' 

10  C 
1 1 
12 

13 

14 

15 

16 

17  C 

18 

19 

20 

21  C 

22 

23 

24 

25  C 

26  C 

27  C 

28  READ  100,NPRINT,NPUNCH,NCOPY 

29  C 

30  C » NPRINT  IS  A PRINT-OUT  OPTION  CODE 

31  C • NPRINT  = 0 • ONLY  BEST-FIT  PARAMETERS  AND  STATISTICS  PRINTEO 

32  C * NPR I NT  = 1 * PI  REPORT  * DATA  FOR  NPRINT. 0 AND  TABLE 

33  C < TABLE  LISTS  INPUT  DATA, CALC  STRESS  & DEVIATION  FOR  BEST-FIT  ONLY 

34  C * NPRINT *2  * ALL  TRIAL  PARAMETERS  AND  STATISTICS  PRINTED 

35  C * NPR I NT=3  * PRINT  OUT  OF  NPRINT. 2 DATA  AND  TABLE:  LISTING  NEXT  LINE 

36  C • TABLES  LIST  INPUT  DATA.  CALC-S,  AND  PCT  DEV  OF  CALC-S  FROM  INPUT-S 

37  C 

3B  C * NPUNCH  Is  punch-out  option  code 

39  C * NPUNCH. 1:  PUNCH  OUT  DATA  ♦ NPUNCH.O:  NO  PUNCH 

40  C 

41  C • NCOPY  INDICATES  NUMBER  OF  COPIES  OF  NPR INT » 1 REQUIRED 

42  C • IF  NO  ENTRY  (NCOPY. 0)  PROGRAM  DEFAULTS  TO  NCOPY-1 

43  C 

44  C 

45  C 

46  C * INITIAL  VALUES  FOR  DATA  BANK  PREPARATION  VARIABLES 

47  C * FILE  10  IS  THE  LOCATER  FILE  FOR  DATABANK 

48  C 

49  I F (NPUNCH , EQ . 0 ) GO  TO  1 

50  DEFINE  FILE  1 0 ( 50 , 80 , E . J J ) 

51  DO  BO  JX« 1 .7 

52  READ( 1 0 1 JX , 502 ) , WA(UX) , ( LOC( JX ,JK) , JK.1 ,19) 

53  BO  CONTINUE 

54  DO  81  JP-1 ,99 

55  JPCH ( UP ) >0 

56  WQUT ( JP ) *0 

57  XOUT ( JP ) «0 

58  I SOU  T ( JP ) .0 

59  81  CONTINUE 

60  C 

61  C 

62  C 

63  1 NRUN.NRUN+1 

64  NPX.NRUN-1 

65  KPRT.O 

66  C 

67  C 

68  C 


DIMENSION  S(99) ,E(99) ,P(99) ,EP(99) ,X(99) ,Y(99) ,CS( 10,99) , 

1 DBAR ( 10,5) ,SBAR( 10,5) , DE L ( 1 0 , 5 ) , CC ( 1 0 , 5 ) , SEE ( 1 0 , 5 ) , SEEM( 1 0 , 5 ) , 
2SY(10) ,SYM( 10) ,XN(10) ,ISPTS(10,5) ,IEP(10) ,NEP(10) ,IAPTS(  10)  , 
3J0UTI 10) ,SPL(10) ,SY1 (10) ,SPLM(10) ,SY1M(10) ,DEV( 10,99) ,KEP(99). 
4 IOUT ( 99 ) ,WA(7) ,XA( 19) , LOC ( 7 , 1 9 ) , WOUT ( 99) , XOUT (99) , IS0UT(99) , 
5JPCH( 99 ) 

DATA  WA/'ST ' , ' TI  1 , 'UR' , 'CU‘ , 'WF' , ' AL' , ‘MG'/ 

DATA  XA/'ZZ'  , 1 1 A * , ' 2A 1 , ‘ 3A  1 , '4A‘ , '5A‘ , '6A‘ , '7A‘ , '8A' , 1 9A 1 , 

1 ‘IB1 , ' 2 B ' , 1 3 B 1 , 1 4B 1 , ' SB ' , 1 6B 1 , ‘ 7 B ' , '8B' , '9B1/ 

SC0NV.1 . E-3 
SMETRC=6. 89476 
NRUN=0 
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69  C 

70  C 

71  C 

72  C 

73  C 

74  C 

75  C 

76  C 

77 

78 

79  C 

80  C 

81  C 

82  C 

83  C 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93  2 

94  3 

95  C 

96  C 

97  C 

98  C 

99  C 

100  C 

101 
102 

103  4 

104 

105 

106 

107  C 

108  C 

109  C 

110  C 

111  C 

112  C 

113  C 

114 

115 

116 

117 

118 

119  5 

120 

121  6 
122 

123 

124  C 

125  C 

126  C 

127  C 

128  C 

129  C 

130  C 

131  C 

132  7 

133  8 

134 

135 

136  9 


* READ-IN  OF  MATERIAL  AND  TEST  DESCRIPTORS  AND  LOAD-STRAIN  DATA: 

* SNO*  SPEC  NO;  W,X,Y,ZMTL=  MATL  DESCRIPTORS;  BOOK , NPAGE*NOTEBOOK  REF; 

* TYPEsTENS,  COMP,  OR  SHEAR;  NDATE*TEST  DATE;  KTEMP=TEST  TEMP: 

* KEDOT * LOG  OF  STRAIN  RATE;  N = NO  OF  DATA  POINTS;  PPL=  PROP  LIMIT  LOAD; 

* AREA=SP£C  AREA  IN  SO-IN;  PULT=ULTIMATE  LOAD  IN  LBS; 

* E LONG= E LONG AT  ION  & RA=RED ' N IN  AREA  BOTH  IN  PCT; 

* RCH=ROCKWELL  C HARDNESS;  I SET=DATABANK  SET  NO. 

READ(5. 101  , END=999) , SNO, WMTL.XMTL.YMT L.ZMTL, BOOK, NPAGE. TYPE, 

1 NDATE, KTEMP, KEDOT, N, PPL, AREA, PULT.E LONG, RA.RCH, I SET 


* INITIAL  VALUES  ESTABLISHED  FOR  SELECTED  VARIABLES 

NEND=0 
DO  3 d = 1 , 10 
I EP ( d ) =0 
NEP( J ) =0 
dOUT ( J ) =0 
DO  2 L = 1 ,5 
SBAR ( d , L ) =0 . 

DBAR ( d , L ) =0  . 

DEL ( d , L ) =0 . 

CC(d, L)*0. 

CONTINUE 


* READ-IN  OF  LOAD-STRAIN  DATA  AND  CONVERSION  OF  LOAD  TO  STRESS 

* P = LOAD  IN  POUNDS;  S*  STRESS  IN  PS  I ; E*  STRAIN  IN  IN/IN 


DO  4 1=1 ,N 

READ  102.Pl  I ) ,E( I ) 

S(I)=P(I)/AREA 

PPL=  PP  L/ AREA 

SULT=PULT*SCONV/AREA 

sultm=sult*smetrc 


* S-E  DATA  MUST  BE  ORDERED  IN  INCREASING  VALUES  OF  STRAIN  MAGNITUDE 

* IF  NOT,  SUBROUTINE  ORDER  REARRANGES  S-E  DATA 

* IF  MORE  Than  N ITERATIONS  REQUIRED  FOR  ORDERING  RUN  IS  ABORTED 

M0RDER=0 

NL=N-1 

DO  5 1*1 ,NL 

I F ( E ( I ) ,LE.E( 1+1 ) ) GO  TO  5 
GO  TO  6 
CONTINUE 
GO  TO  7 

CALL  ORDER! N,S,E,P, MORDER ) 

IF(MORDER.EQ.O)  GO  TO  7 
GO  TO  31 


* ELASTIC  DATA  ANALYSIS  FOR  MODULUS  AND  INTERCEPT  USING  LSTSQ  * 

* STRAIN  INTERCEPT  SUBTRACTED  FROM  ALL  STRAINS  SO  5-£  CURVE  GOES  THRU  ZERO 

* RESULTING  NET  STRAIN  USED  IN  CALCULATIONS.  CALLED  ADdUSTED  (ADd)  STRAIN 

* EP=  PLASTIC  STRAIN  COMPONENT;  EP  USED  LATER  TO  CALCULATE  R-0  PARAMETERS 


1*0 
1 = 1 + 1 

I F ( S ( I ) . GT . PPL . OR . I . EQ. N ) GO  TO  9 
GO  TO  8 
I PL=  I 
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137  IE-IPL-1 

138  DO  10  1-1, IE 

139  X( I ) *E( I ) 

140  10  Y(I)«S(I) 

141  LI-1 

143  Ld*IE 

143  CALL  LSTSQ(X,Y,SL,OF,LI.Ld) 

144  EMOD-SL 

145  El  NT --OF/SL 

148  K.IPL 

147  NJ-N 

148  DO  1 1 1*1  ,N 

149  E( I ) =E( I )-EINT 

150  EP(I  )»E(  I)-S(  D/EMOD 

151  I F ( I . LT.K)  GO  TO  1 1 

152  IF(E( 1-1 ) .LE. .02.AND.E(I).GT. .02)  Nd-I-1 

153  IF(EP( I) .LE.O)  K-I+1 

154  I F ( EP ( I ) .LT. .00005)  Id*I+1 

155  11  CONTINUE 

156  IF(EP(K).GE. .00005)  Id-K 

157  INTCPT-EINT-1 . E+6 

158  POF-OF 

159  OFM=OF*SCQNV*SMETRC 

160  C 

161  C 

162  C 

163  C * IDENTIFICATION  OF  DATA  POINT  WHOSE  EP  VALUE  CLOSEST  TO  .002  IN/IN 

164  C * STRESS  VALUE  OF  THIS  POINT  IS  APPROXIMATE  GRAPHIC  YIELD  STRESS 

165  C 

166  C * LENO-  INDICATOR  OF  EXTENT  OF  DATA;  S(N). STRESS  VALUE  OF  LAST  DATA  POINT 

167  C * LEND-3  * S(N)  LE  GRAPHIC  YIELD 

168  C • LEND-4  • NO  MORE  THAN  3 DATA  POINTS  BETWEEN  S(N)  AND  GRAPHIC  YIELD 

169  C • LEND-5  * AT  LEAST  4 DATA  POINTS  GT  GRAPHIC  YIELD 

170  C 

171  I F ( EP ( N ) .GE. .002)  GO  TO  12 

172  I Y *N 

173  LEND-3 

174  GO  TO  16 

175  12  l-IE 

176  13  I-I-t-1 

177  I F ( EP ( I )- . 002 ) 13,14,15 

178  14  IY-I 

179  LEND*  3 

180  IF(IY.EQ.N)  GO  TO  16 

181  15  I Y* I - 1 

182  NPTS*N-IY 

163  I F ( MPT  5 . LT . 4 ) LEND*4 

184  I F ( NPTS . GE . 4 ) LEND*5 

185  C 

186  

187  C 

188  C * TEN  PAIRS  OF  TRIAL  VALUES  OF  XN  AND  SY  COMPUTED  BY  LOG-LEAST  SQUARES 

189  C 

190  C • SETTING  LOWER  BOUND  VALUES  OF  EP  FOR  LOG-LSTSQ  ANALYSES 

191  C * UPPER  60UND  IS  INPUT  VALUE  WHOSE  EP  IS  CLOSEST  TO  .002 

192  C • IF  FEWER  THAN  8 PTS  FOR  LOG-LSTSQ,  UPPER  BOUND  EXTENDED  TO  GET  B PTS 

193  C 

194  16  IJ-Id-1 

195  DO  17  d*1 ,10 

196  IJ*IJ+1 

197  I EP( d ) * Id 

198  NEP(d)*EP(Id)-1 .E6 

199  17  CONTINUE 

200  C 

202  C 

203  C • LOG-LEAST  SQUARES  ANALYSES  DONE  BY  SUBROUTINE  LSTSQ 

204  C 
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2C5 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 


DO  18  I =K , N 
X(I)=ALOG(EP(I)) 

18  Y(I)=ALOG(S( I) ) 

CONST  - ALOG( .002) 

Ld*  I Y 

DO  19  d = 1 ,10 
LI  * I EP  ( J ) 

NPTS=Ld-LI 

IF(NPTS.LT.7)  Ld- LI+7 
IF(Ld.GT.N)  Ld-N 
IAPTSl d)=Ld-LI+1 
CALL  LSTSQ(X,Y,SL,OF,LI,Ld) 

AS=OF+SL*CONST 
SY ( d ) -EXP ( AS ) 

XN ( d ) = 1 ,/SL 

19  CONTINUE 
C 



c 

C * COMPUTATION  OF  ANALYTIC  STRESS  VALUES  FOR  EACH  INPUT  STRAIN  UP  TO  EPS-. 02 
C * DONE  USING  A TRIAL  AND  ERROR  PROCEDURE  BY  SUBROUTINE  RdA75 
C 

C * COMPUTATION  ABORTED  FOR  ANY  ONE  OF  THREE  CONDITIONS 
C • NO  SOLUTION  IN  100  ITERATIONS:  ITER-100  & dOUT-1 
C • INCREMENTED  VALUE  OF  Z TOO  LARGE:  ITER-200  « dOUT=2 
C • INCREMENTED  VALUE  OF  Z LESS  THAN  .01!  ITER-300  & d0uT»3 
C * NEND  COUNTS  THE  NUMBER  OF  d 1 S THAT  ABORT  • 

C 

C • ABORTS  OCCUR  ONLY  IF  MAGNITUDES  OF  INPUT  DATA  ARE  GROSSLY  IN  ERROR 
C 

LI  -K 

DO  21  d = 1 , 1 0 
B-EMOD/SY ( d ) 

A A- . 002  * B 
XP-XN(d) 

DO  20  I = L I , Nd 
Z-1 

BB-B*  E ( I ) 

CALL  RdA75IZ,XP,BB,AA, ITER) 

IFUTER.LT. 100. AND. ITER. GE.1)  GO  TO  20 
IF(ITER.EQ.O)  dOUT(d)-9 
I F ( ITER. EO- 100)  dOUT ( d ) — 1 
I F ( ITER. EQ. 200)  dOUT(d)=2 
IF(ITER.EQ. 300 ) dOUT ( d ) =3 
I OUT ( d ) - I 
NEND-NEND+1 
GO  TO  21 

20  CS ( d , I ) =Z«SY ( d ) 

21  CONTINUE 
C 

c 

C • STATISTICAL  EVALUATION  OF  FIT  OF  TEST  DATA 
C 

C * CC*  CORRELATION  COEFFICIENT  j SEE-  STANOARD  ERROR  OF  ESTIMATE 
C • PREFIX  E*  ELASTIC  DATA 
C 

C • L IS  INDICATOR  OF  REGION  OF  S-E  CURVE  WHICH  IS  STATISTICALLY  EVALUATED 
C » MIDPT-  NUMBER  OF  DATA  POINT  CLOSEST  TO  MEDIAN  STRESS  BETWEEN  SPL  A SY 
C * L*1  * SPL  TO  MIDPT 
C * L-2  * MIDPT  TO  SY 

C • L-3  » SPL  TO  SY  (OR  TO  S(N)  IF  S(N)  LT  SY) 

C * L-4  • SPL  TO  S ( N ) (USED  ONLY  IF  LEND  GE  4) 

C • L-5  • SY  TO  S ( N ) (USED  ONLY  IF  LEND  EO  5) 

C 

C • ELASTIC  MODULUS  EVALUATION 
C 

ECC-0 . 
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273  ESBAR-O. 

274  ED9AR-0. 

275  EDEL.O. 

276  Ld«K-1 

277  DO  22  1.1 .Id 

278  22  ESBAR«ESBAR+S( I )/LJ 

279  DO  23  I«1  , L J 

280  ECS»E( I )»EMOD 

281  EDBAR«EDBAR+(S(I)“ESBAR)»*2 

282  23  EDEL*EDEL+< S( I )-ECS)**2 

283  ESEE*SQRT(EOEL/LJ) 

284  COR-1 .-EDEL/EDBAR 

285  IF(COR.GT.O)  ECC*SQRT(COR) 

286  C 

287  C * R-0  EQUATION  EVALUATION:  CALCULATIONS  FOR  EACH  SET  OF  TRIAL  PARAMETERS 

288  C 

289  C * UPPER  LIMIT  OF  INDEX  I IS  N,  IF  E(N)  LESS  THAN  .02 

230  C * UPPER  LIMIT  OF  I IS  NJ.  IF  E(N)  GREATER  THAN  .02 

291  C • E(NJ)  IS  VALUE  OF  STRAIN  JUST  BELOW  OR  EQUAL  TO  .02 

292  C 

293  I F ( NEND . EQ . 1 0 ) GO  TO  29 

294  DO  27  J«1 ,10 

295  I F ( JOUT ( J ) . GT . 0 ) GO  TO  27 

296  MIDPT«K+( IY-KJ/2 

297  DO  26  L = 1 .LEND 

298  I F ( L . NE . 2 . AND. L • NE . 5 ) LI»K 

299  I F ( L . FQ . 2 ) LUMIDPT+1 

300  I F ( L . EQ . 5 ) LI=IY+1 

301  IF(L.EO.I)  LJ=MIDPT 

302  IF( L.EQ.2.0R. L.EQ.3)  LJ«IY 

303  I F ( L . GE . 4 ) LJ«Nd 

304  ISPTS(d, L)=Ld-LI+1 

305  DO  24  I.LI.Ld 

306  24  SBAR(d.L)»SBAR(d,L)+S(I)/ISPTS(d,L) 

307  CO  25  I.LI.Ld 

3C8  DBAR (J,L)=OBAR(J,L)+(S( I ) -S8AR ( d , L ) ) **2 

309  25  DEL(d,L)»DEL(d,L)+(S(I)-CS(d,I))**2 

310  SEE ( d , L ) »SQRT (DEL(J,L)/ISPTS(d,L)) 

311  C0R.1 .-DEH J,L)/DBAR(d,L) 

312  IF(COR.GT.O)  CC( d , L ) »SQRT (COR ) 

313  26  CONTINUE 

314  27  CONTINUE 

315  C 

316  C 

317  C 

318  C * BEST-FIT  CRITERION:  TRIAL  VALUES  OF  XN  A SY  WHICH  GIVE  MAX  CC  FOR  SPL  TO  SY 

319  C 

320  , I F (NEND . EQ. 1 0)  GO  TO  29 

321  d 1 * 1 

322  C0R=0. 

323  DO  28  J«1 ,10 

324  I F ( JOUT ( d ) • GT . 0 ) GO  TO  28 

325  CR I T*CC( d , 3) 

326  IF(CRIT.LE.COR)  GO  TO  28 

327  COR*CRIT 

328  d1*J 

329  28  CONTINUE 

330  C 

332  C 

333  C * UNITS  CONVERSION  FOR  PRINT-OUT  OF  NPRINT*0  OR  1 AND  PUNCH-OUT 

334  C • S A SEE  CONVERTED  TO  KSI  A MPA;  MODULUS  TO  MEGA-PSI  A GPA 

335  C • TEST  TYPE  CODED  FOR  PUNCH-OUT:  TYPE  TO  NTYPE 

336  C • PROP-LIM  (SPL)  AND  .1%~YLD  ( SY1 ) CALCULATED 

337  C 

338  29  EMOD*EMOD«1 . E-6 

339  EMODM-£mOO*SMETRC 

340  ESEE*ESEE*SCONV 
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367 
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369 
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371 
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373 
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376 
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380 

381 
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383 

384 

385 

386 

387 

388 

389 

390 

391 

392 

393 

394 

395 

396 

397 

398 

399 

400 

401 

402 

403 

404 

405 

406 

407 

408 


eseem«esee*smetrc 

NTYPE=0 

IF (TYPE.EQ. 'TENS  ')  NTYPE-1 
IFITYPE.EQ. 'COMP  ')  NTYPE«2 
IF(TYPE.EQ. 'SHEAR' ) NTYPE*3 
IF(NPRINT.GT. 1 ) CO  TO  37 
I F ( NEND . EQ • 10)  GO  TO  31 
SY ( J 1 >*SY( dl ) *SCONV 
SYM( J1 ) = SY  v dl )»SMETRC 
XP=1 ,/XN( dl ) 

SPL(J1 ) *SY ( d 1 ) • . 005*  *XP 
SPLIK(d1  )=SPL(d1 ) • SMETRC 
SYM  J1  ) *SY ( d 1 )*.5**XP 
SY 1M(J1)*SY1(J1 )* SMETRC 
DO  30  L=  1 .LEND 
SEE ( J 1 , L ) =SEE ( dl , L)*SCONV 

30  SEEM ( d 1 , L ) = SEE ( d 1 ,L)« SMETRC 
C 

C 

c 

C » PRINT-OUT  AND  PUNCH  OUT  ROUTINES 

C * DATA  BANK  PREPARATION  AND  PUNCH-OUT  START  AT  STATEMENT  NUMBER  65 
C 

C 

C 

C * NPR I NT  = 0 : PRINT-OUT  OF  BEST-FIT  PARAMETERS  * FOUR  SETS  PER  PAGE 
C * NPR I NT  = 1 : PI  REPORT  * BEST-FIT  PARAMETERS,  STATISTICS,  AND  TABLE 
C 

31  I F ( MOD (NPX,4).EQ.O.AND.NPRINT.EQ.O)  PRINT  200 

IF(NPRINT.EQ. 1 ) PRINT  208, AREA 
KPRT=KPRT+1 

PRINT  201 ,SN0,WMTL,XMTL,YMTL,2MTL, TYPE, NOATE, BOOK, NPAGE.ISET, 
1KTEMP,KED0T,SULT,SULTM,EL0NG,RA,RCH 
IF(MORDER.EQ.O)  GO  TO  32 
PRINT  300 

IF(MOD(NRUN, 4) .EQ.0.0R.NPRINT.EQ.1 ) PRINT  218 
GO  TO  1 

32  PRINT  202 

PRINT  203, INTCPT, EMOD .ESEE.ECC.IE 

IF(NEND. LT. 10)  PRINT  204 , NEP( J 1 ) , SY ( d 1 ) ,XN(J1 ) , IAPTS(J1 ) ,SPL(J1 ) , 
1 S Y1 ( J 1 ) 

PRINT  205  EMDDM, ESEEM 
IFfNEND. LT. 10)  GO  TO  33 
PRINT  301 

IF(M0D(NRUN,4) . EQ • 0 . OR . NPR INT . EQ . 1 ) PRINT  218 
IF(NPRINT.EO. 1 ) GO  TO  53 
GO  TO  65 


33  KIN=0 

PRINT  206 , S YM( d 1 ) , SPLM( d 1 ) , SY 1 M( d 1 ) 

IF(NEND.EQ.O)  GO  TO  35 
DO  34  d=1 , 10 

IF(dOUT(d) .EQ.O)  GO  TO  34 

IF(KIN.EO.O)  PRINT  302 , d , dOUT ( d ) , IOUT ( d ) 

IF(KIN.EO.I)  PRINT  303,d,dOUT(d) , IOUT(d) 

IF(KIN.E0.2)  PRINT  304 , d , dOUT ( d ), IOUT ( d ) 

I F ( K I N . EO . 3 ) PRINT  305 , d , dOUT ( d ), IOUT ( d ) 

IF(KIN.E0.4)  PRINT  306 , d , dOUT ( d ) , IOUT ( d ) 

IF(KIN.E0.5)  PRINT  307,d,dOUT(d) , IOUT(d) 

IF(KIN.EO.S)  PRINT  308 , d , dOUT ( d ), IOUT ( d ) 
lF(KIN.EO.V)  PRINT  309 , d , dOUT ( d ), IOUT ( d ) 

IF(KIN.EO.B)  PRINT  310,d,dOUT(d) , IOUT(d) 

K I N*K I N+1 

34  CONTINUE 

35  IF(KIN.EO.O)  PRINT  207 
PRINT  2 1 0 

PRINT  21 1 . ( SEE ( d 1 , L) ,CC(d1 . L) , ISPTStdl , L) ,L*1 ,3) 

I F ( LEND . EQ ■ 4 ) PRINT  2 1 2 . SEE ( d 1 , 4 ) , CC ( dl , 4 ) , I SPTS ( dl ,4 ) 

I F ( LEND . EQ . 5 ) PRINT  213, (SEE(d1 , L ) , CC ( d 1 ,L),ISPTS(d1,L), L*4 ,5) 
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PRINT  214. (SEEM(d1 , L) ,L«1 ,3) 

I F ( LENO . EQ • 4 ) PRINT  21 5 . SEEM( J 1 , 4 ) 

I F ( LEND . EQ • 5 ) PRINT  21 6 . ( SEEM( J 1 . L) , L-4 , 5 ) 

PRINT  217 

IF(M0D(NRUN,4) . EQ . 0 . OR . NPRINT . EQ . 1 ) GO  TO  36 
GO  TO  65 

36  PRINT  218 
IF(NPRINT.E0.1  ) GO  TO  53 
GO  TO  65 

C 

C 

c 

C * ALTERNATE  FULL-DATA  PRINT-OUT  AND  OPTIONAL  TABLES:  NPRINT. 2 OR  3 

C 

C 

C • DIAGNOSTIC  MESSAGE  IF  ANALYTIC  STRESSES  CANNOT  BE  CALCULATED 
C 

37  IF(NEND. LT. 10)  GO  TO  38 
PRINT  400 

PRINT  401 

PRINT  201 .SNO.WMTL.XMTL.YMTL.ZMTL, TYPE, NOATE, BOOK, NPAGE.ISET, 

1 KTEMP , KEDOT , SULT , SULTM, ELONG, RA , RCH 
PRINT  402.EMOD.POF, I NTCPT , ESEE , ECC , IE 
PRINT  403 , EMODM , OFM , ESEEM 
PRINT  301 

IF(NPRINT.EQ.3)  PRINT  315 
PRINT  217 
GO  TO  65 
C 

c 

c 

C * UNITS  CONVERSION  FOR  PRINT-OUT  OF  NPRINT.2  OR  3 
C * PER-CENT  DEVIATIONS  CALCULATED  ONLY  FOR  NPRINT-3  • 

C 

38  DO  41  d.1 , 10 

I F ( dOUT ( d ) . GT  . 0 ) GO  TO  41 
SY ( J ) «SY( J ) *SCONV 
SYM( J).SY(J) *SMETRC 
XP»1 . /XN( J ) 

SPL(J)=SY(J)*.005**XP 
SPLM( d).SPL( J)*SMETRC 
SY1 (d)«SY(d)*.5*»XP 
SY 1 M ( J ) «SY1 ( J)*SMETRC 
DO  39  L-1 .LEND 
SEE(d,L)-SEE(d,L)*SCONV 
SEEM(d, L)«SEE(d, L)*SMETRC 

39  CONTINUE 
1F(NPRINT.E0.1 ) GO  TO  41 
DO  40  I.K.Nd 

DEV(d,I)«(CS(d,I)-S(I))*100./S(I) 

CS(d,I)=CS(d,I ) *SCONV 

40  CONTINUE 

41  CONTINUE 
IF(NPRINT,EQ.2)  GO  TO  43 
DO  42  I.K.Nd 
KEP(I)*EP(I)»1 . E6 

S( I )«S( I J.SCONV 

42  CONTINUE 
C 



c 

C • PRINT-OUT  OF  all  TRIAL  PARAMETERS  • ONE  PAGE  PER  DATA  SET:  NPRINT-2 
C 

43  PRINT  400 
PRINT  401 

PRINT  201 .SNO.WMTL.XMTL.YMTL.ZMTL, TYPE, NOATE, BOOK, NPAGE.ISET, 

1 KTEMP. KEDOT. SULT, SULTM. ELONG. RA, RCH 
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478 

479 
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483 
434 

485 
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541 

542 
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544 


PRINT  402,EMOD,POF, INTCPT.ESEE.ECC, IE 
PRINT  403 , EMODM , OFM , ESEEM 
PRINT  404 
DO  45  d«1 ,10 

I F ( JOUT ( J ) . GT . 0 ) GO  TO  44  / 

PRINT  40S,d,NEP(d) ,XN(d) , S'fCa ),SYM(d),SPL(d),SPLM(d), 

1SY1 (d) ,SYlM(d) ,IAPTS(d) 

IF(d. EO.dl ) PRINT  406 
CO  TO  45 

44  PRINT  3l6,d,NEP(d),dOUT(J) 

45  CONTINUE 

c 

C * PRINT-OUT  OF  STATISTICAL  DATA  FOR  ALL  TRIAL  VALUES  OF  PARAMETERS  * 

C 

PRINT  409 
DO  47  d = 1 , 1 0 

I F ( dOUT ( d ) ■ GT . 0 ) GO  TO  46 

PRINT  4lO,d,NEP(d) , ( SEE ( J , L) , SEEM ( d , L) , CC ( d, L) , I SPTS ( d , L ) , L«1 ,3) 

I F ( LEND . EQ • 3 ) PRINT  411 

I F ( LEND . EQ . 4 ) PRINT  4 1 2 . SEE ( d , 4 ) , SEEM( d , 4 ) , CC ( d , 4 ) , I SPTS ( d , 4 ) 

I F ( LEND . EQ . 5 ) PRINT  4 1 3 , ( SEE ( d , L ) , SEEM( d , L ) , CC ( d , L ) , I SPTS ( d . L ) , 

1 L=4 , 5 ) 

GO  TO  47 

46  PRINT  3l7,d,NEP(d) , dOUT (d) 

47  CONTINUE 
C 

C 

C * PRINT-OUT  OF  DIAGNOSTIC  MESSAGES  IF  ANALYTIC  STRESS  CALCULATIONS  ABORTED 
C * MESSAGES  APPEAR  AT  BOTTOM  OF  ONE-PAGE  PRINTOUT  OF  PARAMETERS  * 

C 

IF(NEND.GT.O)  PRINT  318 
PRINT  217 

IF(NPRINT.EQ.2)  GO  TO  65 
C 

C 

C * PRINT-OUT  OF  OPTIONAL  TABLES:  NPRINT*3 
C * TABLES  IN  TWO  GROUPS:  d = 1 TO  d=5  AND  d«6  TO  d«10 
C * INDEX  LIMITS  ARE:  dK=1  TO  dL*5  AND  UK*6  TO  dL»10 
C 

dK  = 1 


48 


49 


50 


51 


dL=dK+4 
PRINT  414, SNO 
IF(dK.EQ.I)  PRINT  415 
I F ( dK . EQ . 6 ) PRINT  4l6 
PRINT  417,(NEP(J) , J = JK , J L ) 
K0U1IT  = 0 


IF(NEND.GT.O)  GO  TO  50 
DO  49  t *K , Nd 
KOUNT  = KGUNT  + 1 


PRINT  4l8.I,E(I),S(I),KEP(I),(CS(d,I),DEV(dtI),d«dK,dL) 

I F ( MOD( KOUNT , 5 ) . EQ . 0 ) PRINT  207 

CONTINUE 

I F ( UK . EQ ■ 6 ) GO  TO  65 
dK*6 

GO  TO  48 
DO  52  I — K , N d 
KOUNT  = KOUNT  + 1 

PRINT  419, I ,E( I ) ,S( I ) ,KEP( I ) 

00  51  d?dK  , d L 

I F ( dOUT ( d ) . GT . 0 ) GO  TO  51 

IF(d.EO.dK)  PRINT  420,CS(d, I ) , DE V ( J , I) 

I F ( d . EQ . ( dK+ 1 ) ) PRINT  421  ,CS(d,I) ,OEV(d,I) 

IFfd.EQ. (dK  + 2) ) PRINT  422 , CS( d , I ) , DEV ( d , I ) 

IFId.EQ. (dK+3) ) PRINT  423 , CS( d , I ) , DEV ( d , I ) 

I F ( d . EQ  . ( dK+4 ) ) PRINT  424 , CS( d > X ) , DEV ( d , I ) 

CONTINUE 
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I F ( MOD ( KOUNT , 5 ) . EO • 0 ) PRINT  207 

52  CONTINUE 
PRINT  425 

I F ( JK . EQ . 6 ) GO  TO  65 
JK*6 

GO  TO  48 
C 

C 

C * PI  REPORT  UNITS  CONVERSION  AND  TABULAR  PRINT-OUT:  NPRINT-1 
C 

53  I F ( KPRT . GT . 1 ) GO  TO  57 
PEMOO=EMOO*1 • E6 

I F ( N . LE . 60  ) I J*N 
I F ( N . GT . 60 ) I J = 60 
DO  56  I > 1 , 1 J 
I F ( I . GE . K ) GO  TO  54 
CS(J1 , I ) * PEMOD*  E l I ) 

KEP ( I )«0 

54  IF(I.GE.K)  KEP ( I ) * EP ( I ) * 1 . E6 

IF( I .GE.K. AND.NEND.EQ. 10)  GO  TO  55 

I F ( I .C.T.Nd)  GO  TO  55 

DEV(U1 ,I)s(CS(J1 ,1 )— S (I))*100./S(I) 

CS(U1 . I )«CS( J1 , I )»SCONV 

55  S(I)iS(I)*SCONV 

56  CONTINUE 

57  PRINT  217 

IF(NENO.EQ. 10)  PRINT  319, SNO 

IF(NEND. LT. 10)  PRINT  426 , SNO. SY( J1 ) , XN( dl ) , EMOD 

I F ( N . GT . 30 ) GO  TO  59 

PRINT  427 

DO  58  I«1 ,N 

IF(I.LT.K)  PRINT  428, 

1I,P(I),E(I).KEP(I),S(I).CS(d1,I),DEV(J1,I) 

IFtl. GE.K. AND. I. LE.NJ.AND.NEND.lt. 10)  PRINT  428, 

1I,P( I)  ,E(I) ,KEP(I) ,S(I) ,CS(J1 ,1) ,DEV(d1 ,1) 

IF( I .GE.K. AND.NEND.EQ. 10)  PRINT  429, 

1I,P(I),E(I),KEP(I),S(I) 

IF( I .GE.K.AND. I .GT .NJ. AND.NEND. LT.  10)  PRINT  429, 
1I,P(I),E(I),KEP(I),S(I) 

58  CONTINUE 
IF(NEND.GT.O)  PRINT  318 
IFfKPRT. LT.NCOPY)  GO  TO  31 
GO  TO  65 

59  PRINT  430 

I F(N  . LE  . 60  ) LJ  = N 
I F ( N . GT . 60 ) Ld*60 
DO  64  I»1 ,30 
I P*30+ I 

IF(IP.LE.Ld)  GO  TO  61 
IF(I.GT.NJ)  GO  TO  60 
IFfNEND.EQ.10. AND . I . GE . K ) GO  TO  60 

PRINT  428,I,P(I),E(I),KEP(I),S(I),CS(J1,I) ,DEV( dl , I ) 

GO  TO  64 

60  PRINT  429,1 ,P( I) ,E(I) ,KEP(I),S(I) 

GO  TO  64 

61  I F ( NEND . EQ . 1 0 . AND . I . GE . K ) GO  TO  62 
IFfNEND.LT.10. AND. I . LE. NJ . AND . I P.GT .Nd ) GO  TO  63 
PRINT  431 .I.P(I) ,E(I),KEP(I),S(I).CS(d1 .1), 

1DEV(J.,I),IP,P(IP),E(IP),KEP(IP),S(IP),CS(J1,IP) , DEVf  d 1 , IP) 

GO  TO  64 

62  PRINT  432,I.P(I),E(I),KEP(I),S(X), 
1IP,P(IP),E(IP),KEP(IP),S(IP) 

GO  TO  64 

63  PRINT  433.1 , P( I ) , E( I ) ,KEP( I ) , S( I ) ,CS( J1 , I ) ,OEV( dl , I ) , IP,P( IP) 
1E(IP),KEP(IP),S(IP) 

64  CONTINUE 

IF(N . GT . 60 ) PRINT  434, N 
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613  IF(NEND. GT.O. AND. N. LE.60. AND. Nd.GE.N)  PRINT  318 

614  IF(N. GT. 60. AND. N.GT.Nd. AND. Nd. LE.60)  PRINT  435 

615  IF(N. LE.60. AND. NJ. LT.N)  PRINT  436 

616  IF(KPRT. LT.NCOPY)  GO  TO  31 

617  C 

618  C 

619  C 

620  C * PREPARATION  OF  DATA  FOR  DATA  BANK 

621  C * ASSIGNMENT  OF  NEXT  VALUES  AND  UPDATE  OF  LOCATER  VALUES 

622  C 

623  65  I F ( N PUNCH . EQ . 0 ) GO  TO  1 

624  NEXT  *0 

625  JW=0 

626  JL»0 

627  DO  66  dX=1 .7 

628  IF(WMTL.NE.WA(JX) ) GO  TO  66 

629  JW  = dX 

630  GO  TO  67 

631  66  CONTINUE 

632  67  IF(JW.EQ.O)  GO  TO  70 

633  DO  68  LL=1 , 19 

634  1 F ( XMT  L.NE.XA(LL))  GO  TO  68 

635  dL  = LL 

636  GO  TO  69 

637  68  CONTINUE 

638  IF(JL.EQ.O)  GO  TO  70 

639  69  NEXT=LOC(JW, JL) 

640  LOC( JW , JL )= I SET 

641  70  I F ( JW . EQ . 0 . OR . JL . EQ . 0 ) JPCH(NRUN)=1 

642  C 

643  C 

644  C 

645  C * PUNCH-OUT  OF  DESCRIPTORS  AND  R-0  PARAMETERS 

646  C 

647  IF(NEND. LT. 10)  PUNCH  500 , SNO , KMTL , XMT L , BOOK , NPAGE , NTYPE , NDATE , 

648  1K1EMP.KED0T.SULT,EL0NG,RA,RCH,EM00,XN(J1 ) ,SY(J1 ) ,SEE(J1 ,3), 

649  1NEXT.1SET 

650  IF(NEND.EQ.IO)  PUNCH  50 1 , SNO , WMT L , XMT L . BOOK , NPAGE , NTYPE , NDATE , 

651  1KTEMP.KED0T . SULT , E LONG , RA , RCH , EMOD , NEXT , I SET 

652  IF(JPCH(NRUN) . EQ . 0 ) GO  TO  1 

653  KPCH=KPCH+1 

654  WOUT ( NRUN ) = WMT  L 

655  XOUT ( NRUN)=XMTL 

656  I SOUT ( NRUN ) = I SET 

657  GO  TO  1 

658  C 

660  C 

661  C * PUNCH  OUT  OF  UPDATED  LOCATER  FILE 

C62  C 

663  999  IF(NPUNCH.EQ.O)  STOP 

664  IF(NPRINT.EQ.O.AND.MOO(NPX.4) .NE.O)  PRINT  218 

665  DO  71  JX*1 , 7 

656  PUNCH  502,WA(dX)t( LOC (JX,JK),JK«1,19) 

667  71  CONTINUE 

668  lF(KPCH.EQ.O)  STOP 

669  C 

670  C 

671  C 

672  C • PRINT-OUT  OF  DIAGNOSTIC  MESSAGE  FOR  DATA  CODE  NOT  IN  DATA  BANK 

673  C 

674  PRINT  320 

675  DO  72  NR=1 ,NPX 

676  IF(dPCH(NR) . EQ . 0 ) GO  TO  72 

677  PRINT  321,1 SOUT ( NR ) , WOUT ( NR ) , XOUT ( NR ) 

678  72  CONTINUE 

679  STOP 

680  C 


20 


681 

682 

683 

684 

685 

686 
687 
683 
609 

690 

691 

692 

693 

694 

695 

696 

697 

698 

699 

700 

701 
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703 
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710 
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C 

C * READ-IN  AND  PRINT-OUT  FORMAT  STATEMENTS  * 

C 

C 

C * DATA  READ-IN  FORMAT  STATEMENTS 
C 

100  FORMAT (312) 

101  FORMAT( A6.2A2.3A4, 12, A5, 16, 15 , 2 12 , F7 . 0 , F6 . 4 , F7 . 0 , 3F4 . 1 ,14) 

102  FORMAT(F7.0,F7.5) 

C 

C 

C * FORMAT  STATEMENTS  FOR  BEST-FIT  PRINTOUT:  NPRINT-0  OR  1 
C 

C * TEST  DATA  AND  R-0  PARAMETERS:  FORMAT  STATEMENTS  200  TO  207 
C 

200  FORMAT! 1 HI ,T22,  1 RAMBERG-OSGOOD  BEST  FIT  PARAMETERS  * 0.2X  OFFSET  Y 
HELD  AND  EXPONENTS' ,/,T4,124( '-') ) 

201  F0RMAT(T12,  'SPECIMEN  NO:  ' , A6 , 2X , ' MATER  I AL : ',A2t'-'tA2, 

1 '-' ,2A4,2X, 'TEST  TYPE:  ' , A5 , 2X , ' DATE : ' , 1 6 , 2X , ' NTBK-PG : '.A4, 

2'-' , I2.2X, 'SET  NO:  ' , 14 , / , T1 2 , ' TEMP- ' , 15 , ' C ' , 2X . ' EDOT« 1 . 0 ' , 12 , 

3'  PER-SEC' ,2X,  'S-ULT-  '.F7.2,'  KSI  ( ' . F7 . 1 , • MPA )', 2X ,' E LONG- 
4F4.1 , '%' ,2X. 'RED— A=  ' ,F4.1 ,2X, 'R-CH*  ',F4.1,/) 

202  FORMAT(T20, 'UNITS  : INCPT ', 2X ,' MODULUS 4X ,' SEE ', 4X ,' C-CFT ', 2X , 
I'PTS  : ' , 1 X , ‘ EP-I ' , 2X , ' . 2%-YLD ' ,3X, 'EXPNT' ,2X, 'PTS  :'t3X,'SPL't 
23X , ' . 1%-YLD' ) 

203  FORMAT (T20. ‘ENGL  : 1 , 1 X , 15 , 2X , F6 . 2 , 3X , F5 . 2 ,2X , F7 .4 , 2X , 13 , 1 X , ' : ' ) 

204  FORMAT  (1H-F.T20.46X,  14 , 2X  , 2 ( F7 . 2 , IX)  , 1 X,  13,  IX,  ' : 1 , 2(  IX,  F7.2)  ) 

205  FORMAT(T20.  'METRC  : ' , 8X , F6 . 1 , 3X , F5 . 1 , 1 5X , ' : • ) 

206  FORMAT! 1 H+ , T20 , 52X , F7 . 1 , 1 4X , ' : ' , 2 ( 1 X , F7 . 1 ) ) 

207  FORMAT ( 1 H ) 

C 

c 

C ♦ PI  REPORT  HEADING:  NPRINT*1 
C 

208  F0RMAT(1H1 ,T20, '*»PRINCIPAL  INVESTIGATOR  REPORT  ON  STRESS-', 

1 'STRAIN  TEST  AND  RAMBERG-OSGOOD  ANALYSIS  *♦',/, T 1 0 .' BEST-FIT 

2 ' MODULUS , 0.2%  OFFSET-YIELD,  AND  R-0  EXPONENT  CALCULATED  BY  ', 
3’LEAST  SQUARES.  (SPEC  AREA-' . F6.4, ' SQ-IN)',/) 

C 

C 

C * STATISTICS  FORMAT  STATEMENTS  210  TO  210 
C 

210  FORMAT ( T 1 6 4X,'SPL  TO  MPT', 8X, 'MPT  TO  SY',9X,'SPL  TO  SY'.9X,'SPL 

1 TO  SN'.IOX.'SY  TO  SN ' , / , T 1 6 , 1 OX , 5 ( 2X , ' SEE ' , 3X , ' C-CFT ' , 1 X , ' PTS • . 1 X ) 
2) 

211  F0RMAT(T16, ' ENGL  ' , 5X , 3 ( F6 . 2 , F7 . 4 , 1 X , 13 , 1 X ) ) 

212  FORMAT! 1 H+ , T 1 6 . 64X , F6 . 2 , F7.4, IX, 13) 

213  FORMAT!  1H-f,T16,64X,2(F6.2,F7.4,1X,I3,1X)  ) 

214  FORMAT! T 16, ' METRC ' , 5X , 3 ( F6 . 1 , 1 2X ) ) 

215  FORMAT! 1H+,T16,64X,F6. 1 ) 

216  FORMAT!  1H-f,T16,64X,2(F6. 1 ,12X)) 

217  F0RMAT(T4, 1 24 ( '-' ) ) 

218  FORMAT ( T1 0 , ' **  UNITS:  MODULUS  IN  M-PSI  OR  GPA  • ,2%-YLD,  .1X-YLD, 

1 SPL  & SEE  IN  KSI  OR  MPA  * INTCPT  AND  EP-I  IN  MICRO-IN/lN  *• ' ) 

C 

c 

C * DIAGNOSTIC  MESSAGE  PRINT-OUT  FORMAT  STATEMENTS  300  TO  321 
C 

300  FORMAT!////, T28, DATA  FOR  THIS  SPECIMEN  COULD  NOT  BE  ORDERED. 

1 ANALYSES  WERE  ABORTED ,',//////, T4 , 1 24 ('-') ) 

C 

C * STATEMENTS  301  TO  321  ARE  PRINTED  WHEN  R-0  EQUATION  CANNOT  BE  SOLVED 
C 


21 


750  301  FORMAT*//, T24, '***  THEORETICAL  STRESSES  COULD  NOT  BE  CALCULATED  FO 

751  1 R ANY  OF  THE  R-0  PARAMETERS.  **•',/, T24 , ' STAT I ST ICAL  EVALUATION  WA 

752  2S  ABORTED  AND  NO  BEST-FIT  PARAMETERS  COULD  BE  DETERMINED.',////, 

753  3T4 , 1 24 ( ' — 1 ) ) 


754 

302 

FORMAT ( T5 , ' * S-ABORT : 

d='  . 

12, 

- ' 

,11. '-' 

.12) 

755 

303 

FORMAT ( 1 H+ , T5 , 20X , '* 

J*  ' , 

12, 

,11,  '-' 

,12) 

756 

304 

FORMAT* 1H+, T5, 32X , 1 • 

d=  ' . 

12, 

“ ' 

.11.'-' 

,12) 

757 

305 

FORMAT* 1H+, T5.44X,  ' * 

d=  ' , 

12. 

“ 1 

.11.'-' 

,12) 

758 

306 

FORMAT ( 1H+.T5.56X, ' • 

J*’  , 

12. 

“ ' 

.11,'-' 

,12) 

759 

307 

FORMAT ( 1 H+ , T5 , 68X , '* 

d = ' . 

12, 

.11.'-' 

,12) 

760 

308 

FORMAT* 1H+. T5.80X, ' * 

J*  ' . 

12, 

“ ' 

, 11  . 

,12) 

761 

309 

FORMAT ( 1 H+ , T5 , 92X , '• 

J*  ' . 

12, 

,11.'-' 

,12) 

762 

310 

FORMAT ( 1H+, T5, 104X, ' 

* d=  ' 

.12 

• _ 

' .11  , '- 

M2) 

763 

C 

764 

C * 

STATEMENTS  315  TO  318 

REFER 

TO 

NPRINT  *2 

OR  3 

765 

766 

767 

768 

769 

770 

771 

772 

773 

774 

775 

776 

777 

778 

779 

780 

781 

782 

783 

784 

785 

786 

787 

788 

789 

790 

791 

792 

793 

794 
'795 

796 

797 

798 

799 

800 
801 
802 

803 

804 

805 

806 

807 

808 

809 

810 
811 
812 

813 

814 

815 

816 

817 

818 
819 


C 

315  F0RMAT(T5, 'SINCE  ANALYSES  WERE  ABORTED,  REQUESTED  TABULAR  DATA  WER 
IE  NOT  AVAILABLE  AND  COULD  NOT  BE  PRINTED',/) 

316  FORMAT (T25, 13, 4X, 15, 6X,  ' * ABORT * ANALYTIC  STRESSES  COULD  NOT  BE  CAL 
1 CULATED  * ABORT  CODE*', 12) 

317  F0RMAT(T6,I3,IS,3X, ' *ABORT*  ANALYTIC  STRESSES  COULD  NOT  BE  CALCULA 
1 TED  * ABORT  CODE* ' , 12) 

318  F0RMAT(T4,124( ,/,T4, 'ABORT  MESSAGE:  J VA LUE-CODE-PT  NO1, 

1 2X , ' * ABORT  CODE:  1*  NO  CONVERGENCE*  2*  S TOO  LARGE*', 

2'  3=  S TOO  SMALL*  9*  R-0  EXPNT  LESS  THAN  1') 

C 

C * PI  REPORT  TABLE  DIAGNOSTIC  WHEN  R-0  ANALYSIS  ABORT 
C 

319  FORMAT* 1H0.T10. 'SPEC  NO:  1 , A6 , 2X , '‘DIGITIZED  STRESS  AND  STRAIN  ', 

1 'VALUES  ANO  THEORETICAL  STRESS  DEVIATIONS  ', 

2 ' FOR  ELASTIC  DA T A ' , / , T 1 0 , ' THEORET IC A L DATA  BEYOND  THE  PROP', 
3'0RTI0NAL  LIMIT  WERE  NOT  AVAILABLE',/) 

C 

C*  DIAGNOSTIC  WHEN  TEST  MATERIAL  CODE  NOT  ASSIGNED  IN  DATA  BANK 
C 

320  F0RMAT(1H1 ,T30,36( '*', IX) ,//,T30, 'MATERIAL  ANO/OR  ALLOY  CODES', IX, 
1 ' HAVE  NOT  BEEN  ASSIGNED  IN  THE  DATA  BANK ',/, T30 ,' FOR  THE  TEST', IX, 
2 ' SETS  LISTED  BELOW.  AFTER  CODES  ARE  ASSIGNED,  NEXT  VA LUES  1 , / , T30 , 

3 ' SHOULD  BE  ASSIGNED  AND  PUNCHED  AND  THE  LOCATER  FILE  SHOULD  BE'. 

4'  UPDATED' ,//,T30,36( ' * 1 , IX) ,//,T31 , 'SET' , ,3X,  'MATL' ,4X, 'ALLOY' . 
5/.T31 , 'NUM'  , 3X , 2 * 1 CODE'  .4 X) ,/) 

321  FORMAT (T30,I4,4X,2(A2,6X)) 

C 


c 

C * 
C * 

c 

400 


401 

402 

403 

404 


405 

406 

c 

409 


ALL-TRIAL  VALUE  DATA  AND  TABLE  PRINT-OUT  STATEMENTS:  NPRINT.2  OR  3 

FORMAT  STATEMENTS  400  TO  406  ARE  FOR  PARAMETERS  * 409  TO  413  FOR  STATISTICS 

FORMAT* 1H1 ,T42, 'ANALYTIC  APPROXIMATION  OF  STRESS-STRAIN  PROPERTIES 
1' ,/,T33, 'RAMBERG-OSGOOD  PARAMETERS  OBTAINED  FROM  ANALYSES  OF  EXPER 
2IMENTAL  DATA' ,/) 

F0RMAT(T25, 'VALUES  OF  0.2-PCT  OFFSET-YIELD  AND  R-0', 

1'  EXPONENT  FROM  LOG-LEAST  SQUARES  CALCULATIONS',/) 

FORMAT* T50. ' *****  ELASTIC  PROPERTIES  »**•*' , / , T 1 0 , ' E= ' , F6 . 2 , , 1 X , 

1 ’ MEGA-PS  I ' ,3X, 'S-1NT*' ,F6.0, ' PS  I ' , 3X , 1 E-I NT* 1 , 1 5 , • MU-IN/IN' ,3X, 

2 ' SEE* ' , F6. 2 , ' KSI ' ,3X, 'C-CFT* ' .F7.4.2X, 'NO.  DATA  POINTS*', 13) 

FORMAT(T10, 'E*'  ,F6.1  , ' GPA ' , 8X , ' S-I NT* ' , F6 . 2 , ' MPA ' , 26X , ' SEE* ' . 

1F6.1,’  MPA',/) 

FORMAT ( T51 .' *****  PLASTIC  PARAMETERS  **»**',/, 

1 T25 , ' ANA L I N I T-EPS ' , 4X , ' R-0  EXPNT', 4X, 

2' .2K-Y1ELD  ' ,7X, ' PROP-LIM' , 9X , ' . 1%-YIELD' ,5X, 'ANAL' ,/,T25, 

3 'NO. ' ,2X, 'MU-STRAIN' , 7X , ' N ' , 5X , 3 ( 3X , ' KSI ' , 5X , 'MPA' ,3X) ,1X, 

4'PTS. ' ,/) 

FORMAT (T25,1X,I2,4X,I5,5X,F7.2,3X,3(F7.2,1X,F?.1 ,2X),1X,I3) 

FORMAT (1H+.T23, '**' ,83X,  'BEST-FIT' ) 

F0RMAT(1H0,T49, '**••*  STATISTICAL  PARAMETERS  *»*•*',/, 

1 T5 , 1 8X , ' SPL  TO  MPT* ,12X. 'MPT  TO  SY',13X,'SPL  TO  SY',13X. 

2'SPL  TO  SN'.13X,'SY  TO  SN ' , / , T5 , ' ANA L I N I T ' , 3X , 5 ( 6X , 'SEE' ,5X, ' C-C ' 

3,2X, 'NO  ' ) ,/,T5, IX, 'NO' ,3X, 'EPS' ,3X,5(3X, 'KSI ' ,2X, 'MPA' ,4X, ' R ' ,3X, 
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820 

821 

410 

822 

41  1 

823 

412 

824 

413 

825 

C 

826 

C 

827 

C 

828 

C * 1 

829 

C 

830 

414 

831 

832 

415 

833 

416 

834 

417 

835 

836 

837 

B38 

639 

840 

418 

841 

419 

842 

420 

843 

421 

844 

422 

845 

423 

846 

424 

847 

425 

848 

849 

850 

C 

851 

c 

852 

C 

853 

C • 1 

B54 

c 

855 

426 

856 

857 

858 

859 

427 

860 

661 

862 

428 

863 

429 

864 

430 

865 

866 

867 

868 

431 

869 

432 

870 

433 

871 

872 

434 

873 

435 

874 

875 

436 

876 

877 

C 

878 

c 

879 

C 

880 

C • 

881 

C 

882 

500 

883 

501 

884 

502 

885 

C 

B86 

c 

B87 

c 

888 

4 1 PTS ' ) ,/ ) 

FORMAT (T6,I2,1X,I5,3X,3(1X,F5.2,F5.1,F7.4,I3,1X)) 
F0RMAT(1H+,T5,84X, '••INSUFFICIENT  DATA  BEYOND  SY*.') 

FORMAT  (1H+,T5,79X,F5. 2-,  F5. 1 , F7 .4 , 13 , 3X, '• INSUFFICIENT  DATA* 1 ) 
FORMAT ( 1H+,T5,78X,2(1X,F5.2,F5.1 , F7 . 4 , 13 , 1 X) ) 


NPR I NT  = 3 OPTIONAL  TABLES:  FORMAT  STATEMENTS  414  TO  425 

FORMAT(1H1 ,T7, 'SPEC.  '. A6 , 2X .' POINT-BY-POINT  DEVIATIONS  OF  CALCULA 
1 TED  STRESSES  FROM  OBSERVED  VALUES  USING  THE  PARAMETERS  FOR') 

FORMAT ( 1 H+ .T7.110X, 'J*1  TO  d*5‘,/) 

F0RMAT(1H+,T7, 1 10X, 'Jx6  TO  d«10',/) 

FORMAT (T12.5X. ' ADJ' ,6X, 'OBS' ,2X, 'PLASTIC' , 

15( IX, 14,  IX,  'MU-EPS-I ' , IX), /.Til, 

2 ' PT • ,3X, 'STRAIN' ,2X, 'STRESS' , IX, 'STRAIN' ,1X, 

35(2X, 'STRESS' , 3X , 1 DEV ' , 1 X ) , /T1  1 , 

4 1 NO 1 ,3X, 'IN/IN' ,5X, ' KS I ' ,2X, 'MU-EPS' ,1X, 

55 ( 4X  , 'KSI ' ,4X, ' PC T ' ,1X) ./) 

FORMAT (Tio! I3!2X,F7!6,1X,F7.2,2X,I4,2X,5(1X,F7.2,1XtF5.1 ,1X)) 
FORMAT(T10,I3,2X,F7.6,1X,F7.2,2X,I4) 

FORMAT ( 1H+,T12,27X,F7.2, IX, F5. 1 ) 

FORMAT ( 1H+,T12,42X,F7.2, 1X.F5. 1 ) 

FORMAT ( 1 H+ ,T12,57X,F7.2,1X,F5. 1 ) 

FORMAT ( 1 H+ ,T12,72X,F7.2,1X,F5.1 ) 

FORMAT ( 1 H+ , T 1 2 , 87X , F7 . 2 , 1 X , F5 . 1 ) 

FORMAT(1HO.T4, 124( '-' ),/,T12, '••  BLANK  COLUMNS  ARE  FOR  CASES  WHERE 
1 ANALYTIC  STRESSES  COULD  NOT  BE  CALCULATED.  SEE,  DATA  COMPILATION  F 
20R  REASONS' ,/,T4,124('-' )) 


PI  REPORT  (NPRINT-1)  TABLE:  FORMAT  STATEMENTS  426  TO  433 

FORMATf 1H0.T15. 'SPEC  NO:  ' ,A6,2X, 1 ‘DIGITIZED  STRESS  AND  STRAIN  ', 

1 'VALUES  AND  THEORETICAL  STRESS  AND  STRESS  DEVIATIONS*', 

2/, TIO, ' ANALYTIC  DATA  CALCULATED  USING  BEST-FIT  R AMBERG-OSGOOD 
3 1 PARAMETERS : SY«',F7.2,'  KSI,  N.',F7.2,'  A E«',F6.2,'  M-PSI'./) 
FORMAT ( T1 0, 14X, 'ADd' ,4X, 'PLASTIC' ,4X, 'OBS' ,5X, ' CALC ', /.TIO, ' PT ' ,3X, 
1 ' LOAD' , 4X , 2 ( 'STRAIN' ,2X) ,1X,2( 'STRESS' , 3X ) , 1 X , ' DEV ' , / , T1 0 . ' NO • , 

23X, ' LBS' ,5X, ' IN/IN' ,3X, ' MU-EPS ' , 5X , ' KSI ' , 6X, ' KSI ' , 5X, ' PCT ' , / ) 
F0RMAT(T10, 12, 1X,F7.0,2X,F7.5,2X, I5,3X,2(F7.2,2X) ,F5. 1 ) 

FORMAT (T10.I2,1X,F7.0,2X,F7.5,2X,I5,3X,F7.2) 

FORMAT ( T 1 0 , 2 ( 1 4X , ' ADd' ,4X, 'PLASTIC' ,4X, 'OBS' ,5X, 'CALC' , 16X) ,/,TlO, 
12( 'PT' ,3X, ' LOAD' ,4X,2( 'STRAIN' ,2X) ,1X,2( 'STRESS' , 3X  ) , 1 X , ' pEV  ' . 8X ) 

2 , / , T 1 0 , 2 ( ' NO  1 , 3X , ' LBS ' , 5X , ' IN/ IN ' , 3X , 'MU-EPS' , 5X , ' KSI  ' , 6X ,-' KSI  ' , 
35X, ' PCT' ,8X) ,/) 

FORMAT (T10,2( I2,1X,F7.0,2X,F7.5,2X, I5,3X,2(F7.2,2X) ,F5.1 ,8X) ) 

FORMAT (T10,2(I2,1X,F7.0,2X,F7.5,2X,I5,3X,F7.2,24X) ) 

FORMAT (TIO, 12, 1 X , F7 . 0 , 2X , F7 . 5 , 2X , 1 5 , 3X , 2 ( F7 . 2 , 2X ) ,F5.1 ,8X,I2,1X, 
1F7.0,2X,F7.5,2X,I5,3X,F7.2) 

F0RMAT(T4,124( '-'),/, TIO, 'THERE  WERE  ',12,'  INPUT  DATA  POINTS.') 
F0RMAT(1H+, TIO, 33X , 'STRESS  CALCULATIONS  WERE  MADE  ONLY  ', 

1 ' FOR  DATA  WHOSE  TOTAL  STRAIN  WAS  LESS  THAN  . 02-IN/IN') 
F0RMAT(T4,124('-'),/, TIO, 'STRESS  CALCULATIONS  WERE  MADE  ONLY 
1 ' FOR  DATA  WHOSE  TOTAL  STRAIN  WAS  LESS  THAN  .02-IN/IN') 


PUNCH-OUT  FORMAT  STATEMENTS 

FORMAT (A6,2A2,A4,I2,I1,I6,I5,I2,F6.1 , 3F4 . 1 , F6 . 2 , 2F7 . 2 , F4. 1 , 2 14) 
FORMAT (A6.2A2.A4, 12,11 , A6 , 1 5 , 1 2 , F7 . 2 , 3F4 . 1 , F8 . 2 , 1 8X , 2 14 ) 

FORMAT ( A2 ,1914) 


END 


23 


SUBROUTINE  ORDER ( N , S , E , P . MORDER ) 


1 C 

2 

3 C 

4 C 

SC*  SUBROUTINE  ARRANGES  INPUT  DATA  IN  ORDER  OF  INCREASING  STRAIN  MAGNITUDE 

6 C * USED  ONLY  IF  CHECK  IN  MAIN  PROGRAM  SHOWS  INPUT  DATA  NOT  IN  ORDER 

7 C 

9  C 

10  DIMENSION  S ( 1 ) ,E(1 ) ,P(1  ) 

11  KOUNT  = 0 

12  MsN— 1 

13  1 NOUT  = 0 

14  DO  2 1 = 1 ,M 

15  I F 1 E ( I ) . LE • E ( 1+1 ) ) GO  TO  2 

16  N0UT=N0UT+1 

17  SA  = S I I ) 

18  EA  = E ( I ) 

19  PA=P(I) 

20  SB  = S ( 1 + 1 ) 

21  EB=E( 1+1 ) 

22  PB  = P ( I + 1 ) 

23  S ( I ) =SB 

24  E( I ) =EB 

25  P( I ) =PB 

26  S(I+1)*SA 

27  E ( 1 + 1 ) = EA 

28  P( 1+1 ) = P A 

29  2 CONTINUE 

30  KOUNT=KOUNT+1 

31  IF(KOUNT.LT.N)  GO  TO  3 

32  MORDER=1 

33  RETURN 

34  3 IF(NOUT.GT.O)  GO  TO  1 

35  RETURN 

36  C 

37  END 


1 c 

2 SUBROUTINE  LSTSQ( X , Y , SL , OF , LI , Ld ) 

4 C 

.SC*  SUBROUTINE  PERFORMS  A LINEAR  REGRESSION  ANALYSIS  (LEAST  SQUARES) 

6 C * FOR  MODULUS  AND  INTERCEPT:  X=STRAIN  Y=STRESS 

7 C * FOR  RAM3ERG-0SG00D  PARAMETERS:  X = LOG  PLASTIC-STRAIN  Y* LOG  STRESS 

8 C 


10 

C 

1 1 

DIMENSION  X( 1 ) ,Y( 1 ) 

12 

ET  = 0 

13 

ST*0 

14 

EET  = 0 

15 

EST*0 

16 

DO  1 1=LI,LJ 

17 

ET=ET+X( I) 

18 

ST*ST+Y( I) 

19 

EET*EET+X( I )**2 

20 

1 

EST*  EST+X ( I ) * Y ( I ) 

21 

NC*LJ-LI+1 

22 

DEN=NC*  EET-ET  * *2 

23 

ENUM=NC*EST-ET*ST 

24 

ONUM  = ST  *EET-ET*EST 

25 

SL*ENUM/DEN 

26 

OF  = ONUM/DEN 

27 

RETURN 

28 

C 

29 

ENO 
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1 c 

2 SUBROUTINE  R JA75( Z . XP , BB , AA , I TER ) 

3 C 

4 C 

5 C • SUBROUTINE  CALCULATES  NON-D  R-0  STRESS  FOR  A GIVEN  VALUE  OF  STRAIN  * 

6 C • NUMERICAL  SOLUTION  OF  FUNC(Z)  WHERE  Z*S/SY 

7 C 

8 C • IF  NO  SOLUTION  IN  100  ITERATIONS,  CALCULATIONS  ABORTED:  ITER-100 

9 C 

11  C 

12  FUNC ( Z ) ■Z+AA»Z**XP-BB 

13  I TER  * 1 

14  IF(XP. LT.1 . ) GO  TO  6 

15  RM-1.E-5 

16  DZ«.1 

17  RZ*FUNC(Z) 

18  IF(RZ.GT.O. . AN0.XP.GE.200. ) DZ-.01 

19  IF(RZ. LT.O. -AND. XP.GE. 5000. ) GO  TO  6 

20  IF(RZ.LT.O. .AND. XP.LT. 5000.. ANO. XP.GE. 1000.)  DZ«.0001 

21  IF(RZ. LT.O. .AND. XP. LT. 1000. .AND. XP.GE. 500. ) DZ*.001 

22  IFIRZ.LT.O. .AND. XP.LT. 500. .AND. XP.GE. 100.)  D2-.01 

23  1 IF(ABS(RZ).LE.RM)  GO  TO  6 

24  ITER*ITER+1 

25  IF(ITER.EQ.IOO)  GO  TO  6 

26  I F ( RZ  . LT.O.  )G0  TO  3 

27  I F ( 2 . LE.  .01 ) GO  TO  8 

28  I F ( Z . LE. . 1 . ANO.Z.GT.O. ) GO  TO  2 

29  Z*Z-DZ 

30  I F ( 2 . LE. .01 ) GO  TO  B 

31  TEST=XP.AL0G10(Z) 

32  IF(TEST. LE.-20. ) GO  TO  9 

33  IF(TEST.GE.20. ) GO  TO  7 

34  RZ*FUNC(Z) 

35  GO  TO  1 

36  2 DZ* . 1 * DZ 

37  Z.Z-DZ 

38  RZ*FUNC(Z) 

39  GO  TO  1 

40  3 V.Z+DZ 

41  TEST*XP*AL0G10(V) 

42  I F ( TEST . LE . “20 . ) GO  TO  9 

43  IF(TEST.GE.20. ) GO  TO  7 

44  RV«FUNC(V) 

45  IF(ABS(RV).LE.RM)  GO  TO  5 

46  I TER* I TER+1 

47  IF(ITER.EQ.IOO)  GO  TO  6 

48  I F ( RV . LT . 0 . ) GO  TO  4 

49  DZ* . 1 *DZ 

50  I F ( V . GE . 1 • . AND . DZ • LT • RM ) GO  TO  5 

51  GO  TO  3 

52  4 Z*Z+DZ 

53  GO  TO  3 

54  5 Z.V 

55  6 RETURN 

56  7 I TER*200 

57  RETURN 

58  8 ITER- 300 

59  RETURN 

60  9 Z*BB 

61  RETURN 

62  C 

63  C 

64  C 

65  C • I T ER*200  • ABORT:  Z-VALUE  INCREASING  ANO  FUNC(Z)  DIVERGING 

66  C * TESTED  IN  LINES  33  AND  43 

67  C 

68  C * ITER*300  * ABORT:  Z-VALUE  BECOMING  VANISHINGLY  SMALL 

69  C * TESTEO  IN  LINES  27  AND  30 

70  C 

72  C 

73  END 
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1 PASG , A R»EVALRO. 

2 PASG.A  R*  LOCATER . 

3 PUSE  1 0 , R*  LOCATER . 

4 PXQT  R*EVALRO. FIT DATA 


1 PPREP  R*EVALRO. 

2 PMAP.I  CIS,R«EVALRO. FITDATA 

3 IN  R» EVALRO . MAIN 

4 LIB  R»EVALRO. 

5 END 

6 PPACK  R*EVALRO. 

7 PPRT.TL  R*EVALRO. 
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LISTING  OF  R*RETRIEVE 


The  RETRIEVE  file  consists  of  the  following  file  elements.  Unless  otherwise 
noted,  the  subroutines  are  used  with  both  main  programs. 

Element  Function 

PMAIN  Main  program  for  data  retrieval  of  stored  data 

on  teletype  terminals 

TMAIN  Main  program  for  data  retrieval  of  stored  data 

on  cathode-ray-tube  graphics  terminals 

TYCODE  Subroutine  of  PMAIN  for  displaying  material  and 

alloy  code  tables  and  for  input  of  users 
material/alloy  selection 

CODE  Subroutine  of  TMAIN  with  same  function  as  TYCODE 

RANGE  Finds  and  prints  extreme  values  of  selected  re- 

trieval option  parameter  in  the  available  data 
for  the  selected  material/alloy 

COUNT  Counts  and  prints  the  number  of  data  sets  whose 

value  of  the  selected  retrieval  option  parameter 
is  within  the  user-specified  limits 

LONG  Subroutine  of  TMAIN  which  retrieves  desired  data 

from  data  file  and  prints  a long  form  list 

LPRT  Subroutine  of  PMAIN  with  same  function  as  LONG 

for  132-character  line-length  teletype  terminals 

SHORT  Subroutine  of  TMAIN  which  retrieves  desired  data 

from  data  file  and  prints  a short  form  list 

SPRT  Subroutine  of  PMAIN  with  same  function  as  SHORT 

for  132-character  line-length  teletype  terminals 

PORPRT  Subroutine  of  PMAIN  which  retrieves  desired  data 

from  data  file  and  prints  an  expanded  short 
form  list  on  80-character  line-length  teletype 
terminals 

SYSIZE  Subroutine  of  TMAIN  examines  yield  stress  data  for 

sets  selected  for  graphical  display  and  deter- 
mines which  of  the  sets  has  the  highest  yield 
stress  value 

CALC  Subroutine  of  TMAIN  which  calculates  the  plotting 

points  for  a graphical  stress-strain  curve 
display 

BUFFPK  Used  in  forming  the  absolute  of  TMAIN  with  the 

same  function  as  in  PICPTS  (see  page  2). 
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There  are  four  elements  used  for  control  functions  whose  listings  are  also 
given.  These  are: 


Control 

Element  Function 


TELTYP 

CRTUBE 

MAPTYP 

MAPCRT 


Assigns  DATABANK  and 
executes  PMAIN 

LOCATER 

files 

to 

the 

run 

and 

Assigns  DATABANK  and 
executes  TMAIN 

LOCATER 

files 

to 

the 

run 

and 

After  edit  of  PMAIN, 
the  absolute  (name 

used  to 
: TYP) 

PREP 

file 

and 

to 

form 

After  edit  of  TMAIN, 
the  absolute  (name 

used  to 
: CRT) . 

PREP 

file 

and 

to 

form 

Element  TMAIN  uses  the  following  TEKTRONIX  PLOT-10  subroutines:  INITT, 
CHRSIZ,  ANMODE,  NEWPAG,  IOWAIT,  BINITT,  NPTS,  CHECK,  DSP LAY,  LINE,  SYMBL,  SIZES, 
CPLOT,  MOVABS,  and  FINITT.  No  listings  have  been  given  for  these  subroutines. 
The  PLOT-10  subroutines  are  referenced  from  the  file  DAO*TEX  and  from  the  BUFFPK 
package  when  the  absolute  for  TMAIN  is  created. 


RETRIEVE  may  be  executed  interactively  on  any  of  three  types  of  terminals: 
TEKTRONIX  4014  graphics  terminals,  132-character  carriage-length  teletype  ter- 
minals, or  80-character  carriage- length  teletype  terminals.  Data  are  retrieved 
by  the  program  from  the  R*DATABANK  file  in  the  computer  using  the  R*LOCATER  file 
as  described  in  Reference  2. 


Listings  of  the  two  main  programs  and  the  subroutines 
29-53.  The  four  control  elements  are  listed,  in  the  order 
54. 


are  given  on  pages 
given  above,  on  page 


28 


1 c 

2 C * FILE  R'RETRIEVE  IS  FOR  RETRIEVING  DATA  FROM  DATABANK 

3 C * THIS  IS  THE  ELEMENT  PMAIN  (ABSOLUTE:  TYP) 

5 C 

6 C « THIS  ELEMENT  IS  FOR  KEYBOARD  PRINT  TERMINALS 

7 C * 132-CHARACTER  LINE-LENGTH  TERMINALS  HERE  CALLED  TELETYPE 

8 C * 80-CHARACTER  LINE-LENGTH  TERMINALS  HERE  CALLED  PORTABLE 

9 C * ANALYZED  STRESS-STRAIN  DATA  LISTED  FOR  REQUESTED  MATERIAL/ALLOY 
10  C * THIS  IS  AN  INTERACTIVE  PROGRAM  REQUIRING  RESPONSES  FROM  USER 

1 1 C 

12  C * PROGRAM  EXECUTED  WITH  OADD  R*RETRIEVE.TELTYP 

13  C * FILE  20  IS  DATABANK 

14  C 

15  DEFINE  FILE  20 ( 9999 , 80 , E . 1 1 ) 

16  C 

17  C 

18  C 

19  C * IDENTIFY  TYPE  AND  PREPARE  TERMINAL:  EXPAND  BUFFER  TO  FULL  CAPACITY 

20  C 

21  1 WR I TE ( 6 , 300 ) 

22  REA0(5, 201 ) .TERM 

23  I F ( TERM . EQ . 'TEL' ) WRITE(6,301) 

24  I F ( TERM . EQ . 1 POR 1 ) WRITE(6,302) 

25  IF(TERM.EQ. 'TEL' .OR.TERM.EQ. ‘POR‘ ) GO  TO  10 

26  WR I TE ( 6 , 1 02 ) 

27  READ( 5,201). ANS 

28  IF(ANS.EQ. ' YES' ) GO  TO  1 

29  GO  TO  50 

30  10  READ( 5 ,200 ) 

31  C 

32  

33  C 

34  C * READ-IN  OF  MATERIAL/ALLOY  CODE  OF  OATA  TO  BE  LISTED 

35  C 

36  11  CALL  TYCODEf JEND.NI ,W,X) 

37  C 

38  

39  C FOUR  OPTIONS  OF  SFLECIVE  OF  DATA  RETRIEVAL  FOR  CHOSEN  MATL/ALLOY 

40  C I DENT  = 0 : ALL  DATA  SETS 

41  C I DENT  * 1 : ONLY  DATA  WHOSE  ULTIMATE  STRENGTH  IS  WITHIN  A SELECTED  RANGE 

42  C I DENT=2 : ONLY  DATA  WHOSE  YIELD  STRENGTH  IS  WITHIN  A SELECTED  RANGE 

43  C IDENT.3:  ONLY  DATA  WHOSE  ROCKWELL-C  HARDNESS  IS  WITHIN  A SELECTED  RANGE 

44  C 

45  C • TWO  TELETYPE  PRINT-OUT  OPTIONS:  LONG-FORM  OR  SHORT-FORM 

46  C * ONLY  EXPANDED  SHORT-FORM  PRINT-OUT  AVAILABLE  FOR  PORTABLE  TERMINALS 

47  C 

48  IF(JEND.EQ.O)  GO  TO  50 

49  IHOLD-NI 

50  C 

51  

52  C 

53  C • IDENTIFICATION  OF  SELECTIVE  RETRIEVAL  OPTION 

54  C 

55  12  IDENT.O 

56  XMAX*0 

57  XMIN=0 

58  WRITE(6,101 ) 

59  READ(5.201 ) .ANS 

60  IF(ANS.EQ. ' ALL' ) GO  TO  18 

61  I F ( ANS . EQ . ' ULT 1 ) I DENT  * 1 

62  1 F( ANS . EQ. 1 Y I E 1 ) I DENT  <2 

63  I F ( ANS . EQ . 1 HAR 1 ) I DENT*3 

64  I F ( 1DENT.GT.0)  GO  TO  13 

65  WR I TE ( 6 , 1 02 ) 

66  REAO( 5,201), ANS 

67  I F( ANS. NE. 1 YES  1 ) GO  TO  50 

68  GO  TO  12 


} 


i 


«■ 
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69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
1 00 


ENTRY  OF  PARAMETER  VALUE  AND  RANGE 

SUBROUTINE  RANGE  FINOS  ANO  PRINTS  EXTREMES  OF  PARAMETER  VALUES 


13 


14 


1 5 


16 


C 

C— 

C 

c * 
c 

17 


CALL  RANGE(W,X,NI , IDENT ,JENO) 

NI = IHOLD 

IF(dEND.EQ.O)  GO  TO  50 
IF( IDENT. EO.O)  GO  TO  18 
WRITE (6, 103) 

READI5.200,ERR*30) , XMIN, XMAX 

IF(XMAX.GT.XMIN)  GO  TO  15 

YMIN=XMAX 

YMAX=XM1N 

XMAX  = YMAX 

XMIN=YMIN 

READ! 20' NI ,202) , NT YPE , SU LT , RCH , S Y , NEXT 
IFINTYPE.NE. 1 ) GO  TO  16 

IF( IDENT. EO- 1 .AND. SULT. LE. XMAX. AND. SULT.GE.XMIN)  GO  TO  18 
IF( IDENT. EQ. 2. AND. SY.LE. XMAX. AND. SY.GE. XMIN)  GO  TO  18 
IF( IDENT . EO. 3. AND. RCH. LE .XMAX .AND. RCH. GE.XMIN)  GO  TO  18 
IF(NEXT.EQ-O)  GO  TO  17 
NI =NEXT 
GO  TO  15 


CORRECTION  OPTIONS  WHEN  NO  DATA  IN  DESIRED  RANGE  OF  SELECTED  PARAMETER 

NI = IHOLD 
WRITE(6, 105) 


101 

READ  ( 5 , 

201  ) 

, ANS 

102 

I F ( AMS. 

EO.  ' 

RAN  ' ) 

GO 

TO 

14 

103 

I F ( ANS  . 

£0.  1 

OPT  1 ) 

GO 

TO 

12 

104 

I F ( ANS  . 

EO.  1 

NEW'  ) 

GO 

TO 

1 1 

105 

I DENT -0 

106 

IF( ANS. 

EO.  ' 

ALL'  ) 

GO 

TO 

18 

107 

108 
109 
1 10 
1 1 1 
112 

113 

114 

115 

116 
117 
1 18 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 


GO  TO  50 


PRINT-OUT  OPTIONS:  LONG-FORM  OR  SHORT-FORM  FOR  STANDARD  TELETYPE 
ONLY  EXPANDED  SHORT-FORM  AVAILABLE  FOR  PORTABLE  TERMINAL 

18  NI  = I HOLD 

IF(IDENT.GT.O)  CALL  COUNT ( IDENT , NI , XMIN , XMAX ) 

19  1 F ( TERM . EO • ' POR 1 ) WRITE(6,104) 

IF(TERM.NE. 1 POR 1 ) WRITE(6,106) 

READI 5,201 ) , ANS 

IFIANS.EO. 'SHO' ) GO  TO  20 
I F ( ANS . EO . ' LON' ) GO  TO  21 
I F ( ANS . £0 . ' PR1 1 ) GO  TO  22 
I F ( ANS . EO • 'CHA' ) GO  TO  1 2 
I F ( ANS . EO . 'END' ) GO  TO  50 
WRITEI6, 107) 

GO  TO  19 


C 

c— 

c 

c * 
c 

20 

21 

22 


PRINT-OUT  OF  LISTING 

CALL  SPRT(W,X,NI , IDENT, XMIN, XMAX) 

GO  TO  40 

CALL  LPRT ( W , X , N I , I DENT , XMI N , XMAX ) 

GO  TO  40 

CALL  PORPRT(W,X,NI , I DENT , XMIN . XMAX ) 
GO  TO  40 
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I 


137  C 

138  

139  C 

140  C * ERROR  DIAGNOSTIC  AND  OPTION  TO  CORRECT  WHEN  MIN-MAX  ENTRY  ERROR 

141  C 

142  30  WRITE(6,107) 


143  WR I TE ( 6 , 400  ) 

144  READ ( 5 , 20 1 ) > ANS 

145  IF(ANS.EQ. ' YES' ) GO  TO  14  | 

146  GO  TO  50 

147  C 

149  C 

150  C * OPTIONS:  NEW  LISTING  OR  TERMINATE  l 

151  C 

152  40  WR I TE ( 6 , 1 08 ) 

153  READ( 5 , 201 ) , ANS  I 

154  I F ( ANS ■ EO. ' NEW 1 ) GO  TO  1 1 


155  C 

156  C 

157  C 

158  C • TERMINATION 

159  C ' 

160  50  WRITE(6, 109)  ‘ 

161  51  STOP 

162  C I 

1 63  C 

164  C 

165  101  FORMAT(T2, 'TABLE  BELOW  LISTS  RETRIEVAL  OPTIONS  AND  ENTRIES',//, 

166  1T2, 'RETRIEVAL  OPT  I ON ', 1 4X ,' ENTRY ',/, T2 , 35( '-'),/ , 

167  2T2 , 'ALL  DATA  FOR  SELECTED  MAT L ' , 4X , ‘ ALL ' , / , 

168  3T2 , ' DAT  A FOR  SPECIFIC  UTS  RANGE  1 , 3X ULT IMATE ',/ , 

169  4T2 , 'DATA  FOR  SPECIFIC  SY  RANGE ‘ , 4X , ’ Y I ELD 1 , / , ' 

170  5T2 , ' DATA  FOR  SPECIFIC  RCH  RANGE 3X ,' HARDNESS ',/ , 

171  6T2 , ' MAKE  YOUR  ENTRY  NOW',/) 

172  102  FORMAT(T2. ' ENTRY  DID  NOT  CORRESPOND  TO  EITHER  OF  THE  ABOVE',/, 

173  1 T2 , ' I F YOU  WISH  TO  CORRECT,  ENTER  YES',/,  - 

174  2T2, 'OTHERWISE  ENTER  END  AND  RUN  WILL  TERMINATE') 

175  103  FORMAT(/,T2, 'DATA  WILL  BE  RETRIEVED  FOR  A RANGE  OF  VALUES  OF  ', 

176  1 ' THE  CHOSEN  PARAMETER ',//, T2 ,' NEXT  ENTRY  DEFINES  LOWER  ANO  ', 

177  2 ' UPPER  BOUNDS  OF  DESIRED  RANGE', 

178  3//, T2. 'ENTRY  UNITS:  KSI  UNITS  FOR  STRENGTH  AND  ', 

179  4 ' ROCKWE LL-C  FOR  HARDNESS ‘ ,/, T2 ,' SEPARATE  THE  TWO  VALUES  BY  A • , 

180  5 'COMMA.  EXAMPLE:  1 50 ., 1 55 .',//, T2 ,' ENTER  BOUND  VALUES  NOW') 

181  104  FORMAT! /.T2 ,' PRINT  NOW,  CHANGE  RETRIEVAL  OPTION,  OR  TERMINATE', 

182  1/,T2,4B(  '-'  ),/,T2,  ' 1 . TO  PRINT  A DAT A-LI ST  ENTER:  PRINT', /,T2, 

183  2 ' 2 . TO  CHANGE  RETRIEVAL  OPTION  ENTER:  CHANGE ',/. T2 , 

184  3 ‘ 3 • TO  TERMINATE  THE  PROGRAM  ENTER:  END',/,T2, 

185  4 'MAKE  YOUR  ENTRY  NOW') 

186  105  FORMAT(T2, 'THERE  ARE  NO  DATA  IN  THE  DATA  BANK  FOR  THE'/,T2,, 

187  1 'SPECIFIED  RANGE  OF  THE  SELECTED  PARAMETER ',/,  T2 , ; 

168  3 ' TO  CHANGE  THE  RANGE  LIMITS,  ENTER:  RANGE',/,  v 

189  4T2 , ' TO  LIST  ALL  THE  DATA  FOR  THE  MATERIAL,  ENTER:  ALL',/, 

190  5T2 , 'TO  CHOOSE  ANOTHER  OPTION,  ENTER:  OPTION',/, 

191  6T2 , ' TO  CHOOSE  ANOTHER  MATERIAL.  ENTER  NEW',/,  J 

192  7T2, 'OTHERWISE  ENTER  END  AND  RUN  WILL  TERMINATE') 

193  106  F0RMAT(T2, 'PRINT  NOW,  CHANGE  RETRIEVAL  OPTION,  OR  TERMINATE ',/, T2 , 

194  1 48 ('  — '),/, T2 ,' 1 . TO  PRINT  A LONG-FORM  DA T A-LI ST  ENTER:  LONG', 

195  2/ , T2 , ' 2 . TO  PRINT  A SHORT-FORM  DATA-LIST  ENTER:  SHORT' ,/,T2, 

196  3 ' 3 • TO  CHANGE  RETRIEVAL  OPTION  ENTER:  CHANGE ',/, T2 , 

197  4 ' 4 . TO  TERMINATE  THE  PROGRAM  ENTER:  END',/,T2, 

198  5 'MAKE  YOUR  ENTRY  NOW) 

199  107  FORMAT!/, T2, '•*  ERROR:  ENTRY  NOT  RECOGNIZED  **') 

200  108  FORMAT(T2,2X, 'OPTIONS:  MAKE  ONE  OF  TWO  ENTRIES' ,/,T2,2X,31( ' “' ), 

201  1 / , T2 , ' FOR  A DATA  LIST  OF  ANOTHER  MATERIAL  ENTER:  NEW,/,T2, 

202  2 ' TO  TERMINATE  THE  PROGRAM  ENTER:  END') 

203  109  FORMAT (T2, 'RUN  IS  OVER.  WHEN  NEXT  PROMPT  (>)  APPEARS  ENTER  *FIN') 

204  C 
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205 

200 

FORMAT ( ) 

206 

201 

FORMAT ( A3 ) 

207 

202 

FORMAT (16X. I 1 , 13X, F6. 1 ,8X, F4. 1 , 

,13X,F7.2,4X, 

.14) 

208 

C 

209 

300 

FORMAT (T2, 'THIS  PROGRAM  CAN  BE 

RUN  ONLY  ON 

HARD-PRINT  TERMINALS 

210  2/. T2, ' IDENTIFY  THE  TERMINAL  YOU  USING: 1 ,/, T2 , 

211  3' IF  ON  A 132-CHARACTER  TERMINAL  ENTER:  T ELETYPE ' , / , T2 , 

212  4' IF  ON  AN  80-CHARACTER  TERMINAL  ENTER:  PORTABLE  1 T2 , 

213  5 ' MAKE  ENTRY  NOW  AND  KEY  RETURN') 

214  301  F0RMAT(T2, ix, 'PREPARE  TERMINAL  FOR  DATA  RETR I EVAL T2 , 

215  1 ' 1 . ENTER:  PPTTY  W.132  THEN  KEY  RETURN T2 , 

216  2 ' 2 • WAIT  FOR  RESPONSE  FROM  COMPUTER!  -••COMPLETE' ,/,T2, 

217  3 ' 3 . THEN  WHEN  PROMPT  APPEARS  KEY  RETURN') 

218  302  F0RMAT(T2, IX, 'PREPARE  TERMINAL  FOR  DATA  RETR I EVA L !',//, T2 , 

219  1*1.  ENTER:  PPTTY  W,80  THEN  KEY  RETURN ',/, T2 , 

220  2 ' 2 WAIT  FOR  RESPONSE  FROM  COMPUTER:  -P® COMPLETE ' ,/ , T2 , 

221  3 ' 3 ■ THEN  WHEN  PROMPT  APPEARS  KEY  RETURN') 

222  C 

223  400  FORMAT ( T2 , ‘ I F YOU  WISH  TO  CORRECT,  ENTER  YES',/, 

224  1T2, 'OTHERWISE  ENTER  END  AND  RUN  WILL  TERMINATE') 

225  C 

226  END 


1 C 

2 C * FILE  R’RETRIEVE  IS  FOR  RETRIEVING  DATA  FROM , DAT A8ANK 

3 C * ELEMENT  TMAIN  (ABSOLUTE:  CRT)  FOR  TEKTRONIX  GRAPHICS  TERMINAL 

4 C 

5 C 

6 C * ANALYZED  STRESS-STRAIN  DATA  LISTED  FOR  REQUESTED  MATERIAL/ALLOY 

7 C * OPTIONAL  GRAPHICS  ROUTINE  PLOTS  STRESS-STRAIN  CURVE  FOR  ANY  SET 
SC*  THIS  IS  AN  INTERACTIVE  PROGRAM  REQUIRING  RESPONSES  FROM  USER 

9 C 

IOC*  PROGRAM  EXECUTED  WITH  PADD  R*RETR I EVE . CRTUBE 

1 1 C 

12  DIMENSION  LOC ( 1 9 ) , S ( 64 ) , E ( 64 ) , EP ( 64 ) , SS( 3 ) , EE ( 3 ) 

13  C 

15  C 

16  C * FILE  10  IS  LOCATER  AND  FILE  20  IS  DATABANK 

17  C 

18  DEFINE  FILE  10(50, 80, E.JJ) 

19  DEFINE  FILE  20 ( 9999 , 80 , E , 1 1 ) 

20  C 

22  C 

23  C * PREPARE  TERMINAL:  EXPAND  BUFFER  TO  132  CHARACTERS 

24  C 

25  CALL  INITT( 30 ) 

26  CALL  TERM( 2 , 1 024 ) 

27  CALL  CHRSIZ<4) 

28  CALL  ANMODE 

29  KURVE*0 

30  WR I TE ( 6 , 100) 

31  REAO ( 5 , 200 ) 

32  C 

33  C 

34  C 

35  C • REAO-IN  OF  MATERIAL/ALLOY  CODE  OF  DATA  TO  BE  LISTED 

36  C 

37  1 CALL  NEWPAG 

38  CALL  IOWAIT ( 20 ) 

39  CALL  CODE ( LOC , JEND , N I , W , X ) 

40  C 

41  C 

42  C • FOUR  OPTIONS  OF  SELECIVE  OF  DATA  RETRIEVAL  FOR  CHOSEN  MATL/ALLOY 
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43  C * I DENT  = 0 : ALL  DATA  SETS 

44  C * I DENT* 1 : only  DATA  WHOSE  ULTIMATE  STRENGTH  IS  WITHIN  A SELECTED  RANGE 

45  C * IDENT  = 2 : ONLY  DATA  WHOSE  YIELD  STRENGTH  IS  WITHIN  A SELECTED  RANGE 

46  C • I DENT  = 3 : ONLY  DATA  WHOSE  ROCKWELL-C  HARDNESS  IS  WITHIN  A SELECTED  RANGE 

47  C 

48  C * TWO  PRINT-OUT  OPTIONS:  LONG-FORM  PRINT-OUT  OR  SHORT-FORM  PRINT-OUT 

49  C 

50  I F ( U END . EQ . 0 ) GO  TO  50 

51  IH0LD*NI 

52  C 

53  C 

54  C 

55  C • IDENTIFICATION  OF  SELECTIVE  RETRIEVAL  OPTION 

56  C 

57  2 IDENT  = 0 

58  XMAX  = 0 

59  X M I N = 0 

60  WRITE! 6, 101  ) 

61  RE AD ( 5 , 20 1 ) , ANS 

62  IFfANS.EO. 1 ALL' ) GO  TO  8 

63  IFfANS.EO.  'ULT'  ) I DENT  = 1 

64  IFfANS.EO. ’ Y I E * ) I DENT  = 2 

65  IFfANS.EO. 'HAR' ) IDENT  *3 

66  IFf IDENT. GT.O)  GO  TO  3 

67  WR I T E ( 6 , 102) 

63  READ(5,201 ) .ANS 

69  IFfANS.NE. 'YES' ) GO  TO  50 

70  GO  TO  2 

71  C 

72  C 

73  C 

74  C • ENTRY  OF  PARAMETER  VALUE  AND  RANGE 

75  C • SUBROUTINE  RANGE  FINDS  AND  PRINTS  EXTREMES  OF  PARAMETER  VALUES 

76  C 

77  3 CALL  RANGEfW.X.NI . IDENT, UEND) 

78  NI  = I HO  LD 

79  IFfdEND. EQ.O)  GO  TO  50 

80  IFf IDENT. EO.O)  GO  TO  8 

81  4 WR I TE ( 6 , 1 03 ) 

82  MX= 1 

83  READI5, 200, ERR=45) ,XMIN, XMAX 

84  IF(XMAX.GT.XMIN)  GO  TO  5 

85  YMIN=XMAX 

86  YMAX  = XMI N 

87  XMAX=YMAX 

88  XM I N=  YM I N 

89  5 READ(20'NI, 202) , NTYPE , SULT , RCH , SY, NEXT 

90  IFfNTYPE.NE. 1 ) GO  TO  6 

91  IFf IDcNT.EQ. 1 .AND. SULT. LE. XMAX. AND. SULT. GE.XMIN)  GO  TO  8 

92  IFf IDENT. ED. 2. AND. SY.LE. XMAX. AND. SY. GE.XMIN)  GO  TO  8 

93  IFfIDENT. EQ • 3. AND. RCH. LE. XMAX. AND. RCH. GE.XMIN)  GO  TO  8 

94  6 I F ( NEXT . EO ■ 0 ) GO  TO  7 

95  N I =NE  XT 

96  GO  TO  5 

97  C 

98  C 

99  C 

100  C • CORRECTION  OPTIONS  WHEN  NO  DATA  IN  DESIRED  RANGE  OF  SELECTED  PARAMETER 

101  C 

102  7 NI * I HOLD 

103  WR I T E ( 6 , 1 05 ) 

104  RE AD (5, 201 ) ,ANS 

105  IFfANS.EO. 'RAN' ) GO  TO  4 

106  IFfANS.EO. 'OPT' ) GO  TO  2 

107  IFf ANS. EO. 'NEW' ) GO  TO  1 

108  IDENT  *0 

109  IFfANS.EO. 'ALL' ) GO  TO  8 

110  GO  TO  50 
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s 

1 1 1 

c 

113 

c 

114 

c * 

PRINT-OUT  OPTIONS;  LONG-FORM  OR  SHORT-FORM 

115 

c 

L 

116 

8 

NI * IHOLD 

117 

IF(IDENT.GT.C)  CALL  COUNT ( IDENT , NI , XMIN , XMAX ) 

118 

9 

WRITE16, 106) 

119 

READ! 5 , 201 ) . ANS 

120 

I F ( ANS . EO . ' SHO 1 ) GO  TO  10 

121 

IFIANS.EO. ' LON1 ) GO  TO  11 

122 

I F ( ANS . EO . 'CHA' ) GO  TO  2 

123 

I F ( ANS . EO • 'END' ) GO  TO  50 

124 

WRITE(6, 107) 

125 

GO  TO  9 

126 

C 

127 

C * 

PRINT-OUT  OF  LISTING 

128 

C 

1 29 

10 

CALL  SHORT ( W , X , N I , I DENT , XMIN , XMAX ) 

130 

GO  TO  20 

131 

1 1 

CALL  LONG(W,X,NI , I DENT , XMIN , XMAX ) 

132 

c 

133 

c 

134 

c » 

OPTIONS:  NEW  LISTING  * PLOT  OF  S-E  CURVE  » TERMINATE 

135 

c 

1 35 

20 

CALL  NEWPAG 

137 

CALL  I OWA I T ( 20 ) 

138 

WRITE (6, 1 18) 

139 

READ (5,201 ) .ANS 

140 

I F ( ANS . EO . 'NEW1 ) GO  TO  1 

141 

IFIANS.EO. 'CUR1 ) GO  TO  30 

142 

GO  TO  50 

f 

143 

C 

► 

145 

C 

146 

C * 

GRAPHIC  DISPLAY  OF  STRESS  STRAIN  CURVE 

147 

C 

148 

c • 

PREPARATION  OF  TERMINAL  FOR  PLOTTING:  REDUCE  BUFFER  TO  72  CHARACTERS 

149 

c 

150 

c * 

ASSIGNMENT  OF  PLASTIC  STRAIN  VALUES  USED  IN  COMPUTATIONS 

151 

c * 

STATEMENTS  30  TO  33  DONE  ONCE  IN  RUN  (KURVE-O) 

152 

c 

* 

153 

30 

WRI TE(6,1 19) 

154 

READ! 5.201 ) , ANS 

155 

IF(KURVE.GT.O)  GO  TO  34 

156 

KURVE*KURVE+1 

157 

NC  = 64 

158 

EP( 1 )*0 

159 

EP( 2 ) « 1 . E-5 

160 

DO  31  1*3, NC 

161 

I F ( I . LE . 5 ) EP( I )«( 1-2) *5 . E-5 

162 

I F ( I . GT . 5 ) EP ( I ) *2 . E-4+ ( I -6 ) * 1 . E-4 

163 

31 

CONTINUE 

164 

C 

166 

Vrf”’ 

c 

167 

c * 

IDENTIFICATION  OF  CURRENT  END-OF-FILE  SET  NUMBER  OF  DATABANK 

163 

c 

169 

IENDcO 

170 

DO  33  d-1 ,7 

171 

READ* 1 0 1 J , 203 ) ,AA, ( LOC( L ) , L» 1 , 1 9 ) 

l 

172 

DO  32  L*1 ,19 

173 

I F ( LOC ( L ) . LT . I END ) GO  TO  32 

174 

I END*  LOC ( L ) 

175 

32 

CONTINUE 

176 

33 

continue 

177 

1 78 

c 

34 

i . . - 

179 

C 

180 

C * 

181 

C * 

182 

C 

183 

34 

184 

185 

186 

187 

188 

189 

190 

191 

35 

1 92 

193 

1 94 

195 

196 

197 

198 

1 99 

200 

201 

202 

203 

2C4 

C 

205 

C— 

206 

C 

207 

C * 

2C8 

C * 

209 

c * 

210 

c 

211 

36 

212 

213 

214 

215 

216 

217 

218 

219 

220 

221 

222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

C 

232 

C • 

233 

C 

234 

235 

37 

236 

237 

238 

239 

240 

241 

242 

243 

38 

244 

245 

246 

39 

READ-IN  OF  SETS  TO  BE  GRAPHED  AND  CHECK  THAT  SET  NOS.  NOT  BEYOND  EOF 
THREE  CURVES  ON  ONE  GRAPH  PERMITTED 

WRITE(6, 120) 

MX  = 2 

READ! 5.200, ERR=45) , IANS 

IF( IANS.GT.O. AND. IANS. LE.3)  GO  TO  35 

WRITE(6, 121 ) 

RE AD (5. 201 ) . ANS 

I F ( ANS . EQ . 'YES' ) GO  TO  34 

GO  TO  51 

I SET  1 =0 

I SET2=0 

I SET3=0 

WRITE(6, 122) 

MX  = 3 

IF( IANS.EQ. 1 ) READ(5,200,ERR=45) , ISETl 

IF( IANS.E0.2)  READ(5,200.ERR=45) , ISETl , ISET2 

IF( IANS. EO-3)  READ(5.200.ERR=45) , ISETl , ISET2, ISET3 

IF( IEND.GE. ISETl . AND. IEND.GE. ISET2. AND. IEND.GE. ISET3)  GO  TO  36 

WR I TE ( 6 , 123) . I END 

RE AD ( 5 , 20 1 ) .ANS 

I F ( ANS . EO . 'YES' ) GO  TO  34 

GO  TO  51 


COMPUTATIONS  AND  DISPLAY  OF  CURVE  NUMBER  1 

CURVE  NUMBER  1 IS  ALWAYS  THE  DATA-SET  WITH  MAXIMUM  YIELD  STRESS 
SUBROUTINE  SYSIZE  ASSIGNS  ISETl  TO  SET  WITH  MAX  SY 

S( 1 ) =0 
E ( 1 ) =0 

CALL  SYS  I ZE( IANS. ISETl , ISET2. ISET3) 

I P= I SETl 

READ( 20 ' IP, 204) ,W1 , X 1 , EMOD, XN , S Y 
EMAX=SY/(EM0D*1 .E3)+.004 

CALL  CALC(NC,EMOD,XN,SY,S,E,EP.SS,EE,EMAX,NP) 

CALL  ERASE 
CALL  CHRS I Z ( 3 ) 

CALL  BINITT 
CALL  NPTS(NP) 

CALL  CHECK(E.S) 

CALL  DSPLAY ( E , S ) 

CALL  NPTS ( 3 ) 

CALL  LINE ( “I ) 

CALL  SYMBL(4) 

CALL  SIZES! .8) 

CALL  CPLOT(EE.SS) 

CALL  ANMODE 

IF(ISET2.EO.O.OR.IANS.E0.1 ) GO  TO  40 
COMPUTATIONS  AND  DISPLAY  OF  CURVES  2 ANO  3 FROM  ALLOY  LIST 
IP*ISET2 

READ! 20' I P.204) , WW , XX , EMOD , XN , SY 

CALL  CA LC( NC , EMOD .XN.SY.S.E.EP.SS.EE, EMA X , NP ) 

CALL  NPrS(NP) 

IF( IP.EQ. ISET3)  GO  TO  38 
W2  = WW 
X2  = XX 

CALL  LINE (12) 

GO  TO  39 
W3  = WW 
X3xXX 

CALL  L I NE ( 52 1 2 ) 

CALL  SYMBL(O) 
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247  CALL  CPLOT(E.S) 

248  CALL  NPTS(3) 

249  CALL  LINE(-1 ) 

250  IF(IP.EQ.ISET2)  CALL  SYMBL(1> 

251  IF'tP  nf  '5FT71  CALL  SYMBL ( 3 ) 

252  CALL  SIZES! .5) 

253  CALL  CPLOT(EE.SS) 

254  CALL  ANMODE 

255  IF( I ANS. EQ.2.0R. IP.EQ. ISET3)  GO  TO  40 

256  IP=ISET3 

257  GO  TO  37 

258  40  CALL  MOVABS ( 250 , 60 ) 

259  CALL  ANMODE 

260  WRITE(6,124) ,W1 ,X1 ,ISET1 

261  CALL  MOVABS (250,45) 

262  CALL  ANMODE 

263  I F ( IANS.EQ. 1 .OR. ISET2.EO.O)  GO  TO  41 

264  IF< IANS.EQ. 2)  WR I TE ( 6 , 1 25 ) , W2 , X2 , I SET2 

265  I F ( IANS.EQ. 3)  WR I T E ( 6 , 1 26 ) , W2 , X2 , ISET2 . W3 , X3 , ISET3 

266  CALL  MOVABS ( 250 , 30 ) 

267  CALL  ANMOOE 

268  41  WRITE(6,127) 

269  RE AO ( 5 , 20 1 ) 

270  C 

271  c 

272  C 

273  C « OPTIONS:  EXAMINE  ANOTHER  LIST,  PLOT  ANOTHER  CURVE,  OR  TERMINATE 

274  C 

275  CALL  NEWPAG 

276  CALL  I OWA I T ( 20 ) 

277  WR I TE ( 6 , 1 1 8 ) 

278  RE AD (5, 201). ANS 

279  I F (ANS . EQ. ' NEW ' ) GO  TO  42 

280  I F ( ANS . EQ. ' CUR ' ) GO  TO  34 

281  GO  TO  51 

282  42  WRITE(6,128) 

283  RE AD (5, 201 ) 

284  CALL  CHRSIZ ( 4 ) 

285  GO  TO  1 

286  C 

288  C 

239  C * ENTRY  ERROR  DIAGNOSTIC  AND  CORRECTION  OPTION 

290  C 

291  45  WRITE(6, 107) 

292  WR I TE ( 6 , 300 ) 

293  READ (5,201), AnS 

294  IF(ANS.EQ. 1 YES' ) GO  TO  (4, 34, 35), MX 

295  C 

297  C 

298  C * TERMINATION 

299  C 

300  50  CALL  CHRSIZ(3) 

301  51  CALL  NEWPAG 

302  CALL  1 OWA I T ( 20 ) 

303  CALL  MOVABS (0,700) 

304  CALL  ANMOOE 

305  WRITE(6, 129) 

306  CALL  F IN  I TT ( 0 , 650 ) 

307  STOP 
303  C 

309  

310  C 

311  100  F0RMATIT2. IX, 'PREPARE  TERMINAL  FOR  DATA  RETRIEVAL: 1 ,//,T2, 

312  1 ’ 1 . ENTER:  ®>»TTY  W.132  THEN  KEY  RETURN ',/,  T2 , 

313  2 1 2 . WAIT  FOR  RESPONSE  FROM  COMPUTER:  -••COMPLETE' ,/,T2, 

314  3 ' 3 . KEY  RETURN  AND  WAIT  FOR  PROGRAM  TO  BEGIN') 
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i — r 


u 


► 


! 

■ 


315 

316 

317 
31  9 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

332 

333 

334 

335 

336 

337 

338 
.139 

340 

341 

342 

343 
' 344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 

375 

376 

377 

378 

379 

380 

381 

382 


101  F0RMAT(T2, 'TABLE  BELOW  LISTS  RETRIEVAL  OPTIONS  AND  ENTRIES',//. 

1 T2,  ' RETRIEVAL  OPTION  1 , 14X,  ' ENTRY ' , / , 

2T2 , 'ALL  DATA  FOR  SELECTED  MAT L ' , 4X , ' ALL ' . / , 

3T2 , 'DATA  FOR  SPECIFIC  UTS  RANGE ', 3X .' ULT I MATE ',/ , 

4T2 . 1 DAT  A FOR  SPECIFIC  SY  R ANGE ' , 4X , ' Y I E LD ' , / , 

5T2 , 1 DAT  A FOR  SPECIFIC  RCH  RANGE  1 , 3X ,' HARDNESS ',/ , 

6T2 . ' MAKE  YOUR  ENTRY  NOW',/) 

102  F0RVAT(T2, 'YOUR  ENTRY  DID  NOT  CORRESPOND  TO  ANY  OF  THE  ABOVE',/, 

1 T 2 , ' IF  YOU  WISH  TO  CORRECT,  ENTER  YES',/, 

2T2, 'OTHERWISE  ENTER  END  AND  RUN  WILL  TERMINATE') 

103  FORMAT!/. T2, 'DATA  WILL  BE  RETRIEVED  FOR  A RANGE  OF  VALUES  OF  ■ , 

1 ' THE  CHOSEN  PARAMETER' ,/.T2, 'NEXT  ENTRY  DEFINES  LOWER  AND  UPPER  ' 

2 ' BOUNDS  OF  DESIRED  RANGE', 

3/, T2. 'ENTRY  UNITS:  KSI  FOR  STRENGTH  AND  ', 

4'rockwell-c  for  hardness' ,/,t2, 'separate  the  two  values  by  a •, 

5 ' COMMA . EXAMPLE:  1 50 . , 1 55 . ' , / . T2 . ' ENTER  BOUND  VALUES  NOW',/) 

105  FORMAT ( T2, ' THERE  ARE  NO  DATA  IN  THE  DATA  BANK  FOR  THE'/,T2,, 

1 'SPECIFIED  RANGE  OF  THE  SELECTED  PAR AMET ER ' , / , T2 , 

3 ' TO  CHANGE  THE  RANGE  LIMITS,  ENTER:  RANGE',/, 

4T2 , ' TO  LIST  all  THE  DATA  FOR  THE  MATERIAL,  ENTEP:  ALL',/, 

5T2 , ' TO  CHOOSE  ANOTHER  OPTION,  ENTER:  OPTION',/, 

6T2 , ‘TO  CHOOSE  ANOTHER  MATERIAL.  ENTER  NEW',/, 

7T2. 'OTHERWISE  ENTER  END  AND  RUN  WILL  TERMINATE') 

106  F0RMATIT2. ' PRINT  NOW,  CHANGE  RETRIEVAL  OPTION,  OR  TERMINATE ‘ ,/, T2 , 
1 ' 1 . TO  PRINT  A LONG-FORM  DATA-LIST  ENTER:  L0NG',/,T2, 

2'2.  TO  PRINT  A SHORT-FORM  DATA-LIST  ENTER:  SHORT', /,T2, 

3 ' 3 • TO  CHANGE  RETRIEVAL  OPTION  ENTER:  CHANGE ',/, T2 , 

4*4.  TO  TERMINATE  THE  PROGRAM  ENTER:  END',/,T2, 

5 ' MAKE  YOUR  ENTRY  NOW') 

107  FORMAT!/, T2, '»*  ERROR:  ENTRY  NOT  RECOGNIZED  **') 

118  FORMAT(T2,2X, 'OPTIONS:  MAKE  ONE  OF  THREE  ENTRIES' , / , T2 , 

1 ' FOR  A data  LIST  OF  ANOTHER  MATERIAL  ENTER  NEW',/,T2, 

2 ' FOR  A GRAPH  OF  THE  STRESS  STRAIN  CUPVE  OF  A DATA  SET', 

31X, 'ENTER  CURVE' ,/,T2, 'TO  TERMINATE  THE  PROGRAM  ENTER  END') 

119  FORMAT (T2,1X,//,T2, 'PROCEDURE  FOR  GRAPHIC  STRESS/STRAIN  CURVE', 
11X, 'DISPLAY'  ,/,T2,  '1  . ENTER  @@TTY  W,72  AND  KEY  RETURN ',/, T2 , 

2 ' 2 . WAIT  FOR  RESPONSE  FROM  COMPUTER:  -INCOMPLETE' ,/,T2, 

3 ‘ 3 . KEY  RETURN  AND  WAIT  FOR  NEXT  MESSAGE') 

120  FORMAT!/, T2, 'UP  TO  3 CURVES  MAY  BE  PLOTTED  ON  ONE  GRAPH', /,T2, 

1 'ENTER  THE  NUMBER  YOU  WISH  TO  PLOT  (ENTER  1,  2,  OR  3)') 

121  FORMAT! T2, ' YOU  ASKED  FOR  MORE  THAN  3 CURVES ',/, T2 , 

1 ' TO  MAKE  A CORRECTION  ENTER  YES',/,T2, 

2'0THERWISE  KEY  RETURN  AND  PROGRAM  WILL  TERMINATE') 

122  FORMAT!/, T2, ' 1 . ENTER  DATA  SET  NO(S).  TO  BE  GRAPHED  A KEY  RETURN', 
1 / , T2 , ' USE  COMMAS  BETWEEN  SET  NOS.  IN  ENTRY.  EXAMPLE:  236.167'. 

2/ , T2 , ' 2 . WAIT  FOR  PLOT  TO  APPEAR ' ,/ , T2 , ' 3 • TO  CONTINUE  THE', 

31X, 'PROGRAM  AFTER  PLOT,  KEY  RETURN ',/, T2 , 

4 ' 4 . THE  SCREEN  WILL  ERASE  AND  OPTIONS  TO  CONTINUE  WILL  APPEAR') 

123  FORMAT ( T2 , ' * SET  NUMBER  ERROR  « ’,/,T2, 

1 ' SET  NUMBERS  HIGHER  THAN  ',14,'  NOT  AvA I LABLE ' , // , T2 , 

2' IF  YOU  WISH  TO  MAKE  A CORRECTION  ENTER  YES',/,T2, 

3 ' OTHERW I SE  KEY  RETURN  AND  PROGRAM  WILL  TERMINATE') 

124  FORMAT! IX, 'MATL  '. A2 A2 , 2X ,' SOLI D CURVE : SET',I4,2X, 

1 'STRESS  IN  KSI  * STRAIN  IN  IN/IN') 

125  FORMAT! IX, 'DASHED:  MATL  ' , A2 , ' - ' , A2 , ' SET:',I4) 

126  FORMAT! IX, 'DASHED:  MATL  ' , A2 , ' - ' , A2 , ' SET:', 14, 

1 2X , ' BROKEN : MATL  ' , A2 , ' - ' , A2 , ' SET:', 14) 

127  F0RMAT(1X. 'PROP.  LIMIT,  0.1%-YIELD,  AND  . 2X-YI ELD  IDENTIFIED', 

11X, 'ON  CURVES' ) 

128  F0RMAT(T2, 'PROCEDURE  FOR  NEW  L I ST  I NG ' , / , T2 , 

1 ' 1 . ENTER  PPTTY  W.132  AND  KEY  RETURN ',/, T2 , 

2 ' 2 . WAIT  FOR  RESPONSE  FROM  COMPUTER:  -**COMPL£TE ' , / , T2 , 

3 ' 3 . KEY  RETURN'  ) 

129  FORMAT (T2, 'RUN  IS  OVER.  KEY  RETURN ',/, T2 , 

1 ' WHEN  PROMPT  (>)  APPEARS  ENTER  *FIN') 

C 

C 
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383 

200 

384 

201 

385 

202 

336 

203 

387 

204 

388 

C 

389 

c 

390 

C 

391 

300 

392 

393 

C 

394 

1 



2 

3 

C-- — 

4 

c 

5 

C • ! 

6 

C * 1 

7 

C • 1 

8 

C * ' 

9 

c 

10 

c 

1 1 

c 

12 

13 

14 

15 

16 

17 

c 

IS 

c 

19 

c 

20 

c * 1 

21 

c 

22 

23 

24 

c 

25 

c 

26 

c 

27 

C * 1 

28 

c 

29 

2 

30 

31 

32 

33 

34 

35 

3 

36 

37 

38 

39 

40 

41 

c 

42 

c 

43 

c 

44 

C * 1 

45 

C • 1 

46 

c 

47 

4 

48 

49 

50 

51 

52 

FORMAT ( ) 

FORMAT (A3) 

FORMAT ( 1 6X .11 , 13X, F6. 1 ,8X,F4. 1 , 13X, F7.2.4X.I4) 
FORMAT ( A2 , 1 914 ) 

FORMAT (6* .2A2.38X, F6.2.2F7.2) 


FORMAT ( T2 , 1 TO  MAKE  A CORRECTION  ENTER  YES',/,T2, 

1 'OTHERWISE  KEY  RETURN  AND  PROGRAM  WILL  TERMINATE') 

END 


SUBROUTINE  T YCODE ( JEND , N I , W , X ) 


SUBROUTINE  01  SPLAYS  TABLE  OF  CONTENTS  OF  DATA  FILE  AND  CODE  DESIGNATORS 
USED  WITH  TELETYPE  OR  PORTABLE  TERMINALS 

USER  CHOOSES  MATERIAL  AND  ALLOY  FOR  DATA  RETRIEVAL  AND  OISPLAY 
THEN  CHOICES  ARE  CHECKED  TO  ASSURE  THAT  THE  DATA  ARE  IN  THE  FILE 


DIMENSION  WA ( 7 ) , XA ( 1 9 ) t LOC ( 1 9 ) 

DATA  WA/ 'ST ' , ' TI ' , 'UR' . 'CU' , ‘WF> , ' AL‘ , 'MG'/ 

DATA  XA/'ZZ'  , 1 1 A ' , 1 2A ' , '3A' , '4A'  , '5A' , '6A'  , '7A' , '8A‘ , '9A' , 
1 ‘ IB' , ' 2 B ‘ , '3B'  , '4B' , '5B' , '60' , '7B' . 'BB' , '9B'/ 

DEFINE  FILE  10(50, 80, E.Jd) 


DISPLAY  TABLE  OF  CONTENTS  OF  MATERIALS  WITH  CODE  DESIGNATORS 

WRI TE ( 6 , 300 ) 

JEND«0 


READ  REQUESTED  MATERIAL  CODE  AND  CHECK  THAT  ENTRY  IS  CORRECT 

RE AD ( 5 . 400 ) , W 
DO  3 d»1  ,7 

IF(W.NE.WA( d)  ) GO  TO  3 

dEND= 1 

NJ*  J 

GO  TO  4 

CONTINUE 

WR I TE ( 6 , 500 ) , W 

READ( 5 , 401 ) . ANS 

IF(ANS.NE. 'YES' ) GO  TO  1 1 

WRITE (6, 501 ) 

GO  TO  2 


DISPLAY  TABLE  OF  ALLOYS  AVAILABLE  FOR  REQUESTED  MATERIAL 
READ  REQESTED  ALLOY  CODE 

IF(Nd.EQ.I)  WR I TE ( 6 , 301 ) 

I F ( Nd • EQ . 2 ) WR I TE ( 6 , 302 ) 

1 F ( Nd • EQ . 3 ) WR I TE ( 6 , 303 ) 

I F ( Nd . EQ . 4 ) WR 1 TE ( 6 , 304 ) 

I F ( Nd • EQ . 5 ) WRI TE ( 6 , 305 ) 

I F (Nd . EQ . 6 ) WR I TE ( 6, 306 ) 
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53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 
67 
66 

69 

70 

71 

72 

73 

74 

75 

76 

77 
70 

79 

80 
81 
02 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
1 1 1 
112 

113 

114 

115 
1 16 

117 

118 
1 19 
120 


5 


C 

C — 

c 

c • 
c 


I F (N J . EQ . 7 ) WRITE(6,307  ) 
WR ITE ( 6 , 308 ) 

READ ( 5 . 400 ) < X 


CHECK  THAT  ENTRY  FOR  ALLOY  IS  CORRECT  AND  THAT  DATA  ARE  AVAILABLE 

DO  6 L* 1 .19 

I F ( X . NE . XA( L) ) GO  TO  6 

JEND=1 


NL*L 
GO  TO  7 

6 CONTINUE 

WR I T E ( 6 , 502 ) ,X 
RE AD (5, 401 ) , ANS 
IF(ANS.EO. ’YES' ) GO  TO  5 
GO  TO  1 1 

7 READ ( 1 0 ‘ NJ , 402 ) , ( LOC ( L ) , L- 1 , 1 9 ) 

NI  = LOC ( NL ) 

IF(NI.NE.O)  GO  TO  10 
C 

C 

C 

C * OPTIONS:  MAKE  CORRECTIONS  OR  TERMINATE  THE  RUN 
C 

WR I TE ( 6 , 503 ) ,X 

READ( 5,401 ) . ANS 

I F ( ANS • EO . 'ALL' ) GO  TO  8 

I F ( ANS • EQ . 'BOT' ) GO  TO  9 

JEND=0 

GO  TO  11 

8 JEND=0 
GO  TO  5 

9 WR I T E ( 6 , 50 1 ) 

GO  TO  2 

C 

C 

C * WRITE  TOTAL  NUMBER  OF  ENTRIES  IN  DATA  BANK  FOR  CHOSEN  MATERIAL/ALLOY 
C • THESE  TOTALS  STORED  IN  LOCATER  FILE  AT  N-U+30 
C 

10  N.NJ+30 

RE AD ( IO'N.402) , ( LOC( L) , L*1 , 19) 

NUM=  LOC ( NL ) 

WRI TE ( 6 , 309 ) .W.X.NUM 

11  RETURN 


C 

C * 
C 

300 


301 


WRITE  FORMAT  STATEMENTS 

F0RMAT(T2, IX, 'MATERIAL  AND  ALLOY  ENTRIES  MUST  BE  CODED T2 . 

1 ' CODE  TABLE  FOR  MATERIALS  IN  DATA  BANK  1 , // , T2 , 

22 ( 'MATERIAL  CODE' . 4X ) , / , T2 , 2 ( 1 4 ( ' - ' ) , 4X ) , 

3/, T2. 'STEEL  ST 5X ,' TUNGSTEN  WF't 

4/.T2. 'TITANIUM  TI • ,5X .' ALUMINUM  AL', 

5/.T2, 'URANIUM  UR ' , 5X , ' MAGNESI UM  MG', 

6/, T2. 'COPPER  CU', 

7//, T2, 'ENTER  MATERIAL  CODE  NOW') 

F0RMAT(/,T2, 'ALLOYS  AND  ALLOY  CODES  FOR  STEEL T2 , 

1 2 ( ' ALLOY  NAME  CODE ', 6X ),/, T2 , 2( 1 8 ('-'), 6X) ,/. T2 , 

2'ARMUR  UNSPEC' ,3X, 'ZZ' ,7X, 'RARE  EARTH ', 5X ,' 1 8 ',/, T2 ,' 4340  UNSPEC', 
34X, ' 1 A ' ,7X, ' ARMCO ' . 1 OX , ’ 28 ' , / . T2 . ' MARAGE-250 ' , 5X , ' 2A ' , 7X , ' HI  ' . 
4'NICKEL' ,6X, '38' ,/,T2. 'WC-9  AL LOT ' , 5X , ' 3A ' , 7X , ' 1 340  UNSPEC', 4X 
5'4B' ,/,T2, '4140  UNSPEC '. 4X ,' 4A ', 7X .' LOW  MANGANESE ', 2X ,' 5B ' , 

6/ , T2  , 'GILO  METAL'  ,5X, ' 5A ' ,7X,  'HF-1  ALLOY ', 5X ,' 6B ' , 

7/.T2. '4340  ESR ' , 7X , ' 6A ' , 
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121 

122 

123 

124 

125 

126 

127 

128 

129 

130 
1 31 

132 

133 
1 34 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 


8/.T2, ’3%  NICKEL' ,6X, ' 7A 1 11 , 

9/.T2.'S%  NICKEL' ,6X, '8A' ,/,T2. 'AF  1 41 0 ' , 8X , • 9A ' / ) 

302  FORMAT!/, T2, 'ALLOYS  AND  ALLOY  COOES  FOR  T I TANIUM ' , // , T2 , 

1 ' ALLOY  NAME  CODE ', /T2 , 1 8 T2 , 

2 ' UNI  PEC ! F I rO ' *X, '22' ,/,T2, 

3 ' 6-6-2  ALLOY' ,4X, ' 1A' ,/,T2, 'BASE  SECTN * , 5X . • 2A ' . / ) 

303  FORMAT </,T2, 'ALLOYS  AND  ALLOY  CODES  FOR  URANIUM ',//, T2 , 

1 ' ALLOY  NAME  CODE' ,/T2,18( ) ,/,T2, 

2' DEPLETED' ,7X, ' 1 A ' ,/,T2, ' .75-TI  ALLOY ', 3X .' 2A './. T2 , 

3 ' UNALLOYED ' ,6X, ' 3A ' , / ) 

304  FORMAT(/, T2. ‘ALLOYS  AND  ALLOY  COOES  FOR  COPPER T2 , 

1 ' ALLOY  NAME  CODE' ,/T2,18( '-' ) ,/,T2, 

2'UNSPECIFIED' ,4X, 1 2Z ' ,//.T2, 'NO  OTHER  COPPER  DATA  AVAILABLE',/) 

305  FORMAT!/, T2. 'ALLOYS  AND  ALLOY  CODES  FOR  TUNGSTEN ',//, T2 , 

1 ' ALLOY  NAME  CODE' ,/T2,18( '“' ),/,T2, 

2 'UNSPECIFIED' , 4X . ' 22 ' , / . T2 , ' PM  UNSPEC' ,6X,'1A' ,/. T2 ,' SINTERED* ,7X, 
3'2A' ,/.T2, 'W-2  ALLOY' ,6X. ' 3A ' ./) 

306  FORMAT!/, T2, 'ALLOYS  AND  ALLOY  CODES  FOR  ALUMINUM ',//, T2 , 

1 1  ALLOY  NAME  CODE' ,/T2, 18( '-' ) ,/,T2, 

2 ' 7075-T6 ' ,8X, • 1 A ' ,/,T2,  1 WHISKER  COMP'  ,3X, ' 2A ' ,/) 

307  FORMAT!/, T2. 'ALLOYS  AND  ALLOY  COOES  FOR  MAGNESIUM' ,//, T2 , 

1 ' ALLOY  NAME  CODE' ,/T2,18( '-' ) ,/,T2, 

2 ' AZ~6 1 A ' ,9X, ‘ 1A' ,//,T2, 'ONLY  DATA  FOR  MAGNESIUM  AVAILABLE',/) 

308  F0RMAT(T2,  'ENTER  ALLOY  CODE  NOW) 

309  FORMAT!/, T2. 'TOTAL  NUMBER  OF  DATA  SETS  FOR  MATERIAL 
1A2,‘-',A2,'  IS', 14,/) 

C 

C * READ  FORMAT  STATEMENTS 

c 

400  FORMAT ( A2 ) 

401  FORMAT (A3) 

402  FORMAT(2X, 1914) 

C 

C * DIAGNOSTIC  FORMAT  STATEMENTS  FOR  ERRONEOUS  ENTRIES 
C 

500  FORMAT!/, T2, 'THE  MATERIAL  CODE  ENTERED  ( ' , A2 , ' ) DOES  NOT  APPEAR  ', 
1 ' IN  THE  CODE  T ABLE ' , / , T2 , ' TO  CHANGE  YOUR  ENTRY  TYPE  YES',/,T2, 

2  ‘ OTHE  Ri*  I SE  TYPE  END  AND  PROGRAM  WILL  TERMINATE',/) 

501  FORMAT(T2, 'ENTER  CORRECTED  MATERIAL  CODE') 

502  FORMAT!/, T2. 'THE  ALLOY  CODE  ENTERED  ( ' . A2 , ‘ ) DOES  NOT  APPEAR 

1 ' IN  THE  CODE  TABLE' ,/,T2, 'TO  CHANGE  YOUR  ENTRY  TYPE  YES',/,T2, 

2 ' OTHERWISE  TYPE  END  AND  PROGRAM  WILL  TERMINATE',/) 

503  F0RMAT(T2, 'DATA  FOR  ALLOY  CODE  (',A2,')  ARE  NOT  YET  IN  THE  DATA', 
1'  BANK' ,/,T2, 'TO  CHANGE  ALLOY  CODE  ENTRY  ONLY  TYPE  ALLOY', 
2/,T2,'T0  CHANGE  BOTH  MATERIAL  AND  ALLOY  ENTRIES  TYPE  BOTH',/,T2. 

3 ' OTHERWISE  type  end  AND  PROGRAM  WILL  TERMINATE',/) 

C 

END 


1 C 

2 SUBROUTINE  CODE ( LOC . JEND , NI . W , X ) 

3 C 

4 C 

5 c * subroutine  displays  table  of  contents  of  data  file  and  code  designators 

6 C * USED  WITH  CR-TUBE  GRAPHICS  TERMINALS 

7 C ♦ USER  CHOOSES  MATERIAL  AND  ALLOY  FOR  DATA  RETRIEVAL  AND  DISPLAY 

8 C * THEN  CHOICES  ARE  CHECKED  TO  ASSURE  THAT  THE  DATA  ARE  IN  THE  FILE 

9 C 

1 0  

11  c 

12  DIMENSION  WA( 7 ) , XA( 1 9) , LOC( 1 ) 

13  DATA  WA/'ST'.'TI'.'UR'.'CU'.'WF'.'AL'.'MG'/ 

14  DATA  X A/ 1 2Z ' , ' 1 A 1 , 1 2A 1 , 1 3A ' , 1 4A 1 , 1 BA 1 , ' BA ' , 1 7A 1 , 1 BA ' , 1 BA ' , 

15  1 ' IB' , '2B' , '3B' , '4B' , '5B' , '6B' , '7B' , '8B' , *9B'/ 

16  DEFINE  FILE  10(50, 80. E.Jd) 
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f r 


17 

C 

19 

c 

20 

C • 

01  SPLAY  TABLE  OF  CONTENTS  OF  MATERIALS 

WITH  CODE  DESIGNATORS 

21 

C 

22 

1 

WR I TE ( 6 , 300 ) 

23 

JEND=0 

24 

c 

26 

c 

27 

c * 

READ  REQUESTED  MATERIAL  CODE  AND  CHECK 

THAT  ENTRY  IS  CORRECT 

28 

c 

29 

2 

READ ( 5 , 400 ) ,W 

30 

DO  3 J=1  ,7 

31 

IF(W.NE.WA( J) ) GO  TO  3 

32 

<JEND=  1 

33 

NU=J 

34 

GO  TO  4 

35 

3 

CONTINUE 

36 

WR I T E ( 6 , 500 ) ,W 

37 

RE AD ( 5 , 40 1 ) , ANS 

38 

I F ( ANS . NE . 'YES' ) GO  TO  1 1 

39 

WR I T E ( 6 , 50 1 ) 

40 

GO  TO  2 

41 

42 

C 

43 

C 

44 

c * 

DISPLAY  TABLE  OF  ALLOYS  AVAILABLE  FOR 

REQUESTED  MATERIAL 

45 

c * 

READ  REQUESTED  ALLOY  CODE 

46 

c 

47 

4 

IF(NJ.EQ.I)  WR I TE ( 6 , 30 1 ) 

48 

I F ( N J . EQ . 2 ) WR ITE ( 6 , 302 ) 

49 

I F (NJ . EQ . 3 ) WR I TE ( 6 , 303 ) 

50 

I F ( N J . EQ . 4 ) WR I TE ( 6 , 304 ) 

51 

I F ( N J . EQ . 5 ) WR I TE ( 6 , 305 ) 

52 

I F ( NU • EQ . 6 ) WR ITE ( 6 , 306 ) 

53 

I F ( N J . EQ . 7 ) WR I TE ( 6 , 307 ) 

54 

5 

WR I TE ( 6 , 308 ) 

55 

READ( 5 , 400 ) , X 

56 

57 

C 

58 

c 

59 

c * 

CHECK  THAT  ENTRY  FOR  ALLOY  IS  CORRECT 

AND  THAT  DATA  ARE  AVAILABLE 

60 

c 

61 

DO  6 L= 1 ,19 

62 

I F ( X . NE . XA( L ) ) GO  TO  6 

63 

JEND= 1 

64 

NL  = L 

65 

GO  TO  7 

66 

6 

CONTINUE 

67 

WR I TE ( 6 , 502 ) , X 

68 

READ (5,401 ) ,ANS 

69 

I F ( ANS . EO . 'YES' ) GO  TO  5 

70 

GO  TO  11 

71 

7 

READ( " 0 ' NJ , 402 ) , ( LOC (L).L*1,19) 

72 

NI  = LOC ( NL ) 

73 

IF(NI.NE.O)  GO  TO  10 

74 

C 

76 

C 

77 

c * 

OPTIONS:  MAKE  CORRECTIONS  OR  TERMINATE 

THE  RUN 

78 

c 

79 

WR I TE ( 6 , 503 ) ,X 

80 

READ (5,401 ) « ANS 

81 

I F ( ANS ■ EQ . 'ALL' ) GO  TO  8 

82 

I F ( ANS . EQ. 'BOT' ) GO  TO  9 

83 

UENO=0 

84 

GO  TO  11 
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85  8 JENO-O 

86  GO  TO  5 

87  9 CALL  NEWPAG 

88  CALL  I OWAIT ( 20 ) 

89  GO  TO  1 

90  C 

91  

92  C 

93  C * WRITE  TOTAL  NUMBER  OF  ENTRIES  IN  DATA  BANK  FOR  CHOSEN  MATERIAL/ALLOY 

94  C * THESE  TOTALS  STORED  IN  LOCATER  FILE  AT  N-d+30 

95  C 

96  10  N.NJ+30 

97  READ! 10' N, 402) , ( LOC(L) , L«1 ,19) 

98  NUM=LOC(NL) 

99  WR I T E ( 6 , 309 ) , W , X , NUM 

100  11  RETURN 

101  C 

102  

103  C 

104  C * WRITE  FORMAT  STATEMENTS 

105  C 

106  300  F0RMAT(T2, IX. ’MATERIAL  AND  ALLOY  ENTRIES  MUST  BE  CODED T2 , 

107  1 ' CODE  TABLE  FOR  MATERIALS  IN  DATA  BANK  1 , // , T2 , 

108  22 ( ' MATERIAL  CODE ' , 4X ) , / , T2 . 2 ( 1 4 ( 1 - ' ) , 4X ) , 

109  3/, T2, 'STEEL  ST 5X TUNGSTEN  WF', 

110  4/, T2, ’TITANIUM  T I • , 5X , ' ALUM I NUM  AL' , 

111  5/.T2, 'URANIUM  UR ', 5X MAGNESIUM  MG', 

112  6/, T2. 'COPPER  CU* , 

113  7//, T2. 'ENTER  MATERIAL  CODE  NOW') 

114  301  FORMAT ( / , T2 . 'ALLOYS  AND  ALLOY  CODES  FOR  STEEL ',//, T2 , 

115  12( 'ALLOY  NAME  CODE ' , 6X ) , / , T2 , 2 ( 1 8 ( ' - ' ) , 6X ) , / , T2 . 

116  2 1 ARMOR  UNSPEC'  ,3X, 'ZZ* ,7X, 'RARE  EARTH ', 5X .' 1 B' ./, T2 ,' 4340  UNSPEC', 

117  34X, ' 1 A ' ,7X,  1 ARMCO ' , 1 OX , ' 28 ' , / , T2 , 'MARAGE-250' , 5X , ' 2A ' , 7X , 1 HI  • , 

118  4'NICKEL' ,6X. '38' ,/,T2, ' WC-9  ALLOY ', 5X ,' 3A ', 7X ,' 1 340  UNSPEC',4X 

119  5’4B' ,/,T2, '4140  UNSPEC 4X 4A ', 7X LOW  MANGANESE ', 2X ,' SB 1 , 

120  6/, T2, 'GILD  METAL' ,5X, '5A' ,7X, 'HF-1  ALLOY ' . 5X 6B ' , 

121  7/.T2, '4340  ESR ' , 7X , ' 6A ' , 

122  8/.T2, '3%  NICKEL'  ,6X, '7A'  , 

123  9/.T2. '5%  NICKEL' ,6X, '8A' ,/,T2. 'AF  1 4 1 0 ' , 8X , ' 9A ' / ) 

124  302  FORMAT!/, T2.  'ALLOYS  AND  ALLOY  CODES  FOR  T I T AN  I UM ' , // , T2 , 

125  1' ALLOY  NAME  CODE ', /T2 , 1 8 ('-'),/, T2 , 

126  2 'UNSPECIFIED' . 4X . ' ZZ ' , / . T2 , 

127  3 ' 6*6-2  ALLOY' ,4X, ' 1A' ,/,T2, 'BASE  SECTN ' , 5X , ' 2A ' , / ) 

128  303  FORMAT!/, T2. 'ALLOYS  AND  ALLOY  CODES  FOR  URANIUM' ,//,T2, 

129  1 ' ALLOY  NAME  CODE ', /T2 , 1 8< '-'),/, T2 , 

130  2'DEPLETED' ,7X, ' 1 A ' ,/,T2, ' .75-TI  ALLOY ', 3X ,' 2A ',/, T2 , 

131  3'UNALLOYED' ,6X, ' 3A ’ ,/) 

132  304  FORMAT!/, T2, 'ALLOYS  AND  ALLOY  CODES  FOR  COPPER ',//, T2 , 

133  1 'ALLOY  NAME  CODE ' ,/T2 . 1 8( '-'),/, T2 , 

134  2' UNSPECIFIED' ,4X, ‘ ZZ ' ,//,T2, ' NO  OTHER  COPPER  DATA  AVAILABLE',/) 

135  305  FORMAT!/, T2, 'ALLOYS  AND  ALLOY  CODES  FOR  TUNGSTEN ',//, T2 , 

136  1'ALLOYNAME  CODE', /T2.18( '-'),/. T2, 

137  2 'UNSPECIFIED' ,4X, 'ZZ' ,/,T2. ' PM  UNSPEC' ,6X, • 1A' , / , T2 , ' S INTERED ' ,7X, 

138  3 ' 2A ' , / , T2 , ' W-2  ALLOY '. 6X ,' 3A ',/ ) 

139  306  FORMAT!/. T2, 'ALLOYS  AND  ALLOY  CODES  FOR  ALUMINUM' ,//,T2, 

140  1'ALLOYNAME  CODE' ,/T2,18( '-' ) ,/,T2, 

141  2'7075-T6' ,8x, ' 1A' ,/,T2, 'WHISKER  COMP ' , 3X , ' 2A ' , / ) 

142  307  FORMAT!/, T2. 'ALLOYS  AND  ALLOY  COOES  FOR  MAGNES I UM ' , // , T2 , 

143  1'ALLOYNAME  CODE ' ,/T2 . 1 8( '-'),/, T2 , 

144  2' AZ-61A' ,9X, ' 1 A ' ,//,T2,  'ONLY  DATA  FOR  MAGNESIUM  AVAILABLE',/) 

145  308  FORMAT (T2, 'ENTER  ALLOY  CODE  NOW') 

146  309  FORMAT!/. T2 . 'TOTAL  NUMBER  OF  DATA  SETS  FOR  MATERIAL 

147  1 A2 . ' ” 1 , A2 , ' IS1, 14,/) 

148  C 

149  C • READ  FORMAT  STATEMENTS 

150  C 

151  400  FORMAT(A2) 

152  401  FORMAT(A3) 
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7 


163  402  FORMAT (2X,  1914) 

154  C 

155  C * DIAGNOSTIC  FORMAT  STATEMENTS  FOR  ERRONEOUS  ENTRIES 

156  C 

157  500  F0RMAT(/,T2, 'THE  MATERIAL  CODE  ENTERED  ( ' , A2 , 1 ) DOES  NOT  APPEAR 

150  1 ' I N THE  CODE  T ABLE ' , / , T2 , ' TO  CHANGE  YOUR  ENTRY  TYPE  YES,,/,T21 

159  2'0THERWISE  TYPE  END  AND  PROGRAM  WILL  TERMINATE1,/) 

160  501  F0RMAT(T2, ' ENTER  CORRECTED  MATERIAL  CODE') 

161  502  FORMAT ( / , T2 , ' THE  ALLOY  CODE  ENTERED  ( ' . A2 . ‘ ) DOES  NOT  APPEAR  ', 

162  1 ' IN  THE  CODE  T ABLE ' , / , T2 , ' TO  CHANGE  YOUR  ENTRY  TYPE  YES,,/,T2, 

163  2'0THERWISE  TYPE  END  AND  PROGRAM  WILL  TERMINATE',/) 

164  503  F0RMAT(T2, 'DATA  FOR  ALLOY  CODE  (',A2,')  ARE  NOT  YET  IN  THE  OATA', 

165  1*  BANK' ,/,T2, 1 TO  CHANGE  ALLOY  CODE  ENTRY  ONLY  TYPE  ALLOY', 

166  2/ , T2 , ‘TO  CHANGE  BOTH  MATERIAL  AND  ALLOY  ENTRIES  TYPE  B0TH',/,T2. 

167  3 1 OTHERWISE  TYPE  END  AND  PROGRAM  WILL  TERMINATE',/) 

168  C 

169  END 


J 


1 

2 

C 

SUBROUTINE  R ANGE ( W , X , N I , 1 1 

4 

c 

5 

c * 

SUBROUTINE  FINDS  MIN-MAX  OF 

6 

c * 

FILE  20  IS  DATABANK  FILE 

7 

c 

8 

DEFINE  FILE  20 ( 9999 , 80 , E , 

9 

IH0LD=NI 

10 

1 

NUM  = 0 

1 1 

GO  TO  (2,10,20), I DENT 

12 

c 

14 

c 

15 

c * 

ULTIMATE  STRENGTH  RANGE 

16 

c 

17 

2 

RE AD < 20 'N I, 200) .SULT.NEXT 

18 

N I =NE  XT 

19 

BIG=SULT 

20 

SMA=SULT 

21 

NUM=NUM+1 

22 

I F ( NEXT . EQ . 0 ) GO  TO  30 

23 

3 

READ l 20' NI ,200) .SULT.NEXT 

24 

NI=NEXT 

25 

IF(SULT.GT.BIG)  BIG=SULT 

26 

IF(SULT. LT.SMA)  SMA.SULT 

27 

num=num+i 

28 

I F ( NEXT . EQ . 0 ) GO  TO  30 

29 

GO  TO  3 

30 

c 

31 

c * 

YIELD  STRENGTH  RANGE 

32 

c 

33 

10 

READ( 20'NI ,201 ) , SY , NEXT 

34 

ni*ne,-;t 

35 

BIG*SY 

36 

SMA.SY 

37 

NUM.NUM+1 

38 

IF(NEXT.EQ.O)  GO  TO  30 

39 

1 1 

READ*  20'NI ,201 ) , SY , NEXT 

40 

ni*next 

41 

IF(SY.GT.BIG)  BIG=SY 

42 

I F(SY. LT.SMA)  SMA*SY 

43 

NUM=NUM+1 

44 

IF(NEXT.EQ.O)  GO  TO  30 

45 

GO  TO  11 

46 

c 

47 

c • 

ROCKWELL-C  HARDNESS  RANGE  * 

43 


48 

C 

49 

20 

RE AD ( 20 ‘NI ,202) .RCH.NEXT 

50 

NI*NEXT 

51 

IF(RCH. LE.O. AND. NEXT. GT.O)  GO  TO  20 

52 

IF(°CH.  _E.<1.  AND. NEXT. EQ.O)  GO  TO  22 

53 

BIG.RCH 

54 

SMA.RCH 

55 

NUM»NUM*1 

56 

I F (NEXT . EQ. 0 ) GO  TO  30 

57 

21 

READ( 20 ' NI , 202 ) .RCH.NEXT 

58 

NI .NEXT 

59 

IF(RCH. LE.O. AND. NEXT. GT.O)  GO  TO  21 

60 

IF(RCH. LE.O. AND. NEXT. EO.O)  GO  TO  30 

61 

IF(RCH.GT.BIG)  BIG.RCH 

62 

IF(RCH. LT.SMA)  SMA.RCH 

63 

NUM.NUM+1 

64 

IF(NEXT.EQ.O)  GO  TO  30 

65 

GO  TO  21 

66 

67 

68 

C 

C 

69 

C * 

CHANGE  OF  OPTIONS  IF  NO  HARDNESS  DATA  AVAILABLE 

70 

c 

71 

22 

WRITE'S, 100) ,W,X 

72 

IOENT  = 0 

73 

RE AD ( 5 , 203 ) ,ANS 

74 

I F ( ANS . EO . 'ALL' ) GO  TO  40 

75 

I F ( ANS . EQ . ' ULT ' ) I DENT  « 1 

76 

IFUNS.EQ.  'VIE' ) IDENT-2 

77 

I F ( I DENT. GT.O)  GO  TO  23 

78 

UEND=0 

79 

GO  TO  40 

80 

23 

NI»IHOLD 

81 

GO  TO  1 

82 

C 

c— 

84 

C 

85 

c • 

WRITE  RANGE  VALUES 

66 

c 

87 

30 

I F ( I DENT . EO • 1 ) WRITE (6, 101 ) , W , X , SMA , BIG, NUM 

88 

I F ( I DENT. EO- 2)  WR I TE ( 6 . 1 02 ) , W , X , SMA , BIG , NUM 

69 

I F ( I DENT. EQ. 3)  WRITE (6, 103) . W , X , SMA , BIG , NUM 

90 

40 

RETURN 

91 

C 

93 

C 

94 

100 

FORMAT ( T2 , ' »*  THERE  ARE  NO  HARDNESS  DATA  AVAILABLE  FOR  MATERIAL 

95 

1 A2 , ' - ' , A2  , ' **'  ,/,T2, 'ENTER  ANOTHER  OPTION  FROM 

THE  LIST  ABOVE'. 

96 

2/,T2,'0R  ENTER  END  TO  TERMINATE  RUN ',/, T2 ,' MAKE 

YOUR  ENTRY  NOW) 

97 

101 

FORMAT ( T2 , 'RANGE  OF  ULTIMATE  STRENGTH  FOR  MATL  ' 

,A2,'-',A2, 

98 

1'  IS  ' , F6. 1 , ' TO  ' , F6. 1 , ' KSI ' , / , T2 , 

99 

2 ' THE  TOTAL  NUMBER  OF  ENTRIES  WITH  NON-ZERO  DATA 

IS' ,14) 

100 

102 

FORMAT (T2, 'RANGE  OF  YIELD  STRENGTH  FOR  MATL  ',A2 

,A2, 

101 

1*  IS  ' , F7 . 2 , ' TO  ' , F7 . 2 , ' KSI',/,T2, 

102 

2 ' THE  TOTAL  NUMBER  OF  ENTRIES  WITH  NON-ZERO  DATA 

IS' ,14) 

,A2, '-' ,A3. 

103 

103 

F0RMAT(T2, 'RANGE  OF  ROCKWELL  HARONESS  FOR  MATL  ' 

104 

1 1 IS  ' ,F4.1 , ' TO  ' ,F4.1 ,/,T2, 

105 

2' THE  TOTAL  NUMBER  OF  ENTRIES  WITH  NON-ZERO  DATA 

IS' ,14) 

106 

200 

F0RMAT(30X,F6.1 ,36X,I4) 

107 

201 

F0RMAT(61X,F7.2,4X,I4) 

108 

202 

FORMAT ( 44X, F4. 1 ,24X, 14) 

109 

203 

FORMAT ( A3 ) 

110 

C 

111 

ENO 
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2 SUBROUTINE  COUNT ( I DENT , N I , XMIN , XMAX ) 

3 

4 C 

5 C * COUNTS  NUMBER  OF  DATABANK  ENTRIES  IN  SELECTED  MIN-MAX  RANGE 

6 C * FILE  20  IS  DATABANK  FILE 

7 C 

B DEFINE  FILE  20(9999,80, E, II) 

9 NUM=0 

IQ  I HOLD=N I 

11  GO  TO  (1 ,2,3), IDENT 

12  C 

14  C 

15  C • ULTIMATE  STRENGTH 

16  C 

17  1 READ(20'NI, 200), SULT, NEXT 

18  IF(SULT.GE.XMIN. AND. SULT. LE. XMAX)  NUM.NUM+1 

19  I F ( NEXT . EQ . 0 ) GO  TO  4 

20  N I *NEXT 

21  GO  TO  1 

22  C 

23  C * 0.2%-OFFSET  YIELD  STRESS 

24  C 

25  2 READ(20'NI,20l ),5Y,NEXT 

26  IF(SY. GE. XMIN. AND. SY.LE. XMAX)  NUM.NUM+1 

27  IF(NEXT.EQ.O)  GO  TO  4 

28  NI=NEXT 

29  GO  TO  2 

30  C 

31  C * rockwell-c  hardness 

32  C 

33  3 READ(20’NI, 202), RCH, NEXT 

34  IF(RCH. GE. XMIN. AND. RCH. LE. XMAX)  NUM-NUM+t 

35  IF(NEXT.EQ.O)  GO  TO  4 

36  NI.NEXT 

37  GO  TO  3 

38  C 

39  C 

40  C 

41  C « WRITE  NUMBER  OF  ENTRIES 

42  C 

43  4 WRITE(6, 100) ,NUM 

44  N I a 1 HOLD 

45  RETURN 

46  C 

48  C 

49  100  FORMAT ( T2 , ' THE  NUMBER  OF  ENTRIES  WITHIN  THE  CHOSEN  RANGE  IS', 14) 

50  200  FORMAT(30X. F6. 1 ,36X, 14) 

51  201  FORMAT (61X,F7.2,4X,I4) 

52  202  FORMAT ( 44X ,F4.1,24X,I4) 

53  C 

54  END 


1 C 

2 SUBROUTINE  LONG( W , X , NI , I DENT , XMI N , XMAX ) 

3 C 

4 C 

5 C * THIS  SUBROUTINE  LISTS  LONG-FORM  DATA  ON  CRT  OF  GRAPHICS  TERMINAL 

6 C * THERE  ARE  FOUR  OPTIONS  FLAGGED  BY  VALUE  OF  IDENT 

7 C • I DENT  a 0 : ALL  DATA  FOR  SELECTED  MATL/ALLOY 

8 C ♦ IDENT.l:  only  DATA  FOR  A USER-SELECTED  RANGE  OF  ULTIMATE  STRENGTH 

9 C • I0ENT.2:  ONLY  DATA  FOR  A USER-SELECTED  RANGE  OF  YIELD  STRENGTH 

10  C * I DENT * 3 : ONLY  DATA  FOR  A USER-SELECTED  RANGE  OF  ROCKWELL-C  HARDNESS 
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11  c 

1 2 

13  C 

14  C * FILE  20  IS  DATABANK 

15  C 

16  DEFINE  FILE  20 ( 9999 , 80 , E , 1 1 ) 

17  C 

19  C 

20  1 CALL  NEWPAG 

21  CALL  I OWA I T ( 20 ) 

22  WRITE(6,100),W,X 

23  IF(IDENT.EO.O)  GO  TO  2 

24  I F ( I DENT . EG • 1 ) WR I TE ( 6 , 1 01 ) , XMI N . XMAX 

25  IF( I DENT. ED. 2)  WR I TE( 6 . 1 02 ) , AMIN , XMAX 

26  I F ( IDENT.EQ.3)  WRI TE ( 6 . 1 03 ) , XMI N , XMAX 

27  2 K=0 

28  K I =0 

29  WR ITE(6,104) 

30  3 READ! 20 1 NI , 200 ) , SN , BOOK , NPAGE , NTYPE , NDATE , KTEMP , KEDOT , SULT, 

31  1 E LONG, RA, RCH, EMOD.XN.SY, SEE, NEXT, I SET 

32  NI=NEXT 

33  IF(NTYPE.NE. 1 ) GO  TO  4 

34  IF(IDLNT.EQ.O)  GO  TO  5 

35  IFdDENT.EQ.1 -AND. SULT. LE. XMAX. ANO. SULT. GE.XMIN)  GO  TO  5 

36  IF(IOENT.EO. 2. AND. SY.LE. XMAX. AND. SY. GE.XMIN)  GO  TO  5 

37  IF(IDENT. EQ. 3. AND. RCH. LE. XMAX. AND. RCH. GE.XMIN)  GO  TO  5 

38  4 IF(NEXT.EQ.O)  GO  TO  9 

39  GO  TO  3 

40  5 XP= 1 . /XN 

41  SPL=SY«5.E-3**XP 

42  SY1«SY»5.E-1**XP 

43  WRITE (6, 105) , I SET , KTEMP , KEDOT , EMCD , XN , SP L , SY1 , SY , SULT , SEE . 

44  1 E LONG , RA , RCH .NDATE , BOOK , NPAGE , SN 

45  K.K+1 

46  KI«KI+1 

47  6 IF(NEXT.EO.O)  GO  TO  9 

48  IF(K.E0.56)  GO  TO  8 

49  I F ( MOO ( K I , 5 ) . EQ. 0 ) GO  TO  7 

50  GO  TO  3 

51  7 WR I TE ( 6 , 106) 

52  K*K+ 1 

53  GO  TO  3 

54  8 WRITE(6,107) 

55  READ( 5 , 201 ) . ANS 

56  GO  TO  1 

57  9 WRITE<6, 103) ,W.X 

58  RE AD( 5 , 20 1 ) , ANS 

59  RETURN 

60  C 

61  C 

62  C 

63  100  F0RMAT(T2, ' LONG-FORM  LIST  OF  OATA  FOR  MATL  CODE:  ,,A2,'-,.A2) 

64  101  F0RMAT(T2, ' LIST  INCLUDES  ONLY  DATA  FOR  TESTS  WHERE  THE 

65  1 1 ULT IMATE  STRENGTH  WAS  IN  THE  RANGE:  '.F6.1,'  TO  '.F6.1,'  KSI') 

66  102  F0RMAT(T2, ' LIST  INCLUDES  ONLY  DATA  FOR  TESTS  WHERE  THE  ', 

67  1 ‘YIELD  STRENGTH  WAS  IN  THE  RANGE:  '.F7.2,'  TO  ',F7.2,'  KSI') 

68  103  F0RMAT(T2, ' LIST  INCLUDES  ONLY  DATA  FOR  TESTS  WHERE  THE  ', 

69  1 1 ROCKWELL-C  HARDNESS  WAS  IN  THE  RANGE:  *^4.1,'  TO  ',F4.1,) 

70  104  FORMAT ( T2 , 1 * * MAKE  NOTE  OF  ANY  SET  NUMBERS  TO  BE  GRAPHED  *•',//. 

71  1 T2 , ' SET  TEMP  LOG  MODULUS ', 4X ,’ R-0 5X ,' P- L IM •, 2X , 

72  2“ . 1%-YLD' ,2X. ' .2%-YLD' . 4X , 1 UTS ' , 4X , ' SEE ' , 2X , 1 ELNG* ,3X, 'R-A'  ,2X, 

73  3 1 ROCK  TEST'  ,4X, ‘NTBK-'  ,4X, 'SPEC' ,/,T2, ' NO  DEG-C  EDOT  ' 

74  4'M-PSI 1 ,4X, ' EXPNT ' , 2X , 3 ( 3X , 'KSI 1 ,3X) , 1 X , 2 ( 1 X , 'KSI 1 ,3X), 1 PCT ' , 

75  53X , 1 PCT' ,2X, 'C-HD1 ,3X, ’DATE' ,4X, 'PAGE’ , 6X, 'NO* ,/) 

76  105  F0RMAT(T2, 14, 2X. 15, 3X, 1 2 . 3X , F6 . 2 , 2X , 4 ( F7 . 2 , 2X ) , F6 . 1 ,2X, 

77  14(F4. 1 ,2X). I6.2X.A4, , I2.2X.A6) 

78  106  FORMAT ( 1 H ) 
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79  107  FORMAT)/, T2, ‘‘-THERE  ARE  MORE  DATA  FOR  THIS  ALLOY  IN  THE  FILE1, 

80  1 1 X . ‘ • MAKE  NOTES  NOW  AND  TO  CONTINUE  KEY  RETURN1) 

81  108  F0RMAT(/,T2, '••LISTING  IS  COMPLETE  FOR  MATL  •1A2,'-',A2t 

82  11X,'  • MAKE  NOTES  NOW  AND  TO  CONTINUE  KEY  RETURN') 

83  C 

84  C 

85  C 

86  200  FORMAT ( A6 , 4X , A4 , 1 2 , 1 1 , 1 6 , 1 5 , 1 2 . F6 . 1 ,3F4. 1 , F6 . 2 , 2F7 . 2 , F4 . 1 ,2 14 ) 

87  201  FORMAT ( A3 ) 

88  C 

89  END 


1 C 

2 SUBROUTINE  LPRT ( W , X , NI , I DENT , XMI N , XMAX ) 

3 C 

4 C 

5 C • THIS  SUBROUTINE  LISTS  LONG-FORM  DATA  ON  TELETYPE  TERMINAL 

6 C • THERE  ARE  FOUR  OPTIONS  FLAGGED  BY  VALUE  OF  IDENT 

7 C * I DENT  * 0 : ALL  DATA  FOR  SELECTED  MATL/ALLOY 

8 C • I OENTs 1 : only  data  for  a user-selected  range  of  ultimate  strength 

9 C • IDENT *2 : ONLY  DATA  FOR  A USER-SELECTED  RANGE  OF  YIELD  STRENGTH 

10  C • IDENT  = 3 : ONLY  DATA  FOR  A USER-SELECTED  RANGE  OF  ROCKWELL-C  HARDNESS 

1 1 C 

12  C 

13  C 

14  C * FILE  20  IS  DATABANK 

15  C 

16  DEFINE  FILE  20 ( 9999 , 80 , E , 1 1 ) 

17  C 

18  C- 

19  C 

20 
21 
22 

23 

24 

25  2 

26 

27  3 

28 

29 

30 

31 

32 

33 

34 

35  4 

36 

37  5 

38 

39 

40 

41 

42 

43  6 

44 

45 

46  7 

47 

48 

49  C 

50  C- 

51  C 

52  100  FORMAT (//, T2, ' LONG-FORM  LIST  OF  DATA  FOR  MATL  CODE:  ',A2,'-',A2) 

53  101  F0RMAT(T2, ' LIST  INCLUDES  ONLY  DATA  FOR  TESTS  WHERE  THE 


WRITE (6, 100) ,W,X 
I F ( IDENT. EQ.O)  GO  TO  2 
I F ( IDENT. EQ. 1 ) WRITE(6, 101 ) .XMIN.XMAX 
I F ( I DENT. EQ. 2)  WRITE(6, 102) .XMIN.XMAX 
I F ( I DENT. EO. 3)  WRITE (6. 103) .XMIN.XMAX 
K I =0 

WRITE16, 104) 

READ( 20'NI .200) , SN , BOOK . NPAGE , NT YPE , NOATE , KTEMP , KEDOT , SULT , 
1 ELONG , RA , RCH . EMOD , XN , SY , SEE .NEXT , I SET 
NI =NEXT 

IF(NTYP=.NE. 1 ) GO  TO  4 
I F ( IDENT. EQ.O)  GO  TO  5 

IF(  IDENT. EO.  1 .AND. SULT. LE. XMAX. AND. SULT.GE.XMIN)  GO  TO  5 
IF(IDENT.E0.2. AND. SY.LE. XMAX. AND. SY.GE.XMIN)  GO  TO  5 
IF(IDENT.EQ.3.AND.RCH.LE.XMAX.AND.RCH.GE.XMIN)  GO  TO  5 
I F( NEXT .EQ.O)  GO  TO  7 
GO  TO  3 
XP=1 . / XN 

SPL*SY»5.E-3»*XP 

SY1*SY*5.E-1**XP 

WRITE (6, 105) . I SET, KTEMP. KEDOT, EMOD , XN , SP L , SY 1 , SY , SULT , SEE . 

1 ELONG . R A . RCH , NDAT  E , BOOK , NPAGE , SN 
K I *K I + 1 

IF(NEXT.EQ-O)  GO  TO  7 
I F I MOD ( K I ,5) .EQ.O)  WRITE(6,106) 

GO  TO  3 

WRITE(6, 107) ,W,X 
RE AD ( 5 , 20 1 ) ,ANS 
RETURN 
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54 

55 

102 

56 

57 

103 

58 

59 

104 

60 

61 

62 

63 

64 

65 

105 

66 

67 

106 

68 

107 

69 

70 

C 

71 

C 

72 

c 

73 

200 

74 

201 

75 

C 

76 

1 'ULTIMATE  STRENGTH  WAS  IN  THE  RANGE:  \F6.1,'  TO  • . F6 . 1 . * KSZ ‘ ) 
F0RMAT(T2, ' LIST  INCLUDES  ONLY  DATA  FOR  TESTS  WHERE  THE 
1 'YIELD  STRENGTH  WAS  IN  THE  RANGE:  '.F7.2,'  TO  '.F7.2,'  KSI' ) 
FORMAT(T2,  1 LIST  INCLUDES  ONLY  DATA  FOR  TESTS  WHERE  THE 
1 ' R0'"KWCI  l.-<~  HARDNESS  WAS  IN  THE  RANGE:  * , F4 . 1 . * TO  ',F4.1.) 
FORMAT (/,T2,112( ),/, 

1 T2 . 1 SET  TEMP  LOG  MODULUS 4X R-0 5X , 1 P- L IM 2X , 

2' . 1%-YLD'  ,2X,  1 . 2%-YlD ' , 4X , ' UT S ' , 4X , ' SEE ' , 2X , ' ELNG ' . 3X , ' R-A ' , 2X , 
3 1 ROCK  TEST' ,4X,  'NTBK-'  ,4X,  'SPEC'  ,/,T2,  1 NO  OEG-C  EDOT  ' 
4'M-PSI ' ,4X. ' EXPNT ' , 2X , 3 ( 3X . ' KS I ' , 3X ) , 1 X , 2 ( 1 X . ' KS I ' , 3X ) , 1 PCT ' , 
53X , ' PCT' ,2X. 'C-HD* ,3X, 'DATE' ,4X, 'PAGE' ,6X, 'NO' ./) 

FORMAT ( T2. 14, 2X, I5.3X, 1 2 , 3X , F6 . 2 , 2X , 4 ( F7 . 2 , 2X ) , F6 . 1 ,2X, 

1 4 ( F4. 1 ,2X) , I6.2X, A4, '-' , I2.2X, A6) 

FORMAT ( 1 H ) 

FORMAT (/,T2, '»*LISTI NG  IS  COMPLETE  FOR  MATL  ',A2,'-',A2, 

1 1 X . ' * TO  CONTINUE  THE  PROGRAM  KEY  RETURN') 


FORMAT (A6.4X.A4, 12,11 , 16, 15, 12 , F6. 1 ,3F4. 1 , F6.2, 2F7.2, F4. 1 .214) 
FORMAT( A3) 

END 


1 C 

2 SUBROUTINE  SHORT ( W , X , N I . I DENT , XMIN , XMAX ) 

3 C 

4 C 

5 c * this  subroutine  lists  short-form  data  on  crt  of  graphics  terminal 

6 c * there  are  four  options  flagged  by  value  of  ident 

7 C * I DENT  = 0 : ALL  DATA  FOR  SELECTED  KATL/ALLOY 

8 C * IDENT. 1:  ONLY  DATA  FOR  A USER-SELECTED  RANGE  OF  ULTIMATE  STRENGTH 

9 C * IDENT  = 2 : ONLY  DATA  FOR  A USER-SELECTED  RANGE  OF  YIELD  STRENGTH 

10  C • IDENT  = 3 : ONLY  DATA  FOR  A USER-SELECTED  RANGE  OF  ROCKWELL-C  HARDNESS 

1 1 C 

12  

13  C 

14  DIMENSION  R AC ( 2 ) , XPC ( 2 ) , UTC ( 2 ) , ELC ( 2 ) , HDC ( 2 ) , EMC ( 2 ) , SYC ( 2 ) , 

15  1 I SC ( 2 ) ,KTP(2) 

16  C 

18  C 

19  C * FILE  20  IS  DATABANK 

20  C 

21  DEFINE  FILE  20(9999, 80, E, II) 

22  C 

24  C 

25  1 CALL  NEWPAG 

26  CALL  I QwA I T ( 20 ) 

27  WRITE (6, 100) ,W,X 

28  IF( IDENT. EO.O)  GO  TO  2 

29  I F ( I DENT . EO • 1 ) WR I TE ( 6 , 1 0 1 ) , XMI N , XMAX 

30  I F ( IDENT. EO. 2)  WR I TE ( 6 , 1 02 ) , XMIN , XMAX 

31  I F ( IDENT. EO. 3)  WR I TE ( 6 , 1 03 ) , XM IN , XMAX 

32  2 WR I TE ( 6 , 1 04 ) 

33  K.O 

34  KI*0 

35  3 IX«0 

36  4 IX*IX+1 

37  5 READ! 20'NI .200) ,NTYPE,KTP( IX) ,UTC( IX) ,ELC( IX) ,RAC( IX) , 

38  1HDC( IX) ,EMC( IX) ,XPC( I X ) , SYC ( I X ) , NEXT , I SC ( I X ) 

39  NI .NEXT 

40  IF(NTYPE.NE. 1 ) GO  TO  6 

41  IF(IDENT.EQ.0.AN0.IX.E0.2)  GO  TO  8 
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42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 


IF(IDENT. EQ.O. AND. NEXT. EQ.O)  GO  TO  8 
IF( IDENT.EQ.O)  GO  TO  4 

I F ( IDENT.EQ. 1 . AND.UTC( IX) . LE . XMAX . AND . UTC ( IX) .GE.XMIN)  GO  TO  7 
IF( IDENT.EQ.2.AND.SYC( IX) . LE.XMAX.AND.SYC( IX) .GE.XMIN)  GO  TO  7 
IF( IDENT.EO- 3. AND. HDC( IX) . LE. XMAX. AND. HDC( IX) .GE.XMIN)  GO  TO  7 

6 IF( IX . EQ. 1 . AND. NEXT . NE . 0)  GO  TO  3 
I F ( IX . EQ.2. AND. NEXT . NE . 0 ) GO  TO  5 
IFUX.EQ.1  ) GO  TO  11 

I X=  1 

GO  TO  8 

7 IF( IX. EQ.2)  GO  TO  8 

I F ( I X . f Q . 1 .AND. NEXT. EQ.O)  GO  TO  8 
GO  TO  4 

8 I F ( I X . EQ . 2 ) WRITE(6, 105) . ( ISC( IX) ,KTP( IX) ,EMC( IX) , 

1SYC( IX) ,UTC( IX) ,XPC( IX) ,ELC( IX) ,RAC( IX) ,HDC( IX) , IX* 1 ,2) 

IF( IX. EQ. 1 ) WRITE! 6, 106) , I SC ( IX) ,KTP( IX)  ,EMC( IX) , 

1SYC( IX) ,UTC( IX) ,XPC( IX) , ELC( IX) ,RAC( IX) ,HDC( IX) 

K = K+  1 
K I =K I + 1 

I F ( NEXT . EQ . 0 ) GO  TO  11 
I F ( K . EQ . 56 ) GO  TO  10 
I F ( MOD ( KI .5) .EQ.O)  GO  TO  9 
GO  TO  3 

9 WR I TE ( 6 . 1 07 ) 

K = K+1 

GO  TO  3 

10  WR I TE ( 6 , 1 08 ) 

RE AD ( 5 , 20 1 ) , ANS 
GO  TO  1 

11  WRITE! 6, 109) ,W,X 
RE AD ( 5 , 201 ) . ANS 
RETURN 

C 

c 

c 

100  FORMAT ( T2, 1 X LIST ING  OF  SHORT-FORM  TEST  DATA  FOR  MATL/ALLOY 
1A2, ,A2) 

101  FORMAT! T2,  1 LIST  INCLUDES  DATA  FOR  TESTS  WHERE  THE 

1 'ULTIMATE  STRENGTH  WAS  IN  THE  RANGE:  1 , F6 . 1 , ' TO  ',F6.1.'  KS I ' ) 

102  F0RMATIT2. 1 LIST  INCLUDES  DATA  FOR  TESTS  WHERE  THE 

1 'YIELD  STRENGTH  WAS  IN  THE  RANGE:  ',F7.2,'  TO  ‘,F7.2,'  KSi') 

103  F0RMAT(T2, ' LIST  INCLUDES  DATA  FOR  TESTS  WHERE  THE 

I'ROCKWELL-C  HARDNESS  WAS  IN  THE  RANGE:  '.F4.1,'  TO  ',F4.1) 

104  FORMAT ( T2 , 1 * MAKE  NOTE  OF  ANY  SET  NUMBERS  TO  BE  GRAPHED',//, 

1T2.2! IX, 'SET  TEMP  EMOD  . 2%-SY ' ,5X, 'UTS' ,5X, 'R-Q  ELNG 1 , 

2'  R-A  R-C ,3X),/,T2,2(1X, 'NO  DEG-C  M-PSI ' , 2X , 2 ( 3X , ' KSI ' , 
33X ) , IX, ’ EXPNT ' ,3X,2( 1 PCT ' ,3X) , 'hRD' ,3X) ,/) 

105  FORMAT (T2, 2( 14 , 1 X , 1 5 , 2X , F6 . 2 , 2X , F7 . 2 , 2X , F6 . 1 , 2X , F7 . 2 . 2X , 

1 3 ( F4. 1 ,2X) , IX)  ) 

106  FORMAT (T2, 14, IX, 1 5 , 2X , F6 . 2 , 2X , F7 . 2 , 2X , F6 . 1 ,2X,F7.2,2X, 

1 3 ( F4 . 1 , 2X ) ) 

107  FORMAT! 1H  ) 

108  FORMAT!/, T.1, '••THERE  ARE  MORE  DATA  FOR  THIS  ALLOY  IN  THE  FILE', 
1 1 X , ' • MAKE  NOTES  NOW  AND  TO  CONTINUE  KEY  RETURN') 

109  FORMAT!/, T2, ' **LISTING  IS  COMPLETE  FOR  MATL  ',A2,'-',A2, 

11X,'  • MAKE  NOTES  NOW  AND  TO  CONTINUE  KEY  RETURN') 

C 

C 

200  FORMAT! 16X, II , 6X , 1 5 , 2X , F6 . 1 , 3F4. 1 , F6 . 2 . 2F7 . 2 , 4X , 

1214) 

201  FORMAT ( A3 ) 

C 

ENO 
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1 C 

2 SUBROUTINE  SPRT ( W , X , NI , I DENT , XMIN , XMAX ) 

3 

4 C 

5 C * THIS  SUP'IOU’'  I NF  LISTS  SHORT-FORM  DATA  ON  TELETYPE  TERMINAL 

6 C * THERE  ARE  FOUR  OPTIONS  FLAGGED  BY  VALUE  OF  IDENT 

7 C • I DENT » 0 : ALL  DATA  FOR  SELECTED  MATL/ALLOY 

8 C * IDENT-i:  ONLY  DATA  FOR  A USER-SELECTED  RANGE  OF  ULTIMATE  STRENGTH 

9 C • I DENT » 2 : ONLY  DATA  FOR  A USER-SELECTED  RANGE  OF  YIELD  STRENGTH 

10  C * I DENT  = 3 : ONLY  DATA  FOR  A USER-SELECTED  RANGE  OF  ROCKWELL-C  HARDNESS 

1 1 C 

12  

13  C 

14  DIMENSION  RAC(2) , X PC ( 2 ) ■ UTC ( 2 ) , E LC ( 2 ) , HOC ( 2 ) , EMC (2) , SYC ( 2 ) , 

15  1ISC(2),KTP(2) 

16  C 

17  C 

18  C 

19  C • FILE  20  IS  DATABANK 

20  C 

21  DEFINE  FILE  20 ( 9999 , 80 , E , 1 1 ) 

22  C 

23  C 

24  C 

25  1 WRITE(6, 100) ,W,X 

26  IF( IDENT. EQ.O)  GO  TO  2 

27  IF( IDENT. EO- 1 ) WR I TE ( 6 , 1 01 ) , XM IN , XMAX 

28  IF( IDENT. EO-2)  WR I TE ( 6 , 1 02 ) , XMIN , XMAX 

29  I F ( IDENT. EO- 3)  WR I TE ( 6 , 1 03 ) , XMI N , XMAX 

30  2 WR I T E ( 6 , 1 04 ) 

31  K I =0 

32  3 IX=0 

33  4 IX=IX+1 

34  5 READ( LO 1 NI ,200) ,NTYPE,KTP( IX) ,UTC( IX) ,ELC( IX) , 

35  1RAC( IX) ,HDC( IX) ,EMC( IX) ,XPC(IX) ,SYC( IX) , NEXT . ISC( IX) 

36  NI «NEXT 

37  IF(NTYPE.NE. 1 ) GO  TO  6 

38  IF( IDENT. EO.O. AND. IX. EO. 2)  GO  TO  B 

39  IF( IDENT. EO.O. AND. NEXT. EO.O)  GO  TO  8 

40  I F ( IDENT. EO-O)  GO  TO  4 

41  I F ( IDENT. EO. 1 .and. UTC( IX) .LE. XMAX. AND. UTC(IX) .GE. XMIN)  GO  TO  7 

42  IF( IDENT. EO. 2. AND. SYC(IX) . LE. XMAX. AND. SYC( IX) .GE. XMIN)  GO  TO  7 

43  I F ( IDENT. EO. 3. AND. HOC ( IX) . LE. XMAX . AND . HOC ( I X ) . GE . XMIN ) GO  TO  7 

44  6 I F ( IX. EO. 1 • AND. NEXT. NE.O)  GO  TO  3 

45  IF(IX.E0.2.AND.NEXT.NE.0)  GO  TO  5 

46  I F ( IX.EO. 1 ) GO  TO  9 

47  IX*1 

48  GO  TO  8 

49  7 I F ( IX.EO. 2)  GO  TO  0 

50  IFdX.EQ.I -AND. NEXT. EO-O)  GO  TO  8 

51  GO  TO  4 

52  8 I F ( I X . EO . 2 ) WRITE(6, 105) . ( I SC ( IX) ,KTP( IX ) ,EMC( IX) , 

53  1 S Y C ( IX) ,UTC( IX) ,XPC(IX) . E LC ( I X ) , RAC (IX) ,HDC(IX) , IX-1 ,2) 

54  IF(IX.EO.I)  WRITE(6, 106) , I SC ( IX) , KTP( IX) , EMC( IX) f 

55  1SYC( IX) ,UTC( IX) ,XPC( IX) .ELC(IX) , RAC ( IX) ,HDC(IX) 

56  KI-KI+1 

57  I F ( NEXT . EO • 0 ) GO  TO  9 

58  IF(f.10D(KI  ,5)  .EO.O)  WRITE(6,107) 

59  GO  TO  3 

60  9 WRITE(6, 108) ,W.X 

61  READ( 5.201 ) .ANS 

62  RETURN 

63  C 

64  

65  C 

66  100  F0RMAT(T2, IX, 'LISTING  OF  SHORT-FORM  TEST  DATA  FOR  MATL/ALLOY 

67  1 A2 , 1 - ' , A2 ) 

60  101  F0RMAT(T2, ' LIST  INCLUDES  DATA  FOR  TESTS  WHERE  THE 
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69  1'ULTIMATE  STRENGTH  WAS  IN  THE  RANGE:  '.F6.1,'  TO  ',F6.1,’  KSI'I 

70  102  FORMAT ( T2 , 'LIST  INCLUDES  DATA  FOR  TESTS  WHERE  THE  1 , 

71  1 'YIELD  STRENGTH  WAS  IN  THE  RANGE:  '.F7.2,'  TO  '.F7.2.'  KSI') 

72  103  F0RMAT(T2, ' LIST  INCLUDES  DATA  FOR  TESTS  WHERE  THE 

73  1'ROCKWElL-C  HARDNESS  WAS  IN  THE  RANGE:  * , F4 . 1 , * TO  ’,F4.1) 

74  104  FORMAT!/. T2, 62 ( ,3X,62( ,/, 

75  1T2.2( IX, 'SET  TEMP  EMOD  . 2K-SY ' , 5X , ' UTS ' , 5X , ' R-0  ELNG', 

76  2'  R-A  R-C' ,3X),/,T2,2(1X, 'NO  DEG-C  M-PSI ' , 2X , 2 ( 3X. ' KSI ' , 

77  33X) . IX. ‘ EXPNT 1 , 3X . 2 ( ' PCT ' . 3X ) , 1 HRD ' ,3X) ,/) 

78  105  F0RMAT(T2,2(I4,1X,I5.2X,F6.2,2X,F7.2.2X,F6.1 .2X.F7.2.2X. 

79  1 3 ( F4 . 1 , 2X ) , 1 X ) ) 

80  106  F0RMAT(T2. 14. IX, I5.2X.F6.2.2X, F7.2.2X.F6. 1 .2X.F7.2.2X, 

81  13(F4.1 ,2X)) 

82  107  FORMAT ( 1 H ) 

83  108  FORMAT!/, T2, '**LISTING  IS  COMPLETE  FOR  MATL  ,,A2,'-',A2. 

84  1 1 X ' * TO  CONTINUE  THE  PROGRAM  KEY  RETURN') 

85  C 

87  C 

88  200  FORMAT (16X.I1 ,6X,I5,2X,F6.1 , 3F4 . 1 .F6.2.2F7.2.4X, 

69  1214) 

90  201  FORMAT ( A3) 

91  C 

92  END 


1  

2 SUBROUTINE  PORPRT ( W , X , N I , I DENT , XMIN , XMAX ) 

3 

4 C 

5 C « this  SUBROUTINE  LISTS  EXPANDED  SHORT-FORM  DATA  ON  PORTABLE  TERMINAL 

6 C * THERE  ARE  FOUR  OPTIONS  FLAGGED  BY  VALUE  OF  IDENT 

7 C * IDENT  = 0 : ALL  DATA  FOR  SELECTED  MATL/ALLOY 

8 C • IDENT=i:  ONLY  DATA  FOR  A USER-SELECTED  RANGE  OF  ULTIMATE  STRENGTH 

9 C * I DENT=2 : ONLY  DATA  FOR  A USER-SELECTED  RANGE  OF  YIELD  STRENGTH 

10  C • IDENT=3:  ONLY  DATA  FOR  A USER-SELECTED  RANGE  OF  ROCKWELL-C  HARDNESS 

1 1 C 

13  C 

14  C * FILE  20  IS  DATABANK 

15  C 

16  DEFINE  FILE  20 ( 9959 , 80 , E . 1 1 ) 

17  C 

19  C 

20  WR1TE(6, 100) ,W,X 

21  IF( IDENT. EQ.O)  GO  TO  2 

22  IF( IDENT. EO. 1 ) WR I TE ( 6 , 1 01 ) , XMIN , XMAX 

23  IF! IDENT. EO- 2)  WR I TE ( 6 , 1 02 ) , XM IN . XMAX 

24  IF( IDENT. EO- 3)  WR I TE ( 6 , 1 03 ) , XMI N , XMAX 

25  2 K I =0 

26  WR I TE ( 6 , 1 04  ) 

27  3 RE AD ( 20 ' N I , 200 ) , NT YPE , KTEMP , KEDOT , SULT , 

28  1 E LONG , RA , RCH , EMOD , XN , SY , NEXT , I SET 

29  NI -NEXT 

30  IF(NTYPE.NE.I)  GO  TO  4 

31  IF! IDENT. EQ.O)  GO  TO  5 

32  IF! IDENT. EO. 1 .AND. SULT. LE. XMAX. ANO. SULT. GE. XMIN)  GO  TO  5 

33  IF( IDENT. EO- 2. AND. SY.LE. XMAX. AND. SY.GE. XMIN)  GO  TO  5 

34  IF(IDENT.EQ.3.AND.RCH.LE.XMAX.AND.RCH.GE.XMIN)  GO  TO  5 

35  4 I F ( NEXT .EQ.O)  GO  TO  6 

36  GO  TO  3 

37  5 XP«1 . /XN 

38  SPL*SY»5.E-3«*XP 

39  SY1«SY«.5»*XP 

40  WRITE! 6. 105 ) , I SET , KTEMP , KEDOT , EMOD. XN.SPL.SY1 . SY , SULT , ElONG, 
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41  1RA.RCH 

42  K I *K I + 1 

43  I F (NEXT . EO ■ 0 ) GO  TO  6 

44  1 F ( MOD ( K 1 , 5 ) . EQ . 0 ) WRITE(6,106) 

45  GO  T0  3 

46  6 WRITE(6, 107) ,W,X 

47  READ( 5,201 ) .ANS 

48  RETURN 

49  C 

50  C— 

51  C 

52  100 

53 

54 

55  101 

56 

57  102 

58 

59  103 

60 

61  104 

62 

63 

64 

65  105 

66 

67  106 

68  107 

69 

70 

71  C 

72  C— 

73  C 

74  200  FORMATf 16X, II .6X.I5, 12. F6. 1 ,3F4. 1 , F6 . 2 . 2F7.2 ,4X , 214) 

75  201  FORMAT(A3) 

76  C 

77  END 


1  C 

2 SUBROUTINE  S YS I ZE ( I ANS , I SET1 . I SET2 , I SETS ) 

3 c 

4 C 

5 C • FINDS  SET  WITH  maximum  SY  and  ASSIGNS  ISET1  TO  THAT  SET 

6 C * ISET1  IS  KEY  SET  TO  ESTABLISH  GRAPH  SIZE  FOR  S-E  CURVE  PLOT 

7 C • PROCEDURE  ASSURES  THAT  ALL  POINTS  OF  ALL  CURVES  WILL  BE  ON  GRAPH 

8 C 

9 

10  C 

lie*  FILE  20  IS  DATABANK  FILE 

12  C 

13  DEFINE  FILE  20(9999, 80. E, II) 

14  C 

15  C 

16  C 

17  IF( IANS. GT. 1 .AND. ISET2.E0.0)  IANS-1 

18  IF( IANS.EQ. 1 ) RETURN 

19  I F ( IANS.EQ. 3. AND. ISET3.EQ.0)  IANS*2 

20  1 F ( IANS.E0.2)  GO  TO  2 

21  REA0(20' ISET1 ,200) ,SY1 

22  READ(20'ISET2,200),SY2 

23  READ120' ISET3.200)  ,SY3 

24  1FISY1 .GE.SY2.AND.SY1 .GE.SY3)  RETURN 

25  IF(SY2.GE.SY3)  GO  TO  1 

26  IH0L01*ISET3 

27  IH0LD2* I SET  1 


FORMAT (//,T2.75( ) ,//.T2. 1 EXPANDED  SHORT-FORM'  , 

1'  LIST  of  DATA  FOR  MATL  CODE:  ' , A2 , ' - ' , A2 . / , T2 , 

2 1 ( FOR  A FULL  LONG-FORM  LIST  USE  TELETYPE  OR  CRT  TERMINAL)',/) 
F0RMAT(T2, ' LIST  INCLUDES  ONLY  TESTS  WITH 
1 1 ULTIMATE  FROM  ',F6.1,‘  TO  '.F6.1,'  KSI',/) 

F0RMAT(T2, ' LIST  INCLUDES  ONLY  TESTS  WITH 
1 'YIELD  FROM  '.F7.2,'  TO  ’.F7.2,'  KSI',/) 

F0PMAT(T2, ' LIST  INCLUDES  ONLY  TESTS  WITH 
1 ' R-C  HARDNESS  FROM  ' , F4 . 1 , ' TO  ' ,F4.1,/) 

FORMAT ( T3 , 

1 ' SET  TEMP  LOG  EMOD ' , 5X , ' R-0 ' , 4X , ' P-L IM  .1%-YLD  ,2%-YLD'. 

23X , 'UTS  ELNG  R-A  ROCK ' , / , T4 , ' NO  DG-C  EDOT  M-PSI  EXPNT ' , 

33 (4X, 'KSI ' . IX) ,3X, 'KSI ' , 1X,2(2X, 'PCT ' ) , IX, 'C-HD' ,/) 

FORMAT (T2. 14, I5.2X, 1 2 , 2X , F6 . 2 , 4 ( IX, F7.2) , IX, F6. 1 , 

1 3( 1 X , F4 . 1 )) 

F0RMAT(1H  ) 

F0RMAT(/.T2,75( '-'),//, 

1T2, ' '‘LISTING  IS  COMPLETE  FOR  MATL  '.A2,'-',A2. 

2/ , T2 , 'TO  CONTINUE  THE  PROGRAM  KEY  RETURN') 
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28 

I SET  1 * IHOLD1 

29 

ISET3* IHOLD2 

30 

return 

31 

1 

IHOLD1 = 1 SET2 

32 

IH0LD2* I SET  1 

33 

ISET1*IH0LD1 

34 

I SET2* 1H0LD2 

35 

RETURN 

36 

2 

READ ( 20 1 I SET  1 .200) ,SY1 

37 

RE AD ( 20 ' I SET  2 , 200 ) ,SY2 

38 

I F ( S Y 1 .GE.SY2)  RETURN 

39 

GO  TO  1 

40 

C 

42 

C 

43 

200 

FORMAT ( 61 X , F7 . 4 ) 

44 

C 

45 

END 

1 c 

2 SUBROUTINE  CALC ( NC , EMOD , XN , SY , S , E , EP , SS , EE , EMAX , NP ) 

3 

4 C 

5 C » SUBROUTINE  CALCULATES  STRESS-STRAIN  DATA  FROM  WHICH  CURVE  IS  PLOTTED 

6 C • PLASTIC  STRAIN  TERM  OF  R-0  EQUATION  USED:  EP* . 002* ( S/SY ) **XN 

7 C * EQUATION  SOLVED  FOR  S FOR  GIVEN  INCREMENTED  VALUES  OF  EP 

8 C * LAST  POINT  CALCULATED  IS  VALUE  OF  S FOR  EP*0.006  IN/IN 

9 C * PROPORT  L LIMIT  & 0.1%  YIELD  STRESS-STRAIN  & 0.2%  YIELD  STRAIN  CALCULATED 


10 

C 

12 

L” 

c 

13 

DIMENSION  S( 1 ) ,E( 1 ) ,EP( 1 ) ,SS(1 ),EE( 1 ) 

14 

EMOD* EMOD* 1 ■ E3 

15 

XP  = 1 ./XN 

16 

S(2)=SY»5.E-3**XP 

17 

E ( 2 ) * S ( 2 ) /EMOD 

18 

SS< 1 )=S<2) 

19 

EE( 1 ) = E ( 2 ) 

20 

SS(2)*SY*5.E-1**XP 

21 

EE ( 2 ) = 1 . E-3+SS ( 2 ) /EMOD 

22 

SS ( 3 ) * SY 

23 

EE ( 3 ) * 2 . E-3+SY/EM0D 

24 

DO  1 1*3. NC 

25 

S(I)*SY*(EP(I)/2.E-3)**XP 

26 

E(I)*EP( I )+S( I ) /EMOD 

27 

IF(£( I) .LT.EMAX)  GO  TO  1 

28 

GO  TO  2 

29 

i 

CONTINUE 

30 

NP*NC 

31 

RETURN 

32 

2 

NP*  I 

33 

RETURN 

34 

C 

35 

END 
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1 OASG.A  R.LOCATER. 

2 OUSE  1 0 , R-  LOCATER 

3 OASG.A  R.DATADANK. 

4 OUSE  20 i R*DAT ABANK 

5 OXQT  R.RE'RIEVE.TYP 


1 OASG.A  R« LOCATER. 

2 OUSE  1 0 . R*  LOCATER 

3 OASG.A  R*DAT ABANK . 

4 OUSE  20 , R*DAT ABANK 

5 OXQT  R.RETRIEVE.CRT 


1 OPREP  R-RETRIEVE. 

2 OMAP.I  LISS,R*RETRIEVE.TYP 

3 IN  R*RETRIEVE.PMAIN 

4 LIB  R*RETRI EVE . 

5 END 


1 OPREP  R*RET  RI EVE . 

2 OMAP.I  CIS.R-RETRIEVE.CRT 

3 IN  R-RETRIEVE.TMAIN 

4 IN  R*RETRIEVE.BUFFPK 

5 LIB  R*RETRIEVE. 

6 LIB  OAO*TEX . 

7 END 


LISTING  OF  R*FIXBANK 


L 


The  FIXBANK  file  consists  of  two  main  program  elements  and  four  control 
elements  as  follows: 


Element 


Function 


MAIN1 

MAIN2 

START 1 
START 2 
MAPI 

MAP2 


Print  tabular  listing  of  the  DATABANK  and  LOCATER 
file  contents 

Add  data  to  DATABANK  file,  update  LOCATER  file, 
and  print  updated  contents  of  both  data  files 

Control  routine  to  assign  DATABANK  and  LOCATER 
files  and  to  execute  MAIN1 

Control  routine  same  as  START1  except  executes 
MAIN2 

Control  routine  used  when  MAIN1  is  revised  to 
PREP  the  file,  to  form  the  absolute  (name: 
PRINT) , to  PACK  the  file,  and  to  print  a list 
of  file  elements 

Same  function  as  MAPI,  but  for  MAIN2  (absolute 
name:  UPDATE). 


Each  FIXBANK  program  is  a batch-process  execution  from  a UNIVAC  9300  hard- 
copy terminal.  The  input  to  MAIN2  consists  of  the  two  decks  of  output  cards 
punched  by  the  EVALRO  program  (see  Reference  1). 

Listings  of  the  two  main  programs  are  given  on  pages  56-59.  The  four  con- 
trol elements  are  listed,  in  the  order  given  above,  on  page  60. 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 


C 

c • file  r*fixbank  is  for  management  of  r*locater  and  r*databank  files 

C • THIS  IS  ELEMENT  MAIN1  (ABSOLUTE:  PRINT) 

C 

c 

C * MAIN1  PRINTS  CONTENTS  OF  R»  LOCATER  AND  R*DAT ABANK  ON  9300  TERMINAL 

C 

C 

C 

DIMENSION  LOC (19) 

C 

C 

C 

C * FILE  10  IS  THE  LOCATER  FILE 

C • FILE  20  IS  THE  DATABANK  FILE 

C 

DEFINE  FILE  1 0 ( 50 . 80 , E , J J ) 

DEFINE  FILE  20 ( 9999 , 80 , E . 1 1 ) 

C 

C 

C 

C * IDENTIFY  SET  NUMBER  OF  LAST  DATA  BANK  ENTRY  (LAST  SET*NI) 

C * LAST  SET  NUMBER  IS  UPPER  LIMIT  OF  DATA  BANK  PRINT-OUT  DO-LOOP 
C 

C * SECTION  1 OF  LOCATER  FILE  PRINTED 
C 

NI«0 

PRINT  203 
DO  20  J«1 ,7 

READ( 10' J, 102) , ( LOC( L) , L»1 , 19) 

PRINT  204.J.W, (LOC(L) ,L*1 ,19) 

DO  10  L*1 ,19 

I F ( LOC ( L ) . LE . N I ) GO  TO  10 
N I - LOC ( L ) 

10  CONTINUE 
20  CONTINUE 
C 

c 

C * CONTENT  OF  SECTION  2 OF  LOCATER  NOW  PRINTED 

c 

PRINT  205 
DO  30  d«1  ,7 
N.J+30 

READ(10'N,101),W,(LOC(L),L«1,19) 

PRINT  204,N,W,(LOC(L),L«1,19) 

30  CONTINUE 

c 

c 

C • PRINT-OUT  OF  DATA  BANK 
C 

NPAGE* 1 

PRINT  200, NPAGE 
DO  50  I»1 , N I 

RE AD (20' 1,100) ,SN,W,X,BK,NPG,NTY,NDA,KT,KE,SU,ELG,RA,RCH.EMOD. 
1XN, SY, SEE, NEXT, ISET 

PRINT  201 ,sn,w.x,bk,npg,nty,nda,kt,ke,su,elg,ra,rch,emod. 

1XN.SY, SEE, NEXT, ISET 
1F(M00< I .5) .NE.O)  GO  TO  50 
IF(MOD( I .45) .EQ.O)  GO  TO  40 
PRINT  202 
GO  TO  50 

40  IF(N.EQ.NI)  GO  TO  50 
NPAGE«NPAGE+1 
PRINT  200. NPAGE 
SO  CONTINUE 


S6 


1 


69  C 

70  C 

71  C 

72  100  FORMAT! A6.2A2, A4 , 12, II , 1 6 , 15 , 1 2 , F6 . 1 , 3F4 . 1 , F6 . 2 . 2F7 . 2 , F4 . 1 , 2 14 ) 

73  101  FORMAT! A2. 1914) 

74  102  F0RMAT(2X, 1914) 

75  103  FORMAT ( 72X , 2 1 4 ) 

76  C 

77  C 

78  C 

79  200  FORMAT! 1H1 ,T5, 'PAGE  NUMBER' ,13,/. 

80  1 T5.2X, 'SN' ,4X, 'WM' , 2X , ' XM ' ,2X, 'BOOK' , 2X , ' PG ' ,2X, 'T' . 

81  23X, 'DATE' ,4X, 'TEMP'  ,2X, 'ED' ,4X, ' ULT  ' ,3X,  'E  LNG ' , 3X , ' R-A ' ,3X,  ' R-C ' , 

B2  34X , ' EMOD' ,4X, ' EXPNT ' , 5X , ' SY ' , 5X , ' SEE ' , 2X , 'NEXT'  , 2X , ' I SET ' . / ) 

83  201  FORMAT (T5,A6,2X,2(A2,2X) , A4 , 2X , 1 2 , 2X , 1 1 , 2X , 16 , 2X , 15 , 2X , 12 . 2X , 

84  1 F6. 1 ,2X,3(F4. 1 ,2X) ,F6.2,2X,2( F7.2.2X) , F4. 1 ,2(2X, 14) ) 

85  202  FORMAT! 1H  ) 

86  203  F0RMAT(1H1 ,T10, 'CONTENT  OF  LOCATER  F I LE ' , , T 1 0 , 

87  1 ' NO ' ,2X, 'MT ' ,3X, 'ZZ' ,2X, ' 1A' ,2X, *2A' ,2X, '3A‘ , 2X , ' 4 A ' , 2X , ' 5A ' , 2X , 

68  2'6A'  ,2X,  '7A' ,2X, '8A' ,2X,  '9A' ,2X, ' IB'  ,2X, '28*  r2X, '38* ,2X, '4B' , 

! 89  32 X, '5B' ,2X, '68* ,2X, '78' ,2X, '8B' ,2X, '9B' ,/) 

90  204  FORMAT (T10, 1 2 , 2X , A2 , IX, 1914) 

91  205  FORMAT (/,T10,83( ' - ' ) ,/,T10, 1 11 * • NUMBER  OF  ENTRIES  FOR  EACH  MATERIAL 

92  1/ALLOY  ,/,T10,83( '-' ),/) 

93  C 

94  END 


1 C 

2 C • FILE  R » F I XBANK  IS  FOR  MANAGEMENT  OF  R« LOCATER  AND  R*DATABANK  FILES 

3 C • THIS  IS  ELEMENT  MAIN2  (ABSOLUTE:  UPDATE) 

4 C 

5 C 

6 C * MAIN2  UPDATES  CONTENTS  OF  R«  LOCATER  AND  R*DATABANK  FILES 

7 C * IT  ALSO  PRINTS  OUT  THE  CONTENTS  OF  EACH  OF  THE  UPOATEO  FILES 

8 C 

9 C * PROGRAM  INPUT:  CARDS  PUNCHED  BY  ANALYSIS  PROGRAM  R*EVALRO 
IOC*  TWO  DECKS  OF  CAROS:  ONE  FOR  DATABANK  AND  ONE  FOR  LOCATER  FILE 

1 1 C 

1 2 C 

13  C 

14  DIMENSION  LOC ( 1 9 ) , NLOC! 19), NHOLD! 10,19), NUM( 10,19) 

15  C 

16  C 

17  C 

18  C * FILE  10  IS  THE  LOCATER  FILE 

19  C * FILE  20  IS  THE  DATABANK  FILE 

20  C 

21  DEFINE  FILE  1 0 ( 50 , 80 , E , U J ) 

22  DEFINE  FILE  20 ( 9999 , 80 , E , 1 1 ) 

23  C 

24  

25  C 

26  C • UPDATE  OF  DATABANK 

27  C * EACH  LINE  OF  DATA  BANK  IS  A CARD  IMAGE  OF  A DATA  CARO  OF  EVALRO 

28  C * LINE  NUMBER  IS  SAME  AS  DATA  SET  NUMBER  (I*ISET)  AND  IS  IN  COLS  77-80 

29  C • SET  NUMBER  OF  LAST  PREVIOUS  ENTRY  FOR  SAME  MATL/ALLOY  IN  COLS  73-76 

30  C * LAST  PREVIOUS  ENTRY  SET  NO  IS  CALLED  NEXT  ( SET*NEXT ) 

31  C 

32  

33  C 

34  C • READ  EVALRO  DATA  CARDS  AND  WRITE  CARD  IMAGES  IN  DATA  BANK 

35  C * LAST  DATA  CARD  IS  A DUMMY  WITH  SET  NO.  OF  ZERO 

36  C 

37  1 READ! 5,100) ,SN,W, X , BK , NPG , NTY . NDA , KT.KE.SU, ELG.RA.RCH, EMOD, 

38  1XN, SY, SEE, NEXT, ISET 


57 


I 


39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 


I F ( I SET . EQ - 0 ) GO  TO  2 
I-ISET 

WRITE ( 20 ' I i 100) , SN , WiX. BK , NPG, NTY , NDA , KT , KE , SU . ELG , RA , RCH.EMOD. 

1XN.SY, SEE, NEXT, ISET 
GO  ' 

C 

C 

C * READ  ENTIRE  CONTENT  OF  DATA  BANK  AND  PRINT 
C 

2 NI«I 
NPAGE* 1 

PRINT  200, NPAGE 
DO  4 1*1 ,NI 

READ (20' 1,100) , SN, W.X.BK, NPG, NT Y, NDA, KT.KE.SU.ELG.RA, RCH.EMOD. 

1XN.SY, SEE, NEXT, ISET 

PRINT  201 , SN, W, X, BK, NPG, NTY, NDA, KT.KE.SU.ELG.RA, RCH.EMOD. 

1XN.SY, SEE, NEXT, ISET 
I F ( MOD ( I , 5 ) • NE • 0 ) GO  TO  4 
I F ( MOD ( I , 45 ) . EO . 0 ) GO  TO  3 
PRINT  202 
GO  TO  4 

3 NPAGE*NPAG£+1 
PRINT  200, NPAGE 

4 CONTINUE 

c 

c 

c 

C * UPDATE  OF  LOCATER  FILE 
C 

C * LOCATER  FILE  HAS  TWO  SECTIONS:  1.  KEYSET  LOCATION  AND  2.  NUMBER  OF  ENTRIES 
C * SECTION  1 (LINES  1-30):  SET  NO.  OF  LAST  DATABANK  ENTRY  FOR  EACH  MATL/ALLOY 
C * EACH  MATER I AL  ON  A SEPARATE  LINE  (LINE  NO.*J) 

C * COLS  1 & 2 MATERIAL  CODE  FOLLOWED  BY  19  FOUR-COLUMN  GROUPS 
C * EACH  GROUP  (L=GROUP  NO.)  REFERS  TO  SPECIFIC  ALLOY  • ENTRY  IS  LAST  SET  NO. 

C * COLS  79-80  BLANK  AND  RESERVED  FOR  LATER  USE  AS  CONTINUATION  LOCATIONS 
C * SECTION  2 (LINES  31-50):  NUMBER  OF  DATA  ENTRIES  FOR  EACH  MATERIAL/ALLOY 
C * MATL/ALLOY  OF  LINE  U POSITION  L SAME  AS  FOR  LINE  J+30  POSITION  L 
C 

C 

C 

C * UPDATE  OF  SECTION  1 OF  LOCATER 

C * READ  EVALRO  LOCATER  CAROS  * READ  AND  HOLD  LOCATER  VALUES  CURRENTLY  IN  FILE 
C 

PRINT  203 
J*0 

10  READ(5, 101 ,END*20) , W , ( LOC ( L) , L* 1 , 1 9 ) 

U*  J+1 

READ( 10' J, 102) , ( NHOLD( J , L) , L* 1 ,19) 

C 

C 

C ENTER  NEW  LOCATER  VALUES  IN  FILE  AND  PRINT  FILE  CONTENT 
C 

WRITE! 10' U, 101 ),W,( LOC (L),L*1 ,19) 

PRINT  204.J.W, (LOC(L) ,L=1 ,19) 

GO  TO  10 
C 



c 

C * UPDATE  OF  SECTION  2 OF  LOCATER 

C * NUMBER  OF  ENTRIES  OF  EACH  MATL/ALLOY  JUST  ENTERED  IN  DATABANK  COUNTED 
C * THIS  NUMBER  ADDED  TO  CURRENT  TOTALS  IN  SECTION  2 OF  LOCATER 
C • THEN  NEW  NUMBER  OF  ENTRIES  WRITTEN  IN  SECTION  2 OF  LOCATER  TO  UPDATE  IT 
C 

20  NJ* J 

DO  24  J*1  , N J 
N* J+30 


) 


I 
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107 
1 OS 

109 

110 
1 1 1 
112 

113 

114 

115 

116 
117 
1 18 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 
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140 

141 

142 

143 

144 
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147 
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151 

152 

153 
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155 

156 

157 

158 

159 

160 
161 
162 


READ(10'J,101).W,(LOC(L),L«1,19) 

00  22  1*1 ,19 
NUM( J , L ) *0 
NI=LOC(L) 

IF(NI.EQ.O)  GO  TO  22 

21  READ(20'NI, 103), NEXT, ISET 

I F ( I SET . EQ . NHOLD ( U , L) ) GO  TO  22 
NUM(  J,L)*NUM(J,L)-M 

IF(NEXT.EQ.O.OR.NEXT.EQ.NHOLD(J,L))  GO  TO  22 
N I =NEXT 
GO  TO  21 

22  CONTINUE 

READ( 10'N, 102) , (NLOC(L) , L*1 ,19) 

DO  23  L=1 , 19 

NLOCI L)=NLOC( L)+NUM(d.L) 

23  CONTINUE 

WRITE! 10‘N, 101 ) ,W, (NLOC( L) ,1*1 .19) 

24  CONTINUE 
C 

c 

C * CONTENT  OF  SECTION  2 OF  LOCATER  NOW  PRINTED 
C 

PRINT  205 
DO  30  d=1 ,NJ 
N=J+30 

RE AD ( IO'N.101 ) ,W, (NLOC( L) , L*1 , 19) 

PRINT  204 , N , W , ( NLOC ( L ) , L* 1 .19) 

PUNCH  101 ,W. (NLOC(L) ,L*1 ,19) 

30  CONTINUE 
C 

c 

100  FORMAT ( A6 ,2A2,A4,I2,I1, 16 , 15 , 1 2 , F6 . 1 ,3F4. 1 , F6 . 2 . 2F7 . 2 , F4 . 1 , 2 1 4 ) 

101  FORMAT! A2, 1914) 

102  FORMAT ( 2X , 1 9 1 4 ) 

103  F0RMAT(72X,2I4) 

C 



c 

200  F0RMAT(1H1 ,T5,  'PAGE  NUMBER 1 , 1 3 , / , 

1 T5.2X, 'SN' ,4X, 1 WM ' , 2x , ' XM 1 , 2X , 'BOOK' , 2X , ' PG ' , 2X , ' T • , 

23X. ' DATE' ,4X. ' TEMP* . 2X , ' ED ' . 4X , • ULT • , 3X , ' E LNG ' , 3X , ' R-A 1 , 3X . 'R-C* . 
34X, 1 EMOD'  ,4X, ' EXPNT ' , 5X , ' SY 1 , 5X . 'SEE'  ,2X, 'NEXT'  , 2X . ' I SET ' . / ) 

201  FORMAT (T5. A6.2X.2! A2.2X) , A4.2X, 1 2 , 2X , II ,2X, I6.2X, I5.2X.I2.2X, 

1F6.1  ,2X,3(F4.1 ,2X),F6.2,2X,2(F7.2,2X),F4.1 , 2 ( 2X , 1 4 ) ) 

202  FORMAT! 1H  ) 

203  FORMAT ( 1K1 ,T10, 'CONTENT  OF  LOCATER  F I LE ' , // , T 1 0 . 

1 'NO'  ,2X,  'MT ' ,3X,  '22'  ,2X, ' 1 A' ,2X, '2A' ,2X,  ' 3A1 , 2X , ' 4A ' ,2X, 1 5A ' ,2X. 

2 ' 6 A ' , 2X . ' 7A ' , 2X , 1 8 A 1 , 2X , 1 9 A ' , 2X , ' 1 8 ' , 2X , ' 2B* , 2X , ' 3B ' , 2X , ' 4B 1 • 

32 X, ' 5B ' ,2X, ' 6B ' , 2X , ' 7B ' , 2X , 'SB' , 2X , ' 9B 1 ,/) 

204  F0RMAT(T10. I2.2X.A2, IX. 1914) 

205  FORMAT!/, T10, 83( ’-') ,/.T10, NUMBER  OF  ENTRIES  FOR  EACH  MATERIAL 
1/ALLOY  .*•,/,  T10.B3C-')./) 

C 

END 
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1 PASG.A  R»  DAT  A BANK 

2 PUSE  20,R*DATABANK 

3 9ASG , A R*LOCATER 

4 PUSE  1 0 , R»  LOCATER 

5 PXQT  R«FIXBANK. PRINT 


1 PASG.A  R*DATABANK 

2 PUSE  20,R*DATABANK 

3 PASG.A  R. LOCATER 

4 PUSE  1 0 , R*  LOCATER 

5 PXQT  R*FIXBANK. UPDATE 


1 PPREP  R* F I XBANK . 

2 PMAP.I  CISSY, R»FIXBANK. PRINT 

3 IN  R*FIXBANK.MAIN1 

4 END 

5 PPACK  R* F I XBANK . 

6 PPRT.TL  R * F I XBANK . 


1 PPREP  R* F I XBANK . 

2 ©MAP, I CIS, R*FIXBANK, UPDATE 

3 IN  R*FIXBANK.MAIN2 

4 END 

5 PPACK  R*FIXBANK. 

6 PPRT.TL  R* " I XBANK . 
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